莫比乌斯反演新手入门及练习题

参考了以下大佬写出来的详细博客。

莫比乌斯反演入门讲解_tomandjake_的博客-CSDN博客_莫比乌斯反演

狄利克雷卷积与莫比乌斯函数_RBS的专栏-CSDN博客

在自己学习莫比乌斯反演时，自己总是会纠结为什么这一步能跳到下一步，总是希望有一个证明过程，于是每个知识自己都去学习证明过程并写下来方便自己以后回来看。整篇文章可能会比较多证明过程，但是却十分适合初学者，每一步都有详细的过程，如何得到。已经尽自己的力量尽可能详细地写出来了。

前置知识

狄利克雷卷积定义：

​ 设有两个函数g，f。
$$\begin{gathered} f\circ g(n)=\sum _{d|n}^{n}f(d)g(\frac{n}{d}) \\ d|n表示为d是n的因子，即n\%d==0 \end{gathered}$$

---

积性函数：

​ 对于完全互质的整数a，b。函数有性质$f(ab)=f(a)f(b)$。积性函数的卷积仍然为积性函数。

完全积性函数：

​ 对于任意整数a，b。函数有性质$f(ab)=f(a)f(b)$。

---

不变函数

​ 又叫全一函数，用$u$表示。 其值域始终为1，即无论n为何值，$u(n)\equiv 1$ 。 {$1,1,1,1,\cdots$}

---

恒等函数

​ 用$ID$表示。表达式为$ID(n)=n$。 {$1,2,3,4,5,\cdots$}

---

单位函数

​ 用$\epsilon$表示。只有当n=1时，$\epsilon (1)=1$, 否则，$\epsilon (n)=0$。表达式也可写为$\epsilon (n)=[n==1]$。任何函数和单位函数的狄利克雷卷积仍然得到自身。{$1,0,0,0,0,\cdots$}

---

逆函数

​ 定义一个函数$f$的逆函数$f^{-1}$,满足$f\circ f^{-1}=\epsilon$。即两者的狄利克雷卷积等于单位函数。并且积性函数的逆也是积性函数。

---

莫比乌斯函数

