狄利克雷卷积和莫比乌斯反演

积性函数

定义

一个数论函数f,\forall (n,m)=1,有f(nm)=f(n)×f(m),那么称f为积性函数。
一个数论函数f,对于\forall n,m,有 f(nm)=f(n)×f(m),那么称f为完全积性函数
其中数论函数的含义是

在数论上,算术函数(或称数论函数)指定义域为正整数、陪域为复数的函数,每个算术函数都可视为复数的序列。

主要概念就是定义域是正整数。

性质

性质:若f(x),g(x)是积性函数那么下面的函数都是积性函数
h(x) = f(x^p)
h(x) = f^p(x)
h(x) = f(x)g(x)
h(n) = \sum_{d\mid n}f(d)g(\frac n d)
证明方法可以将x拆分为两个数a,b满足a*b=x,(a,b)=1。然后验证h(x) = h(a)*h(b)就可以了。最后一个同样可以拆分。
h(ab) = \sum_{d\mid ab}f(d)g(\frac {ab} d)
h(b)*h(b) = (\sum_{d_1\mid b}f(d_1)g(\frac {b} {d_1}))*(\sum_{d_2\mid a}f(d_2)g(\frac {a} {d_2}))
h(b)*h(b) = \sum_{d\mid ab}f(d)g(\frac {ab} d)
可以这样想d是a*b的正约数,d_1是a的正约数,d_2是b的正约数。(a,b)=1,所以d_1和d_2不包含相同的质因数。a*b的正约数是d=d_1*d_2

举例

常见的积性数论函数
单位函数\epsilon(n) = [n=1]
常数函数1(n)=1
恒等函数id(n) = n
莫比乌斯函数\mu(n)= \begin{cases} 1 & n = 1\\ (-1)^k & n = p_1*p_2*...*p_k \forall p_i \neq p_j\\ 0 & others \end{cases}
欧拉函数\varphi(n) = \sum_{i=1}^n [(i,n)=1]
除数函数:\sigma_{k}(n) = \sum _{d\mid n} d^k当k取0时,记作d(n)当k取1记作\sigma (n)
其中莫比乌斯函数是当n=1时为1,对于任意一个质数p,如果p|n,但是p^2\nmid n也就是说对n进行质因数分解后,不含有相同的质因数的时候莫比乌斯函数为(-1)^k,k是分解后质因数的个数,其他情况为0。除数函数k为0的时候求的是正约数的数量,k为1的时候求的是正约数的和。单位函数[n=1]的含义是n=1的时候为1,否则为0,就是布尔运算。欧拉函数的(i,n)=1是gcd(i,n)=1。

狄利克雷卷积

定义

两个数论函数f(x),g(x)的狄利克雷卷积是
h(x) = (f*g)(x) = \sum _{d\mid n}f(d)*f(n/d) = \sum _{a*b=n}f(a)*g(b)

性质

交换律:f*g= g*f
结合律:(f*g)*h = g*(f*h)
分配律:f*(g+h) = f*g+f*h
单位元:f*\epsilon = f
逆元:对于每个f(1)\neq1的数论函数都有f*g=\epsilon

交换律证明: (f*g)(x) = \sum _{a*b=n}f(a)*g(b) = \sum _{a*b=n}g(a)*f(b)=(g*f)(x)
结合律证明:((f*g)*h)(n) = \sum _{d*c=n}(f*g)(d)*g(c) = \sum _{d*c=n}(\sum _{a*b=d} f(a)*g(b))*g(c)
= \sum_{a*b*c=n}f(a)*g(b)*h(c)
剩下的写一下公式就可以推导出来。

卷积关系

\epsilon = \mu*1
d=1*1
\sigma = id*1
\varphi = \mu * id


\epsilon = \mu*1证明
f = (\mu * 1) (n)=\sum_{d\mid n} \mu(d),根据积性函数的性质得到f也是积性函数,上面已经证明过了,在证明一次因为\forall (a,b)=1,a*b=n,f(n) = \sum_{d\mid a*b} \mu(d) = ( \sum_{d\mid a} \mu(d))*( \sum_{d\mid b} \mu(d))=f(a)*f(b)。那么对n进行质因数分解n=p_1^{c_1}*p_2^{c_2}*...*p_m^{c_m},f(n) = f(p_1^{c_1})*f(p_2^{c_2})*...*f(p_m^{c_m}),单独考虑每一个f(p^c) = \sum _{d|p^c}\mu (d) = \sum_{i=0}^c \mu(p^i)
\sum_{i=0}^c \mu(p^i)单独来看这个公式,\mu(1) = 1,\mu(p)=-1,\mu(p^2) = 0,....,所以\sum_{i=0}^c \mu(p^i) = 0带回f(n)中得到f(n)=0,证明完毕。
得到\epsilon = \mu*1 = \sum_{d\mid n} \mu(d)


