莫比乌斯反演 超详细

前置芝士：数论分块。

莫比乌斯函数

定义莫比乌斯函数：$\mu (n)=\{\begin{array}{l}1n=1\\ 0n含有完全平方因子\\ (-1{)}^{k}k为n的不同质因子个数\end{array}$

根据神神奇奇的定义式，易得它是一个积性函数。那么就可以用线性筛快速求解：

	for(int i=2;i<=n;++i){
		if(!vis[i]) p[++tot]=i,mu[i]=-1;
		for(int j=1;j<=tot&&1ll*p[j]*i<=n;++j){
			int now=i*p[j];
			vis[now]=1;
			if(i%p[j]==0){mu[now]=0;break;}
			mu[now]=mu[i]*mu[p[j]];
		}
	}

它有一个非常重要的性质：

$\sum _{d\mid n}\mu (d)=\{\begin{array}{l}1n=1\\ 0n\ne 1\end{array}\Rightarrow \sum _{d\mid \gcd (i,j)}\mu (d)=\{\begin{array}{l}1\gcd (i,j)=1\\ 0\gcd (i,j)\ne 1\end{array}$

利用这个性质，我们可以通过一系列推柿子来降低一些求 $gcd$、$lcm$ 和的问题的时间复杂度。这就是莫比乌斯反演。

莫比乌斯变换

• 对于两个数论函数 $f(n),g(n)$，若 $f(n)=\sum _{d\mid n}g(d)$，则 $g(n)=\sum _{d\mid n}f(d)\mu (\frac{n}{d})$。

一些式子

其实都是一些比较套路的转化，多做几道就会了qwq

• 求 $\sum _{i=1}^n\sum _{j=1}^m(\gcd (i,j)=k)$

首先，套路性的将 $d$ 转化为 $1$：$\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }(\gcd (i,j)=1)$

然后，再套路性的引入莫反：$\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }\sum _{d\mid \gcd (i,j)}\mu (d)$

接着，再套路性的将 $d$ 提前（假设 $n<m$）：$\sum _{d=1}^n\mu (d)\sum _{i=1}^{\lfloor \frac{n}{d\times k}\rfloor }\sum _{i=1}^{\lfloor \frac{m}{d\times k}\rfloor }1$

于是我们就得到了最终的柿子：$\sum _{d=1}^n\mu (d)\lfloor \frac{n}{d\times k}\rfloor \lfloor \frac{m}{d\times k}\rfloor$，只要预处理莫比乌斯函数的前缀和，再数论分块，就可以做到 $O(\sqrt{n})$ 的复杂度了。

• 求 $\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)$

$lcm$ 不好求，我们可以转化为 $gcd$：$\sum _{i=1}^n\sum _{j=1}^m\frac{i\times j}{\gcd (i,j)}$

还是套路性的，枚举 $gcd$（依然假设 $n<m$）：$\sum _{k=1}^n\sum _{i=1}^n\sum _{j=1}^m\frac{i\times j}{k}(\gcd (i,j)=k)$

接着将 $d$ 化为 $1$：$\sum _{k=1}^n\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }\frac{i\times j}{k}\times k^2(\gcd (i,j)=1)$

引入 $\mu$：$\sum _{k=1}^n\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }i\times j\times k(\sum _{d\mid \gcd (i,j)}\mu (d))$

将 $d$ 提前：$\sum _{k=1}^n\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\sum _{i=1}^{\lfloor \frac{n}{k\times d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k\times d}\rfloor }i\times j\times k\times d^2$

这里有是一个常见套路：将 $k\times d$ 用 $T$ 代换：$\sum _{T=1}^{n}T\times \sum _{d\mid T}d\times \mu (d)\sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{T}\rfloor }i\times j$

再变一下形：$\sum _{T=1}^{n}T\times (\sum _{d\mid T}d\times \mu (d))(\sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }i)(\sum _{j=1}^{\lfloor \frac{m}{T}\rfloor }j)$