​ 莫比乌斯函数用$\mu (n)$表示。是一个积性函数。
$$\mu (n)=\{\begin{array}{ll}1,& n=1\\ (-1{)}^k,& n=\prod _{i=1}^k{p}_i,p互不相同\\ 0,& otherwise\end{array}$$

莫比乌斯函数完整定义的通俗表达（人话：

1）莫比乌斯函数$\mu (n)$的定义域是N；

2）$\mu (1)$=1；

3）当n存在平方因子时，$\mu (n)$=0；

4）当n是素数或奇数个不同素数之积时，$\mu (n)$=-1；

5）当n是偶数个不同素数之积时，$\mu (n)$=1。

#include <bits/stdc++.h>
using namespace std;

const int N = 1e5 + 10;
int prime[N], mobius[N], sum[N];
bool vis[N];

void isprime(int end)
{
	mobius[1] = 1;
	for(int i = 2; i <= end; i ++)
	{
		if(!vis[i])
		{
			prime[++prime[0]] = i;
			mobius[i] = -1;
		}
		for(int j = 1;j <= prime[0] && i * prime[j] <= end; j ++)
		{
			vis[i * prime[j]] = 1;
			if(i % prime[j] == 0)
				break;
			mobius[i*prime[j]] = - mobius[i]; 
		}
	}
	for(int i = 1; i <= end; i ++)
	{
		sum[i] = sum[i-1] + mobius[i];
	}
}

int main()
{
    isprime(N-10);
    return 0;
}

其具有积性的简易证明:

$\mu (n)=\mu (a)\cdot \mu (b)$，a，b互质。若a或者b有平方因子，易得两边都为0。若两者都没有平方因子，a，b提供的质因数个数可以累加至n上。

---

莫比乌斯函数性质

性质一：

$$\sum _{d|n}\mu (d)=\{\begin{array}{ll}1,& n=1\\ 0,& n\ne 1\end{array}$$

证明：

①显然当n=1时成立

②当n $\ne$ 1，将n分解成 $n=p_1^{a_1}{p}_2^{a_2}\cdots p_k^{a_k}$ 。那么在n的所有因子中，$\mu$值不为0的只有所有质因子的次数为1的因子，其中质因数个数为r的因子有$C_k^r$个。显然：
$$\sum _{d|n}=C_k^0-C_k^1+C_k^2-\cdots +(-1{)}^k{C}_k^k=0$$

---

性质二：

莫比乌斯函数的逆元为全一函数$u$。即$\mu \circ u=\epsilon$。

证明：

①当n=1时，$\mu (1)=1,u(1)=1，易得\mu \circ u=1$。即 $=\epsilon (1)$ 。

②当n $\ne$ 1时，则表达式为$\sum _{d|n}\mu (d)\cdot u(\frac{n}{d})$。将全一函数带入值，表达式化简为$\sum _{d|n}\mu (d)$。由莫比乌斯函数的性质易得此时值为0。综上，得证。

---

性质三：

恒等函数与莫比乌斯函数的卷积等于欧拉函数。

​ 欧拉函数记作$\psi (n)$，表示不超过n的与n互质的数的数量。即：$ID(n)\circ \mu (n)=\sum _{d|n}ID(d)\cdot \mu (\frac{n}{d})=\psi (n)$

证明：
$$\begin{gathered} 根据欧拉函数的定义可写出 \\ \psi (n)=\sum _i[gcd(i,n)==1] \\ 又因为只有n=1时,\mu (n)=1,故可以换成 \\ \sum _i[gcd(i,n)==1]=\sum _i\sum _{d|gcd(i,n)}\mu (d) \\ 式子左边考虑取出所有i,求和所有d。右边考虑取d,再考虑i。两者等价 \\ 考虑每一个可能的d,即d取值1,2,\cdots ,有多少对(i,j)满足i有因子d,同时n也有因子d呢？ \\ \sum _i\sum _{d|gcd(i,n)}\mu (d)=\sum _i\sum _{d|i,d|n}\mu (d)=\sum _{d|n}\frac{n}{d}\cdot \mu (d) \end{gathered}$$

---

定理一：(以下假定$\theta (n)$为积性函数，$\theta (n)=\theta (p_1^{a_1})\theta (p_2^{a_2})\cdots \theta (p_k^{a_k})$)

$$\begin{array}{rl}\psi (n)& =\sum _{d|n}\theta (d)\psi (n)为积性函数\\ & =(1+\theta (p_1)+\theta (p_1^2)+\cdots )(1+\theta (p_2)+\theta (p_2^2)+\cdots )\cdots ②\end{array}$$

证明：
$$\begin{gathered} 设gcd(m,n)=1，则mn的任意因子可写为d=d_1{d}_2。其中d_1|m,d_2|n。易得,d_1,d_2互质 \\ 若\{\begin{array}{l}{d}_1|m\\ d_2|n\end{array}\Rightarrow d_1{d}_2|mn \\ \Rightarrow \psi (mn)=\sum _{d|mn}\theta (d)=\sum _{d_1|m}\sum _{d_2|n}\theta (d_1{d}_2)=\sum _{d_1|m}\sum _{d_2|n}\theta (d_1)\theta (d_2)=\sum _{d_1|m}\theta (d_1)\sum _{d_2|n}\theta (d_2) \end{gathered}$$
由于$\sum _{d_1|m}\theta (d_1)\sum _{d_2|n}\theta (d_2)=\psi (m)\psi (n)$，故$\psi (n)$为积性函数得证。

由于$\theta (n)=\theta (p_1^{a_1})\theta (p_2^{a_2})\cdots \theta (p_k^{a_k})$，带入得到②式成立。

---

定理二：若$\theta$为积性函数，$\mu \theta$也为积性函数。

证明：定理一可知，$\sum _{d|n}u(d)\theta (d)$也是积性函数。并且
$$\begin{array}{rl}\sum _{d|n}\mu (d)\theta (d)& =(1+\mu (p_1)\theta (p_1)+\mu (p_1^2)\theta (p_1^2)+\cdots )(1+\mu (p_1)\theta (p_1)+\mu (p_1^2)\theta (p_1^2)+\cdots )\cdots \\ & \stackrel{将莫比乌斯函数的值带入得到}{=====}(1-\theta (p_1))(1-\theta (p_2))\cdots (1-\theta (p_k))\end{array}$$

莫比乌斯反演

莫比乌斯反演：设函数
$$\begin{gathered} g(n)=\sum _{d|n}f(d)\Rightarrow f(n)=\sum _{d|n}\mu (d)g(\frac{n}{d})① \\ g(n)=\sum _{n|d}f(d)\Rightarrow f(n)=\sum _{n|d}\mu (\frac{d}{n})g(d)② \end{gathered}$$
①式的证明：

​ 等价于证明 $g = f \circ u $ ，两边同时乘以莫比乌斯函数$\mu$，得到$\mu \circ g=f\circ u\circ \mu =f\circ \epsilon =f$。得证。

​ 或者
$$\sum _{d|n}\mu (d)g(\frac{n}{d})\stackrel{将已知条件带入}{==}\sum _{d|n}\mu (d)\sum _{x|\frac{n}{d}}f(x)\stackrel{d和n无区别}{=}\sum _{x|n}f(x)\sum _{d|\frac{n}{x}}\mu (d)=\sum _{x|n}f(x)[\frac{n}{x}==1]=f(n)$$

②式的证明：

​ 令$k=\frac{d}{n}$,则
$$\sum _{n|d}\mu (\frac{d}{n})g(d)=\sum _k\mu (k)g(nk)=\sum _k\mu (k)\sum _{(nk)|t}f(t)=\sum _{t}f(t)\sum _{(nk)|t}\mu (k)=\sum _{t}f(t)\epsilon (\frac{n}{t})=f(n)$$

---

常见应用

问题一：在1~N的范围里，互质的数对一共有多少对？

解：这个问题问的就是
$$\begin{gathered} \sum _i\sum _j[gcd(i,j)==1] \\ 仿照性质三的推导 \\ =\sum _i\sum _j\sum _{d|gcd(i,j)}\mu (d) \\ =\sum _d^N\mu (d)\cdot \lfloor \frac{N}{d}\rfloor \cdot \lfloor \frac{N}{d}\rfloor \end{gathered}$$
或者（有关这类的问题都可以这么设函数。）
$$\begin{gathered} 设f(n)=\sum _i^N\sum _j^N[gcd(i,j)==n] \\ g(n)=\sum _i^N\sum _j^N[n|gcd(i,j)] \\ 那么g(n)=\sum _{n|d}^{N}f(d)且g(n)=\lfloor \frac{N}{n}\rfloor \lfloor \frac{N}{n}\rfloor \\ 根据反演,得到f(n)=\sum _{n|d,d⩽N}\mu (\frac{d}{n})f(d)=\sum _{n|d,d⩽N}\mu (\frac{d}{n})\lfloor \frac{N}{d}\rfloor \lfloor \frac{N}{d}\rfloor \\ 题目要求是f(1)=\sum _d^N\mu (d)\cdot \lfloor \frac{N}{d}\rfloor \cdot \lfloor \frac{N}{d}\rfloor \end{gathered}$$
问题二：$i$在1～N范围内，$j$在1～M范围内，问互质的数对的数量有多少对？
$$\begin{gathered} 问题问的是\sum _i^N\sum _j^M[gcd(i,j)==1]=\sum _i^N\sum _j^M\sum _{d|gcd(i,j)}\mu (d) \\ =\sum _d^{d⩽min(M,N)}\mu (d)\lfloor \frac{N}{d}\rfloor \cdot \lfloor \frac{M}{d}\rfloor \end{gathered}$$
问题三：在1～N的范围内，gcd为质数的数对一共有多少对？即问：
$$\begin{array}{rl}\sum _{p_k}\sum _i^N\sum _j^N[gcd(i,j)==p_k]& =\sum _{p_k}\sum _i^N\sum _j^N[gcd(\frac{i}{p_k},\frac{i}{p_k})==1]\\ & =\sum _{p_k}\sum _i^N\sum _j^N\sum _{d|gcd(\frac{i}{p_k},\frac{j}{p_k})}\mu (d)\\ & =\sum _{p_k}\sum _d\mu (d)\lfloor \frac{\lfloor \frac{N}{p_k}\rfloor }{d}\rfloor \lfloor \frac{\lfloor \frac{N}{p_k}\rfloor }{d}\rfloor \end{array}$$

习题

能量采集

YY的GCD

余数求和

公约数的和

奇妙数论题

数字表格

Sky Code - POJ