浅谈莫比乌斯反演

Latex炸了，可以看这里

α.莫比乌斯函数的性质

性质1

莫比乌斯函数是积性函数，即：

若 $x,y$ 互质，则$\mu (xy)=\mu (x)\mu (y)$

证：

考虑讨论 $\mu (x)$ 和 $\mu (y)$ 的值

1.若 $\mu(x)=$0或$\mu(y)=0$，则 $x,y$ 两数中有平方因子，所以 $xy$ 必然存在平方因子，所以 $\mu (xy)=0$。

2.若 $\mu (x)\ne 0且\mu (y)\ne 0$，则 $x,y$ 的质因子都只出现了一次，又因为 $x,y$ 互质，所以 $xy$ 的每个质因子都只出现了一次，所以 $\mu (xy)=\mu (x)\mu (y)$

利用莫比乌斯函数是不完全积性函数的特点，我们可以用线性筛来快速求出莫比乌斯函数的值，代码如下：

const int N = 1e6+5;

int cnt,fl[N],pr[N],mu[N];

void init(int n){
	mu[1]=1;
	for(int i=2;i<=n;i++){
		if(!fl[i]) pr[++cnt]=i,mu[i]=-1;//质数只有一个质因子，所以mu[i]为-1
		for(int j=1;j<=cnt&&i*pr[j]<=n;j++){
			vis[i*pr[j]]=1;
			if(i%pr[j])mu[i*pr[j]]=-mu[i];//i*pr[j]相较于i多了一个质因子，所以值取反
			else{
				mu[i*pr[j]]=0;//pr[j]是i的因子，所以i*pr[j]就至少含有两个pr[j]，即存在平方因子 
				break;
			} 
		} 
	}
}

性质2

若 $n=1$，则 ${\sum }_{d|n}\mu (d)=1$

若 $n>1$，则 ${\sum }_{d|n}\mu (d)=0$

证：

当 $n=1$ 时，显然成立

当 $n>1$ 时，令 $n=\prod p_k^{k_k}$，$n$ 的正因子 $d=\prod p_k^{q_k}$

因为质因子幂次数大于 $1$ 的 $\mu$ 值为 $0$, 所以只需考虑幂为 $0,1$ 的质因子项,所以：

$$\sum _{d|n}\mu (d)=\sum _{k=0}^m\left(\begin{array}{c}m\\ k\end{array}\right)\times (-1{)}^k=\sum _{k=0}^m\left(\begin{array}{c}m\\ k\end{array}\right)\times (-1{)}^k\times 1^{m-i}$$

根据二项式定理，上述式子是 $(1-1{)}^m$ 的二项式展开，所以原式为 $0$

性质3

${\sum }_{i=1}^n{\sum }_{d|n}\mu (d)={\sum }_{d=1}^n\lfloor \frac{n}{d}\rfloor \mu (d)$

证：

反向思考，满足是 $d$ 的倍数的 $i$ 在 $[1,n]$ 中只出现了 $\lfloor \frac{n}{d}\rfloor$ 次。

性质4

${\sum }_{d|n}\frac{\mu (d)}{d}=\frac{\phi (n)}{n}$

β.莫比乌斯函数的经典运用

运用1

试求 ${\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=1]$的值

不妨设，$n<m$，根据性质2，将 $[gcd(i,j)=1]$ 替换为 ${\sum }_{d|gcd(i,j)}\mu (d)$，则原式化为：

$$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|gcd(i,j)}\mu (d)$$

我们又发现，$d|gcd(i,j)$的充分必要条件是 $d|i$ 且 $d|j$，所以原式又可以化作：

$$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|i}\sum _{d|j}\mu (d)$$

前后调整求和符号的位置，得：

$$\sum _{i=1}^n\sum _{d|i}\sum _{j=1}^m\sum _{d|j}\mu (d)$$

根据性质3，化简原式得：

