莫比乌斯函数
定义莫比乌斯函数:
μ ( n ) = { 1 , n = 1 ( − 1 ) k , n = p 1 p 2 p 3 … p k 且 p i 为互异质数 0 , 其他情况 \mu(n)=\begin{cases} 1,n=1\\(-1)^k,n=p_1p_2p_3\dots p_k且p_i为互异质数\\ 0,其他情况\end{cases} μ(n)=⎩ ⎨ ⎧1,n=1(−1)k,n=p1p2p3…pk且pi为互异质数0,其他情况
性质一
对于任意正整数 n n n 有:
∑ d ∣ n μ ( d ) = { 1 , n = 1 0 , n > 1 \sum\limits_{d|n}\mu(d)=\begin{cases}1,n=1\\0,n>1 \end{cases} d∣n∑μ(d)={1,n=10,n>1
证明:
-
当 n = 1 n=1 n=1 时显然正确
-
当 n ≠ 1 n \ne 1 n=1 时,可以将 n n n 分解得到 n = p 1 a 1 p 2 a 2 … p k a k n=p_1^{a_1}p_2^{a_2}\dots p_k^{a_k} n=p1a1p2a2…pkak
在 n n n 的所有因子中, 值不为零的只有所有质因子次数都为 1 1 1 的因子,那么就有:
∑ d ∣ n μ ( d ) = C k 0 − C k 1 + C k 2 − C k 3 + ⋯ + ( − 1 ) k C k k = ∑ i = 0 k ( − 1 ) i C k i \begin{align*} \displaystyle \sum_{d|n} \mu(d) = &C_{k}^{0} -C_k^1+C_k^2-C_k^3+\dots+(-1)^kC_k^k \\ = & \displaystyle \sum_{i=0}^{k} (-1)^iC_k^i \end{align*} d∣n∑μ(d)==Ck0−Ck1+Ck2−Ck3+⋯+(−1)kCkki=0∑k(−1)iCki
容易发现其形式类似于二项式定理,可以由:
( x + y ) n = ∑ i = 0 n C n i x i y n − i , n ∈ N ∗ (x+y)^n=\displaystyle \sum_{i=0}^n C_n^ix^iy^{n-i},n \in N^* (x+y)n=i=0∑nCnixiyn−i,n∈N∗
令 x = − 1 , y = 1 x=-1,y=1 x=−1,y=1 得到:
( − 1 + 1 ) n = ∑ i = 0 n C n i ( − 1 ) i 1 n − i = 0 , n ∈ N ∗ (-1+1)^n=\displaystyle \sum_{i=0}^n C_n^i(-1)^i1^{n-i}=0,n \in N^* (−1+1)n=i=0∑nCni(−1)i1n−i=0,n∈N∗证毕.
此处可以将该性质化为下面这种形式(这才是真正意义上的莫比乌斯函数),以便于下文利用:
∑ d ∣ n μ ( d ) = [ n = 1 ] \displaystyle \sum_{d|n}\mu(d)=\left [ n=1 \right ] d∣n∑μ(d)=[n=1]其中等号右边的的方括号是艾弗森括号,若括号内的表达式成立则值为 1 1 1,否则为 0 0 0。
性质二
莫比乌斯函数是不完全积性函数,即:
μ
(
x
∗
y
)
=
μ
(
x
)
∗
μ
(
y
)
,
g
c
d
(
x
,
y
)
=
1
\mu(x*y)=\mu(x)*\mu(y),gcd(x,y)=1
μ(x∗y)=μ(x)∗μ(y),gcd(x,y)=1
积性函数的性质:
- μ ( 1 ) = 1 \mu(1)=1 μ(1)=1.
- 积性函数的前缀和也是积性函数.
性质三
对于任意正整数
n
n
n 有:
∑
d
∣
n
μ
(
d
)
d
=
ϕ
(
n
)
n
\displaystyle \sum_{d|n}\frac{\mu(d)}{d}=\frac{\phi(n)}{n}
d∣n∑dμ(d)=nϕ(n)
通过下面的莫比乌斯反演证明,可以先看下文,因为这个结论不重要 ,由:
F
(
n
)
=
∑
d
∣
n
f
(
d
)
f
(
n
)
=
∑
d
∣
n
μ
(
d
)
F
(
n
d
)
F(n)= \displaystyle \sum_{d|n}f(d)\\ f(n)= \displaystyle \sum_{d|n}\mu(d)F(\frac{n}{d})
F(n)=d∣n∑f(d)f(n)=d∣n∑μ(d)F(dn)
令
f
(
n
)
=
ϕ
(
n
)
,
F
(
n
)
=
n
f(n)=\phi(n),F(n)=n
f(n)=ϕ(n),F(n)=n,则:
ϕ
(
n
)
=
∑
d
∣
n
μ
(
d
)
∗
n
d
=
n
∗
∑
d
∣
n
μ
(
d
)
d
\begin{align*} \phi(n) = &\displaystyle \sum_{d|n}\mu(d)*\frac{n}{d}\\ = &n*\displaystyle \sum_{d|n}\frac{\mu(d)}{d}\\ \end{align*}
ϕ(n)==d∣n∑μ(d)∗dnn∗d∣n∑dμ(d)
最后移项得:
∑
d
∣
n
μ
(
d
)
d
=
ϕ
(
n
)
n
\displaystyle \sum_{d|n}\frac{\mu(d)}{d}=\frac{\phi(n)}{n}
d∣n∑dμ(d)=nϕ(n)
证毕.
至于为什么中间过程可以设
f
(
n
)
=
ϕ
(
n
)
,
F
(
n
)
=
n
f(n)=\phi(n),F(n)=n
f(n)=ϕ(n),F(n)=n,则需要证明对于任意正整数
n
n
n 有:
∑
d
∣
n
ϕ
(
d
)
=
n
\displaystyle \sum_{d|n}\phi(d)=n
d∣n∑ϕ(d)=n
证明如下:
容易得到
∑
d
∣
n
ϕ
(
d
)
=
∑
n
d
∣
n
ϕ
(
d
)
=
∑
d
∣
n
ϕ
(
n
d
)
\displaystyle \sum_{d|n}\phi(d)=\displaystyle \sum_{\frac{n}{d}\mid n}\phi(d)=\displaystyle \sum_{d|n}\phi(\frac{n}{d})
d∣n∑ϕ(d)=dn∣n∑ϕ(d)=d∣n∑ϕ(dn)
对于任意的 d ∣ n d|n d∣n, 1 ≤ a ≤ n 1\leq a \leq n 1≤a≤n,设 g c d ( a , n ) = d gcd(a,n)=d gcd(a,n)=d,则 g c d ( a d , n d ) = 1 gcd(\frac{a}{d},\frac{n}{d})=1 gcd(da,dn)=1.
一方面数 a a a 有 n n n 个,另一方面按照 d = g c d ( a , n ) d=gcd(a,n) d=gcd(a,n) 分类计数,满足 g c d ( a , n ) = d gcd(a,n)=d gcd(a,n)=d 的数 a a a 有 ϕ ( n d ) \phi(\frac{n}{d}) ϕ(dn) 种取法.
因此有 ∑ d ∣ n ϕ ( n d ) = n \displaystyle \sum_{d|n}\phi(\frac{n}{d})=n d∣n∑ϕ(dn)=n
证毕.
莫比乌斯反演定理
形式
对于任意正整数
n
n
n,若:
F
(
n
)
=
∑
d
∣
n
f
(
d
)
F(n)=\displaystyle \sum_{d|n}f(d)
F(n)=d∣n∑f(d)
则:
f
(
n
)
=
∑
d
∣
n
F
(
d
)
μ
(
n
d
)
f(n)=\displaystyle \sum_{d|n}F(d)\mu(\frac{n}{d})
f(n)=d∣n∑F(d)μ(dn)
证明方法:狄利克雷卷积
在这里只给出简单证明,详细的可以查看CJX的博客了解。
我们知道狄利克雷卷积的形式:
∑
d
∣
n
f
(
d
)
ε
(
n
d
)
\displaystyle \sum_{d|n}f(d)\varepsilon(\frac{n}{d})
d∣n∑f(d)ε(dn)
单位元
什么是单位元?对于一种运算
⊕
\oplus
⊕ ,若存在函数
ε
\varepsilon
ε 使得:
f
⊕
ε
=
g
⊕
ε
=
f
f \oplus \varepsilon=g \oplus \varepsilon=f
f⊕ε=g⊕ε=f
那我们就称函数 ε \varepsilon ε 为 f f f 的单位元。
例如加法运算的单位元为
0
0
0,因为
0
0
0 加上任何数都等于原数。
而乘法运算的单位元为
1
1
1,因为
1
1
1 乘上任何数都等于原数。
那么狄利克雷卷积的单位元
ε
\varepsilon
ε 是什么? 看以下推导:
要使
∑
d
∣
n
f
(
d
)
ε
(
n
d
)
=
f
(
n
)
\displaystyle \sum_{d|n}f(d)\varepsilon(\frac{n}{d})=f(n)
d∣n∑f(d)ε(dn)=f(n)
将
n
n
n 这一项提取出来,则
f
(
n
)
ε
(
1
)
+
∑
d
∣
n
d
≠
n
f
(
d
)
ε
(
n
d
)
=
f
(
n
)
f(n)\varepsilon(1)+\displaystyle \sum_{\substack{d|n\\d \ne n}}f(d)\varepsilon(\frac{n}{d})=f(n)
f(n)ε(1)+d∣nd=n∑f(d)ε(dn)=f(n)
显然对于狄利克雷卷积的单位元函数 ε ( n ) \varepsilon(n) ε(n),当 n = 1 n=1 n=1 时, ε ( n ) = 1 \varepsilon(n)=1 ε(n)=1;当 n ≠ 1 n \ne 1 n=1 时, ε ( n ) = 0 \varepsilon(n)=0 ε(n)=0,才能使得上式成立。
因此我们可以得到狄利克雷卷积的单位元函数是 ε ( n ) = [ n = 1 ] \varepsilon(n)=[n=1] ε(n)=[n=1].
逆元
对于单位元函数 ε \varepsilon ε,如果有 f ⊕ g = ε f \oplus g=\varepsilon f⊕g=ε,则称 f f f 和 g g g 互为逆元,可以把 g g g 写成 f − 1 f^{-1} f−1.
那么就有 f ⊕ g = ε ⇒ f = ε ⊕ g − 1 f \oplus g=\varepsilon \Rightarrow f = \varepsilon \oplus g^{-1} f⊕g=ε⇒f=ε⊕g−1.
证明
考虑运用上述性质一: ∑ d ∣ n μ ( d ) = [ n = 1 ] \displaystyle \sum_{d|n}\mu(d)=\left [ n=1 \right ] d∣n∑μ(d)=[n=1].
我们发现它的形式和狄利克雷卷积的形式很像,我们把它写成狄利克雷卷积的形式,就有:
∑
d
∣
n
μ
(
d
)
I
1
(
n
d
)
\displaystyle \sum_{d|n}\mu(d)I_1(\frac{n}{d})
d∣n∑μ(d)I1(dn)
这里为了符合狄利克雷卷积的形式,我们定义了一个函数 I 1 I_1 I1,把它写在式子的右边。
我们惊奇的发现莫比乌斯函数 ∑ d ∣ n μ ( d ) \displaystyle \sum_{d|n}\mu(d) d∣n∑μ(d) 其实就是狄利克雷卷积的单位元,那么什么样的函数可以与单位元互为逆元呢?
显然是常函数 I 1 ( x ) = 1 I_1(x)=1 I1(x)=1.
给
F
(
n
)
=
∑
d
∣
n
f
(
d
)
F(n)=\displaystyle \sum_{d|n}f(d)
F(n)=d∣n∑f(d) 乘上一个常函数
I
1
I_1
I1(对
F
(
n
)
F(n)
F(n) 没有影响),那么就有:
F
(
n
)
=
∑
d
∣
n
f
(
d
)
∗
I
1
(
n
d
)
⇒
f
(
n
)
=
∑
d
∣
n
F
(
d
)
∗
I
1
−
1
(
n
d
)
⇒
f
(
n
)
=
∑
d
∣
n
F
(
d
)
μ
(
n
d
)
F(n)=\displaystyle \sum_{d|n}f(d)*I_1(\frac{n}{d}) \\ \Rightarrow f(n)=\displaystyle \sum_{d|n}F(d)*I_1^{-1}(\frac{n}{d}) \\ \Rightarrow f(n)=\displaystyle \sum_{d|n}F(d)\mu(\frac{n}{d})
F(n)=d∣n∑f(d)∗I1(dn)⇒f(n)=d∣n∑F(d)∗I1−1(dn)⇒f(n)=d∣n∑F(d)μ(dn)
得证.
应用情境
对于一些函数 f ( n ) f(n) f(n),如果我们很难直接求出它的值,而容易求出倍数和或约数和 F ( n ) F(n) F(n),那么我们可以通过莫比乌斯反演来求得 f ( n ) f(n) f(n) 的值。