发现后两个括号就是前缀和，而第一个括号可以枚举 $d$ 预处理，根据某些神奇的积分知识可以证明这样是 $O(n\log \log n)$ 的。但这样在某些题里会被卡，又发现 $f(n)=\sum _{d\mid n}d\times \mu (d)$ 是一个积性函数，是可以线性筛的，这样就做到了 $O(n)$。

inline void shai(int n){
	val[1]=1;
	rep(i,2,n){
		if(!vis[i]) p[++tot]=i,val[i]=(mod+1-i)%mod;
		for(int j=1;j<=tot&&1ll*i*p[j]<=n;++j){
			int now=i*p[j];vis[now]=1;
			if(i%p[j]==0){val[now]=val[i];break;}
			val[now]=1ll*val[i]*val[p[j]]%mod;
		}
	}
}

还有一类 $lcm$ 相关题型给定数组 $a$，求 $\sum _{i=1}^n\sum _{j=1}^{m}lcm(a_i,a_j)$，这样的问题一般 $a$ 的值域不超过 $1e6$，可以考虑开一个桶 $c$ 记录每个数出现了几次，假设值域上界为 $n$，答案就是 $\sum _{T=1}^{n}T\times (\sum _{d\mid T}d\times \mu (d))(\sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }i\times c_{iT})(\sum _{j=1}^{\lfloor \frac{n}{T}\rfloor }j\times c_{jT})$，发现后两个括号虽然变了，但每次暴力计算也不会超时 （对于这些神奇循环，其实笔者只知道是 O(能过)qwq） ，于是可解。

• 设 $d(x)$ 为 $x$ 的约数个数，求 $\sum _{i=1}^n\sum _{j=1}^{m}d(ij)$

首先得到 $d$ 的一个性质：$d(ij)=\sum _{x\mid i}\sum _{y\mid j}(\gcd (x,y)=1)$
简单证明：假设 $i\times j$ 含有因子 $p^c$，$p^a\mid i$ 且 $p^{a+1}\nmid i$，$p^b\mid j$ 且 $p^{b+1}\nmid j$，那么在 $i\times j$ 中有质因子 $p$ 有 $a+b+1$ 种选择方案，设 $(x,y)$ 表示 $i$ 中选 $x$ 个，$j$ 中选 $y$ 个，则等式右边的选择有 $([1,a],0)$、$(0,[1,b])$、$(0,0)$ 共 $a+b+1$ 种方案，等式成立。

转化完 $d$ 函数就可以按套路推柿子了：$\sum _{i=1}^n\sum _{j=1}^m\sum _{x\mid i}\sum _{y\mid j}(\gcd (x,y)=1)$

形如 $x\mid i$ 的柿子不好求，我们改成枚举 $x$ 和 $y$：$\sum _{x=1}^n\sum _{y=1}^m\sum _{i=1}^{\lfloor \frac{n}{x}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{y}\rfloor }\gcd (x,y)=1$

依然是 $gcd$ 的转化，顺便消一些 $i$ 和 $j$：$\sum _{x=1}^n\sum _{y=1}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor \sum _{d\mid \gcd (x,y)}\mu (d)$

依然是将 $d$ 提前（还是假设 $n<m$）：$\sum _{d=1}^n\mu (d)\sum _{x=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{y=1}^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{n}{x\times d}\rfloor \lfloor \frac{m}{y\times d}\rfloor$

打打括号：$\sum _{d=1}^n\mu (d)(\sum _{x=1}^{\lfloor \frac{n}{d}\rfloor }\lfloor \frac{n}{x\times d}\rfloor )(\sum _{y=1}^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{m}{y\times d}\rfloor )$

这就又变成了线性筛+数论分块的经典柿子。

可以发现其实不同的推柿子都有一些相同的套路性步骤，只要手感出来了就万事大吉了qwq。

这些柿子在洛谷的对应题目：

• [POI2007]ZAP-Queries
• 最小公倍数之和
• 约数个数和

---

记个结论：$[k$ 的最大完全平方因子为 $1]=\sum _{i^2\mid k}\mu (i)$