d=1*1证明
(1*1)(n) = \sum_{d\mid n}1(d)*1(n/d) = \sum_{d\mid n} 1 = d(n)


\sigma = id*1证明
(id*1)(n) = \sum_{d\mid n} id(d)*1(n/d) = \sum_{d\mid n}d = \sigma(n)


\varphi = \mu * id
1*\varphi =1*\mu * id = id
只需要证明1*\varphi = id即可。
我们知道f(n) = (1*\varphi)(n) = (\varphi * 1)(n) = \sum_{d\mid n}\varphi(d)根据积性函数的性质可知,\sum_{d\mid n}\varphi(d)也是积性函数。那么对n进行质因数分解n=p_1^{c_1}*p_2^{c_2}*...*p_m^{c_m},f(n) = f(p_1^{c_1})*f(p_2^{c_2})*...*f(p_m^{c_m}),单独考虑每一个f(p^c) = \sum _{d|p^c}\varphi (d) = \sum_{i=0}^c \varphi(p^i) = \varphi(p^0)+\varphi(p^1)+...+\varphi(p^c),单独考虑\varphi(p^0)=1,\varphi(p^1)=p-1,\varphi(p^2) = p*(p-1),\varphi(p^c) = p^{c-1}*(p-1)所以\sum_{i=0}^c \varphi(p^i)=p^c,f(p^c)=p^c带回f(n)中得到f(n)=n,(1*\varphi)(n) = id(n)证明完毕。

莫比乌斯反演

莫比乌斯函数
\mu(n)= \begin{cases} 1 & n = 1\\ (-1)^k & n = p_1*p_2*...*p_k \forall p_i \neq p_j\\ 0 & others \end{cases}
第二情况是n可以被分解为k个质数,并且这k个质数两两不同。
显然莫比乌斯函数满足p\nmid n,p\in prime的时候\mu(n*p) = -\mu(n)成立。
p\mid n,p\in prime的时候\mu(n*p) = 0显然我们可以用线筛O(n)的时机求1~n,
之间的莫比乌斯函数。绝大部分积性函数都可以O(n)的使用线筛时间内求出。
莫比乌斯推论,我们有\mu * 1 = \epsilon
也就是(\mu * 1)(n) = \epsilon(n)但是我们可以对他进行转换
(\mu * 1)(gcd(a,b)) = \epsilon(gcd(a,b)) = [gcd(a,b) = 1] = \sum_{d\mid gcd(a,b) }\mu(d)
对于某些a,b我们可以求出互质的个数。

莫比乌斯反演定理

假设有两个数论函数f,F满足
F(n) = \sum_{d\mid n}f(d)
那么有
f(n) = \sum_{d\mid n}\mu(d)F(\frac n d)


狄利克雷卷积证明:
F = f*1
F *\mu= f*1*\mu
f = \mu * F

其他证明方法:
\sum_{d\mid n}\mu(d)F(\frac n d)=
F(n) = \sum_{d\mid n}f(d)带入
\sum_{d\mid n}\mu(d)(\sum_{p\mid \frac n d}f(p))=
因为能对结果产生贡献必须满足(d*p)|n也就是\mu(d)*f(p),因此换一种形式也是正确的。
\sum_{p\mid n}f(p)(\sum_{d\mid \frac n p}\mu(d))=
根据上面验证的\sum_{d\mid n}\mu(d) = \epsilon(n)可以得知,只有n=1的时候对最后的结果有贡献,也就是p=n
f(n)


莫比乌斯反演的另一种形式
F(n)=\sum_{n|d}f(d)
f(n)=\sum_{n|d}\mu(\frac{d}{n})F(d)

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 215,634评论 6 497
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,951评论 3 391
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 161,427评论 0 351
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,770评论 1 290
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,835评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,799评论 1 294
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,768评论 3 416
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,544评论 0 271
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,979评论 1 308
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,271评论 2 331
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,427评论 1 345
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,121评论 5 340
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,756评论 3 324
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,375评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,579评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,410评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,315评论 2 352