从零开始的莫比乌斯反演(函数)[详细推导]

也许更好的阅读体验

---

$前置技能$

• 学会莫比乌斯函数必须要先知道狄利克雷函数
• 以及什么是逆元（一本正经胡说八道）

  $狄利克雷卷积$

  两个函数 $f,g$
  $\begin{array}{r}f\ast g(n)=\sum _{d|n}^{n}f(d)g(\frac{n}{d})\end{array}$
  • 性质
    交换律$f\ast g=g\ast f$
    结合律$f\ast g\ast h=f\ast (g\ast h)$
    加法分配率$f\ast (g+h)=f\ast g+f\ast h$

  $几个定理$

  • ①　定义单位函数$I(n)=[n=1]$，即只有当ｎ等于１时，$I(1)=1$，其余情况都为０
  • ②　$I\ast f=f$
  • ③　已知函数$f$，定义函数$f$的逆$f^{-1}$，满足$f\ast f^{-1}=I$，且积性函数的逆也是积性函数
  • ④　由③可得$f^{-1}(1)=\frac{1}{f(1)}$
  • ⑤　令$\xi (n)=1$

---

$莫比乌斯函数$

• 莫比乌斯函数 $\mu ={\xi }^{-1}$
  $n=p_1^{a_1}\cdot p_2^{a_2}\cdot \cdots p_k^{a_k}$
  $\mu (n)={\xi }^{-1}(n)=\{\begin{array}{ccc}1& n=1& \\ (-1{)}^k& a1=a2=\cdots =ak=1& \\ 0& otherwise& \end{array}$

---

$如何推导$

• $如何求逆$

  ③等价于
  
  $\begin{array}{r}\forall n:\sum _{d|n}f(d)\cdot f^{-1}(\frac{n}{d})=[n=1]\end{array}$
  
  $\begin{array}{r}\sum _{d|n}f(d)\cdot f^{-1}(\frac{n}{d})=I(n)(n≠1)=0\end{array}$
  
  $f^{-1}(n)\cdot f(1)=-{\sum }_{d|n,d≠1}f(d)\cdot f^{-1}(\frac{n}{d})$
  
  $\begin{array}{r}{f}^{-1}(n)=\frac{-{\sum }_{d|n,d≠1}f(d)\cdot f^{-1}(\frac{n}{d})}{f(1)}\end{array}$

• $莫比乌斯函数$

  令 $f=\xi$,$p$为质数
  
  $f^{-1}(1)=1$
  
  $f^{-1}(p)=\frac{-f(p)\cdot f^{-1}(1)}{f(1)}=-1$
  
  $f^{-1}(p^2)=-(f(p^2)\cdot f^{-1}(1)+f(p)\cdot f^{-1}(p))=0$
  
  $f^{-1}(p^3)=-(f(p^3)\cdot f^{-1}(1)+f(p^2)\cdot f^{-1}(p)+f(p)\cdot f^{-1}(p^2))=0$
  
  $\begin{array}{r}{f}^{-1}(p^k)=-\sum _{i=1}^{k}f(p^i)\cdot f^{-1}(p^{k-i})\end{array}$
  
  其中只有 $f(p^k)\cdot f^{-1}(1)=1,f(p^{k-1})\cdot f^{-1}(p)=-1$，其余项都是0
  
  $f^{-1}(p^k)=\{\begin{array}{cc}-1& k=1\\ 0& k>1\end{array}$
  
  $n=p_1^{a_1}\cdot p_2^{a_2}\cdot \cdots p_k^{a_k}$
  
  所以莫比乌斯函数可以推出来了
  
  $\mu (n)=f^{-1}(n)=\{\begin{array}{ccc}1& n=1& \\ (-1{)}^k& a1=a2=\cdots =ak=1& \Leftarrow 积性函数推出\\ 0& otherwise& \end{array}$

---

$莫比乌斯函数的性质$

• $\mu \ast \xi =I$
  因为$\mu ={\xi }^{-1}$
  所以有$\begin{array}{r}\sum _{d|n}^n\mu (d)=\sum _{d|n}^n\mu \cdot 1=\mu \ast \xi (n)=I(n)=[n=1]\end{array}$
  记住这句，下面的证明要用到
• $I\ast f=f$
• $f\ast \mu \ast \xi =f$
• $f\ast \xi =g,f=g\ast \mu$

---

$线性筛$

利用积性函数的性质以及$\mu$的表达式

int cnt;
int prime[maxn],mu[maxn];
bool vis[maxn];
void mu_sieve (int n)
{
	mu[1]=1;
	for (int i=2;i<=n;++i){
		if (!vis[i]){	prime[++cnt]=i,mu[i]=-1;}
		for (int j=1;j<=cnt&&i*prime[j]<=n;++j){
			vis[i*prime[j]]=true;
			if (i%prime[j]==0){
				mu[i*prime[j]]=0;
				break;
			}
			mu[i*prime[j]]=-mu[i];
		}
	}
}

---

$莫比乌斯反演$

• 若有 $\begin{array}{r}g(n)=\sum _{d|n}f(d)\end{array}$
  则有 $\begin{array}{r}f(n)=\sum _{d|n}\mu (d)g(\frac{n}{d})=\sum _{d|n}\mu (\frac{n}{d})g(d)\end{array}$

  • $证明$

    法一 $\begin{array}{r}\sum _{d|n}\mu (d)g(\frac{n}{d})=\sum _{d|n}\mu (d)\sum _{x|\frac{n}{d}}f(x)=\sum _{x|n}f(x)\sum _{d|\frac{n}{x}}\mu (d)=\sum _{x|n}f(x)[\frac{n}{x}=1]=f(n)\end{array}$
    法二 $g=f\ast \xi ,f=g\ast \mu$

• 若有$\begin{array}{r}g(n)=\sum _{n|d}f(d)\end{array}$
  则有$\begin{array}{r}f(n)=\sum _{n|d}\mu (\frac{d}{n})g(d)\end{array}$

  • $证明$

    $\begin{array}{r}\sum _{n|d}\mu (\frac{d}{n})g(d)=\sum _{k=1}^{+\infty }\mu (k)g(kn)=\sum _{k=1}^{+\infty }\mu (k)\sum _{kn|d}f(d)=\sum _{n|d}f(d)\sum _{k|\frac{d}{n}}\mu (k)=\sum _{n|d}f(d)[\frac{d}{n}=1]=f(n)\end{array}$

  这一种一般用的比较多

---

$莫比乌斯反演的应用$

莫比乌斯函数一般用于$gcd$中
$\begin{array}{rl}\sum _{i=1}^n\sum _{j=1}^m\gcd (i,j)& =\sum _d^{min(n,m)}d\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=d]\\ & =\sum _d^{min(n,m)}d\sum _{i=1}^{\frac{n}{d}}\sum _{j=1}^{\frac{m}{d}}[gcd(i,j)=1]\\ & =\sum _d^{min(n,m)}d\sum _{i=1}^{\frac{n}{d}}\sum _{j=1}^{\frac{m}{d}}\sum _{k\mid gcd(i,j)}\mu (k)\\ & =\sum _d^{min(n,m)}d\sum _k^{min(\frac{n}{d},\frac{m}{d})}\mu (k)\sum _{k\mid i}^{\frac{n}{d}}\sum _{k\mid j}^{\frac{m}{d}}\\ & =\sum _d^{min(n,m)}d\sum _k^{min(\frac{n}{d},\frac{m}{d})}\mu (k)\lfloor \frac{n}{kd}\rfloor \lfloor \frac{m}{kd}\rfloor \end{array}$
令$T=kd$
$\begin{array}{r}=\sum _T^{min(n,m)}\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor \sum _{k\mid T}\frac{T}{k}\mu (k)\end{array}$
以上是莫比乌斯反演推出来的式子

$\begin{array}{r}=\sum _T^{min(n,m)}\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor \phi (T)\end{array}$
这一步是欧拉反演，证明请移步欧拉函数|(扩展)欧拉定理|欧拉反演
当然，这里这个例子不如直接用欧拉反演
$\begin{array}{rl}\sum _{i=1}^n\sum _{j=1}^m\gcd (i,j)& =\sum _{i=1}^n\sum _{j=1}^m\sum _{d|gcd(i,j)}\phi (d)=\sum _{d=1}^{min(n,m)}\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor \phi (d)\end{array}$

---

如有哪里讲得不是很明白或是有错误，欢迎指正
如您喜欢的话不妨点个赞收藏一下吧