莫比乌斯反演学习笔记

莫比乌斯函数

定义莫比乌斯函数:

μ ( 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=p1p2p3pkpi为互异质数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} dnμ(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=p1a1p2a2pkak

    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*} dnμ(d)==Ck0Ck1+Ck2Ck3++(1)kCkki=0k(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=0nCnixiyni,nN
    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=0nCni(1)i1ni=0,nN

    证毕.

    此处可以将该性质化为下面这种形式(这才是真正意义上的莫比乌斯函数),以便于下文利用:
    ∑ d ∣ n μ ( d ) = [ n = 1 ] \displaystyle \sum_{d|n}\mu(d)=\left [ n=1 \right ] dnμ(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 μ(xy)=μ(x)μ(y),gcd(x,y)=1

积性函数的性质:

  1. μ ( 1 ) = 1 \mu(1)=1 μ(1)=1.
  2. 积性函数的前缀和也是积性函数.

性质三

对于任意正整数 n n n 有:
∑ d ∣ n μ ( d ) d = ϕ ( n ) n \displaystyle \sum_{d|n}\frac{\mu(d)}{d}=\frac{\phi(n)}{n} dndμ(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)=dnf(d)f(n)=dnμ(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)==dnμ(d)dnndndμ(d)

最后移项得:
∑ d ∣ n μ ( d ) d = ϕ ( n ) n \displaystyle \sum_{d|n}\frac{\mu(d)}{d}=\frac{\phi(n)}{n} dndμ(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 dnϕ(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}) dnϕ(d)=dnnϕ(d)=dnϕ(dn)

对于任意的 d ∣ n d|n dn 1 ≤ a ≤ n 1\leq a \leq n 1an,设 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 dnϕ(dn)=n

证毕.

莫比乌斯反演定理

形式

对于任意正整数 n n n,若:
F ( n ) = ∑ d ∣ n f ( d ) F(n)=\displaystyle \sum_{d|n}f(d) F(n)=dnf(d)

则:
f ( n ) = ∑ d ∣ n F ( d ) μ ( n d ) f(n)=\displaystyle \sum_{d|n}F(d)\mu(\frac{n}{d}) f(n)=dnF(d)μ(dn)

证明方法:狄利克雷卷积

在这里只给出简单证明,详细的可以查看CJX的博客了解。
我们知道狄利克雷卷积的形式:
∑ d ∣ n f ( d ) ε ( n d ) \displaystyle \sum_{d|n}f(d)\varepsilon(\frac{n}{d}) dnf(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) dnf(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)+dnd=nf(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 fg=ε,则称 f f f g g g 互为逆元,可以把 g g g 写成 f − 1 f^{-1} f1.

那么就有 f ⊕ g = ε ⇒ f = ε ⊕ g − 1 f \oplus g=\varepsilon \Rightarrow f = \varepsilon \oplus g^{-1} fg=εf=εg1.

证明

考虑运用上述性质一: ∑ d ∣ n μ ( d ) = [ n = 1 ] \displaystyle \sum_{d|n}\mu(d)=\left [ n=1 \right ] dnμ(d)=[n=1].

我们发现它的形式和狄利克雷卷积的形式很像,我们把它写成狄利克雷卷积的形式,就有:
∑ d ∣ n μ ( d ) I 1 ( n d ) \displaystyle \sum_{d|n}\mu(d)I_1(\frac{n}{d}) dnμ(d)I1(dn)

这里为了符合狄利克雷卷积的形式,我们定义了一个函数 I 1 I_1 I1,把它写在式子的右边。

我们惊奇的发现莫比乌斯函数 ∑ d ∣ n μ ( d ) \displaystyle \sum_{d|n}\mu(d) dnμ(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)=dnf(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)=dnf(d)I1(dn)f(n)=dnF(d)I11(dn)f(n)=dnF(d)μ(dn)

得证.

应用情境

对于一些函数 f ( n ) f(n) f(n),如果我们很难直接求出它的值,而容易求出倍数和或约数和 F ( n ) F(n) F(n),那么我们可以通过莫比乌斯反演来求得 f ( n ) f(n) f(n) 的值。

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值