$$\sum _{d=1}^n\lfloor \frac{n}{d}\rfloor \sum _{d=1}^m\lfloor \frac{m}{d}\rfloor \mu (d)$$

又发现前后两个 $\sum$ 中的 $d$ 是一个 $d$，所以只能取到$\min (n,m)$ 即 $n$，将两个 $\sum$ 合并得到(前面我们假定了 $n<m$)：

$$\sum _{d=1}^n\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor \mu (d)$$

运用2：

试求 ${\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=k]$的值

考虑对原式整体除以 $k$，得到 ${\sum }_{i=1}^{\lfloor \frac{n}{k}\rfloor }{\sum }_{j=1}^{\lfloor \frac{m}{k}\rfloor }[gcd(i,j)=1]$

现在我们就已经将问题转化为运用1了，最终的结果为：
$$\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\lfloor \frac{n}{kd}\rfloor \lfloor \frac{m}{kd}\rfloor \mu (d)$$

运用3

试求 ${\sum }_{i=1}^n{\sum }_{j=1}^{m}ij\times [gcd(i,j)=k]$ 的值

还是先同运用2，将原式化简，不过要注意的是，要另外考虑 $i,j$ 的贡献，由于我们对原式整体除以 $k$，$i$ 和 $j$ 都变为了原大小的 $\frac{1}{k}$ 的应对结果乘以 $k^2$，即原式化为：

$$\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }ij{k}^2\times [gcd(i,j)=1]$$

还是老套路，将 $[gcd(i,j)=1]$ 替换为 ${\sum }_{d|gcd(i,j)}\mu (d)$，拆开：

$$\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }ij{k}^2\sum _{d|i}\sum _{d|j}\mu (d)$$

将前面 $\sum$ 拆开，得：

$$\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }ij{k}^2\times \sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }\sum _{d|i}\sum _{d|j}\mu (d)$$

同运用1化简后面的四个 $\sum$，得到：

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

合并，并将前面的两个 $\sum$ 提出两个 $d$ 出来，并将 $d^2$ 置于后面，得到：

$$\sum _{i=1}^{\lfloor \frac{n}{kd}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{kd}\rfloor }ij{k}^2\times \sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\times d^2$$

最后发现，几个 $\sum$ 之间没有依赖关系，将它们分开：

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

运用4

前置知识

${\sum }_{d|n}\phi (d)=n$

证：

先令 $f(x)={\sum }_{d|x}\phi (d)=n$，对于 $\forall i,j$ 满足 $gcd(i,j)=1$，有 $f(ij)={\sum }_{d|ij}\phi (d)={\sum }_{d|i}\phi (d)\times {\sum }_{d|j}\phi (d)=f(i)f(j)$，由此可知 $f$ 是不完全积性函数。

将 $n$ 进行质因数分解，假定 $n=\prod p_k^{q_k}$，考虑单独计算每一个 $f(p_k^{q_k})$：

$f(p_k^{q_k})={\sum }_{i=0}^{q_k}\phi (p_k^i)=1+{\sum }_{i=1}^{q_k}{p}_k^{i-1}\times (p_k-1)$

根据等比数列求和公式可得上述式子等于 $p_k^{q_k}$，即 $f(p_k^{q_k})=p_k^{q_k}$，所以 ${\sum }_{d|n}\phi (d)=f(n)=f(\prod p_k^{q_k})=\prod f(p_k^{q_k})=\prod p_k^{q_k}=n$

正文

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

有了上面的铺垫，我们很容易想到将 $\text{gcd}(i,j)$ 替换成 ${\sum }_{d|\text{gcd}(i,j)}\phi (d)$：

$$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|\text{gcd}(i,j)}\phi (d)$$

同运用1，将 $d|\text{gcd}(i,j)$ 的条件拆开，并调换顺序，得：

$$\sum _{i=1}^n\sum _{d|i}\sum _{j=1}^m\sum _{d|j}\phi (d)$$

继续同运用1，得到：

$$\sum _{d=1}^n\phi (d)\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$$