【学习笔记】莫比乌斯函数

传送门 to PoPoQQQ神犇，膜拜一下。

概述

一个神奇的函数。定义也很神奇。是个积性函数。

定义

$\mu (1)=1$ 。对于 $x(x>1)$，若 $x=p_1^{t_1}{p}_2^{t_2}{p}_3^{t_3}\dots p_m^{t_m}$ （指 $x$ 的唯一素数分解），$$\mu (x)=\{\begin{array}{ll}(-1{)}^m& (\forall y\in [1,m],t_y=1)\\ 0& (\exists y\in [1,m],t_y>1)\end{array}$$

应用

思想是类似 容斥 的玩意儿。（看到负一的多少多少次方，很有可能是容斥嘛！）

莫比乌斯反演

考虑这样的定义：

$$F(x)=\sum _{y|x}f(y)$$

现在我们要想办法用 $F$ 表示 $f$ 。

1. $F(x)$ 中有 $f(x)$ 这一项。先拿一个 $F(x)$ 。
2. 多出来的是$x$的其他因子。只好减掉。

怎么枚举其他因子？只需要在唯一素数分解中任选一个质因数，将其指数减少一即可。

但是，如果有两个质因数，其指数都少了一，怎么办？加回来。

或者用韦恩图来表示这一关系。不妨设 $x$ 只有三个质因数。

$y$ 是 $x$ 的因数，充要条件是：唯一素数分解下，$y$ 的每一个质因数的指数都不大于 $x$ 的对应指数。

所以，$x$ 的全部因数（不包含 $x$ ），就是上图中的全集（至少有一个质因数的指数是严格小于）。此时便可容斥。

最后写结论。$$f(x)=\sum _{y|x}\mu (y)F\left(\frac{x}{y}\right)$$

还有另一种形式，也是容斥：$$F(d)=\sum _{d|n}f(n)\iff f(d)=\sum _{d|n}\mu \left(\frac{n}{d}\right)F(n)$$

其他的函数

比如这样的函数：
$$S(n,p)=\sum _{x=1}^n[\gcd (x,p)=1]$$

设 $p=p_1^{t_1}{p}_2^{t_2}{p}_3^{t_3}\cdots p_m^m(p_x\in \mathbb{P})$ ，$\mathbb{P}$ 是质数集。

考虑 $p_x$ 的倍数，必然不与 $p$ 互质。这样的数有 $\frac{p}{p_x}$ 个。

若减去 $\sum \frac{p}{p_x}$ ，会重复。有些数又是 $p_x$ 的倍数，也是 $p_y(x\ne y)$ 的倍数。需减去。

容斥思想很明显了嘛！

注意 $x\ne y$ 。因为 $p_x^2$ 的倍数没有被多减去。

$n(n\le 1000000)$ 最多有 $7$ 个不同的质因数，复杂度就会是 $\mathcal{O}(2^7)$ 。——事实上不会总是七个质因数的。

最后写结论。$$S(n,p)=\sum _{d|p}\mu (d)\lfloor \frac{n}{d}\rfloor$$

奇怪的玩意儿

传说中的 $$\sum _{d|n}\frac{\mu (d)}{d}=\frac{\phi (n)}{n}$$

证明很简单，注意到 $n={\sum }_{d|n}\phi (d)$ ，代入反演公式即可。

还有一个 $$\sum _{x|n}\mu (x)=ϵ(n)$$

证明仍然很简单，利用 $\mu (n)={\sum }_{d|n}\mu (d)ϵ(\frac{n}{d})$ 进行反演的逆推。

计算

线性筛

可以在线性筛的同时一起计算，因为这是一个积性函数。

mu[1] = 1; // here it is!
for(int i=2; i<=n; ++i){
	if(isPrime[i]){
		primes.push_back(i);
		mu[i] = -1; // here it is!
	}
	for(int j=0,len=primes.size(); j<len; ++j){
		if(1ll*i*primes[j] > n) break;
		isPrime[i*primes[j]] = false;
		mu[i*primes[j]] = -mu[i]; // here it is!
		if(i%primes[j] == 0){
			mu[i*primes[j]] = 0; // here it is!
			break;
		}
	}
}

杜教筛

利用这条性质：$$\sum _{d|n}\mu (d)=[n=1]$$

用狄利克雷卷积写出来就是$$\mu \ast I=ϵ$$

所以

$$\begin{gathered} \sum _{i=1}^nϵ(n)=\sum _{i=1}^{n}I(i)S\left(\lfloor \frac{n}{i}\rfloor \right) \\ \Rightarrow S(n)=1-\sum _{i=2}^{n}S\left(\lfloor \frac{n}{i}\rfloor \right) \end{gathered}$$

map<int,int> mp;
int sumMu(int n){
	if(n < sqrtN) return mu[n];
	if(mp.count(n)) return mp[n];
	int res = 1;
	for(long long l=2,r; l<=n; l=r){
		r = n/(n/l)+1;
		res = (res-(r-l)*sumMu(n/l)%Mod+Mod)%Mod;
	}
	return mp[n] = res;
}