莫比乌斯函数与莫比乌斯反演

【目录】

• 莫比乌斯函数
• 莫比乌斯反演

---

莫比乌斯函数

定义

莫比乌斯函数$\mu (n)$，当$n=1$时，$\mu (n)=1$；当$n>1$时，设$n$的唯一分解式为$n=p_1^{c_1}\cdots p_k^{c_k}$，则$\mu (n)$定义为
$$\mu (n)=\{\begin{array}{l}(-1{)}^k,c_1=c_2=\cdots =c_k=1\\ 0,\exists c_i>1(1\le i\le k)\end{array}$$

性质

• $\sum _{d|n}\mu (d)=[n=1]$
  注：约定方括号[]中为一个命题，其结果为$1$（该命题为真），或为$0$（该命题为假）。例如$[p为质数]$=$\{\begin{array}{l}1,p是质数\\ 0,p不是质数\end{array}$
  证明：$n=1$时，显然成立；现设$n>1$，$n$的唯一分解式为$n=p_1^{c_1}\cdots p_k^{c_k}$，则
  $$\begin{array}{rl}\sum _{d|n}\mu (d)=& \mu (1)+\mu (p_1)+\cdots +\mu (p_k)+\mu (p_1{p}_2)\\ & +\cdots +\mu (p_{k-1}{p}_k)+\cdots +\mu (p_1\cdots p_k)\\ =& 1+\left(\genfrac{}{}{0ex}{}{k}{1}\right)(-1)+\left(\genfrac{}{}{0ex}{}{k}{2}\right)(-1{)}^2+\cdots +(\genfrac{}{}{0ex}{}{k}{k})(-1{)}^k\\ =& (1-1{)}^k=0\end{array}$$
• $\phi (n)=\sum _{d|n}\mu (d)\frac{n}{d}$
  证明：
  因为$\phi (n)=n\left(1-\frac{1}{p_1}\right)\cdots \left(1-\frac{1}{p_k}\right)$，其中$n=p_1^{c_1}\cdots p_k^{c_k}$是$n$的标准分解式，利用$\mu (n)$，可得$\phi (n)=\sum _{d|n}\mu (d)\frac{n}{d}$。

---

莫比乌斯反演

观察这两个等式
$$\begin{array}{rl}n& =\sum _{d|n}\phi (d)=\sum _{d|n}\phi \left(\frac{n}{d}\right)\\ \phi (n)& =\sum _{d|n}\mu (d)\frac{n}{d}=\sum _{d|n}\mu \left(\frac{n}{d}\right)d\end{array}$$
考虑将其推广至一般情况

莫比乌斯变换

对于数论函数$f(n),g(n)$，若
$$f(n)=\sum _{d|n}g(d)=\sum _{d|n}g\left(\frac{n}{d}\right)$$
则称$f(n)$为$g(n)$的莫比乌斯变换，而$g(n)$为$f(n)$的莫比乌斯逆变换。

反演公式

若有两个数论函数$f(n),g(n)$满足
$$f(n)=\sum _{d|n}g(d)(1)$$
则有
$$g(n)=\sum _{d|n}\mu (d)f\left(\frac{n}{d}\right)(2)$$
反过来，若满足$(2)$，则$(1)$也成立。
证明： 若$f(n),g(n)$满足$(1)$，则
$$\begin{array}{rl}\sum _{d|n}\mu (d)f\left(\frac{n}{d}\right)& =\sum _{d|n}\mu (d)\sum _{d^{\prime }|\frac{n}{d}}g(d^{\prime })\\ & =\sum _{d{d}^{\prime }|n}\mu (d)g(d^{\prime })\\ & =\sum _{d^{\prime }|n}\sum _{d|\frac{n}{d^{\prime }}}\mu (d)g(d^{\prime })\\ & =\sum _{d^{\prime }|n}g(d^{\prime })\sum _{d|\frac{n}{d^{\prime }}}\mu (d)\\ & =g(n)\end{array}$$
反过来，设$f(n),g(n)$满足$(2)$，同法可证
$$\begin{array}{rl}\sum _{d|n}g(d)& =\sum _{d|n}g\left(\frac{n}{d}\right)\\ & =\sum _{d|n}\sum _{d^{\prime }|\frac{n}{d}}\mu \left(\frac{n}{d{d}^{\prime }}\right)f(d^{\prime })\\ & =\sum _{d{d}^{\prime }|n}\mu \left(\frac{n}{d{d}^{\prime }}\right)f(d^{\prime })\\ & =\sum _{d^{\prime }|n}f(d^{\prime })\sum _{d|\frac{n}{d^{\prime }}}\mu \left(\frac{n}{d{d}^{\prime }}\right)\\ & =f(n)\end{array}$$

在OI中的应用

通常，在OI竞赛中，应用莫比乌斯反演的关键在于构造如下的式子
$$\sum _{i=1}^n\sum _{j=1}^{m}f(\gcd (i,j))$$
其中$f(n)$是一个积性函数。
构造数论函数$g(n)$满足$f(n)=\sum _{d|n}g(d)$，
由莫比乌斯反演公式得$g(n)=\sum _{d|n}\mu (d)f\left(\frac{n}{d}\right)$。
化简原式
$$\sum _{i=1}^n\sum _{j=1}^{m}f(\gcd (i,j))=\sum _{i=1}^n\sum _{j=1}^m\sum _{d|\gcd (i,j)}g(d)$$
因为$d|\gcd (i,j)\iff d|i,d|j$，所以$d$必须是$i,j$的约数
考虑对每个$d$，枚举$d$的倍数，接着化简
$$\begin{array}{rl}\sum _{i=1}^n\sum _{j=1}^{m}f(\gcd (i,j))& =\sum _{d=1}^{\min (n,m)}\sum _{d|i}^n\sum _{d|j}^{m}g(d)\\ & =\sum _{d=1}^{\min (n,m)}\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor g(d)\end{array}$$
这样只需要枚举$d$，就能求出$\sum _{i=1}^n\sum _{j=1}^{m}f(\gcd (i,j))$，时间复杂度为$O(n)$。
考虑$\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$可以使用数论分块，再预处理一下$g(n)$的前缀和，时间复杂度降至$O(\sqrt{n})$。
至于$g(n)$的计算，因为$g(n)=\sum _{d|n}\mu (d)f\left(\frac{n}{d}\right)$，而$f(n),\mu (n)$为积性函数，所以$g(n)$也是积性函数。参考这篇博客，$g(n)$可以在线性时间内求出。