【数论】再谈线性筛(莫比乌斯反演)

须知

此文是作者写线性筛的第二篇文章，再读之前最好对于线性筛有一定的基础和认识，这里是作者写线性筛的第一篇文章。
前置知识

• 线性筛质数
• $\gcd$
• 一些推式子能力
• 一点数学基础

这里有一些函数定义：

1. $\phi (n)$：
2. $d(n)$
3. $(n)$

基础题(复习)

求 $F(n)={\sum }_{i=1}^n\frac{n}{\gcd (i,n)}$
输入格式：

• 第一行，一个整数 $T$。表示有 $T$ 组测试数据。
• 接下来有 $T$ 行，每行一个整数 $n$。
  全部数据：
• $1\le T\le 10^6,1\le n\le 10^7$
  输出格式：
• 共 $T$ 行，每行一个整数。

输入/输出例子
输入：

10
39
12
37
15
91
63
32
87
28
7

输出：

1099
77
1333
147
6751
2623
683
5691
473
43

样例解释

无

分析
考虑枚举 $\gcd (i,n)$，则：
$$F(n)=\sum _{d|n}\frac{n}{d}\sum _{i=1}^n[\gcd (i,n)=d]=\sum _{d|n}\frac{n}{d}\sum _{i=1}^{\frac{n}{d}}[\gcd (i,\frac{n}{d})=1]=\sum _{d|n}\frac{n}{d}\phi (\frac{n}{d})=\sum _{d|n}d\phi (d)$$
$F(n)$为积性函数
因为 $\phi (n)$ 为积性函数，所以显然 $F(n)$也是。
线性筛

1. $f(1)=1$
2. 若 $n$ 是质数，则 $f(n)=1+{\sum }_{i=1}^{n-1}n=n(n-1)+1$
3. 若 $q$ 是质数且$q\nmid n$，则 $f(qn)=f(q)f(n)$
4. 若 $q$ 是质数且$q\mid n$，则： $\{\begin{array}{c}qn=p^k\to f(qn)=f(n)+qn\phi (qn)\\ qn\ne p^k\to f(qn)=f(x)f(\frac{qn}{x}),\text{其中}x\text{为最小质因子的最高次幂的值}\end{array}$

参考代码：

#include<bits/stdc++.h>
#define int long long
#define IOS ios::sync_with_stdio(false),cin.tie(NULL),cout.tie(NULL)
using namespace std;
const int N=1e7+1,M=2e6+1,mod=1e9+7;
int t,n,m,cnt; 
struct fy{
	int prv[N],phi[N],f[N],mx[N];
	bool pr[N];
	int getphi(int x){
		int res=x;
		for(int i=2;i*i<=x;i++){
			if(x%i==0)
				res=res*(i-1)/i;
			while(x%i==0)
				x/=i;
		}
		if(x>1)
			res=res*(x-1)/x;
		return res;
	}
	void ola(int x){
		pr[1]=phi[1]=f[1]=1;
		for(int i=1;i<=x;i++){
			if(!pr[i])
				prv[++cnt]=i,phi[i]=i-1,f[i]=i*(i-1)+1,mx[i]=i;
			for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
				int u=i*prv[j];
				pr[u]=1;
				if(i%prv[j]==0){
					phi[u]=phi[i]*prv[j];
					mx[u]=mx[i]*prv[j];
					f[u]=(mx[u]==u?f[i]+phi[u]*u:f[mx[u]]*f[u/mx[u]]);
					break;
				}
				else{
					phi[u]=phi[i]*(prv[j]-1);
					mx[u]=prv[j];
					f[u]=f[i]*f[prv[j]];
				}
			}
		}
	}
}A;
signed main(){
	IOS;
	A.ola(N-1);
	cin>>t;
	while(t--)
		cin>>n,cout<<A.f[n]<<"\n";
	return 0;
}

莫比乌斯反演

莫反我自认为还是一个神奇的东西，本文对其的证明可能不好，如有错误请指出。

莫比乌斯函数

定义

考虑两个函数 $f(x),g(x)$ 满足 $f(n)={\sum }_{d|n}g(d)$
易证，一定存在某个函数 $h(x)$，使得 $g(n)={\sum }_{d|n}f(d)h(d)$
然后我们发现，对于一个定值 $d$，如果 $n$ 不同，$h(d)$ 的值也不同。（如 $g_1=f_1,g_2=f_2-f_1$）
但是呢，$h(d)$ 的值与 $\frac{n}{d}$ 有关，所以我们可以把 $h(d)$ 看成 $\mu (\frac{n}{d})$，即 $h(d)=\mu (\frac{n}{d})$，其中 $\mu$ 为一个函数。
所以可得 $g(n)={\sum }_{d|n}f(d)\mu (\frac{n}{d})$

---

这个 $\mu$ 就是莫比乌斯函数了。
它是这么定义的：
$\mu (x)=\{\begin{array}{cc}1& x=1\\ 0& \exists p\in Prime,p^2|x\\ (-1{)}^k& x={\prod }_{i=1}^k{p}_i(\text{质因数分解})\end{array}$

性质

1. $\mu$ 是积性函数。
2. $\forall n\in N_{+}$，${\sum }_{d|n}\mu (d)=[n=1]$

证明：
1 ：

• 若 $x,y$ 中有一个 $1$，易证。
• 若 $\mu (x),\mu (y)$ 中有一个 $0$，易证。
• 否则 $\mu (xy)=(-1{)}^{cn{t}_x}(-1{)}^{cn{t}_y}=\mu (x)\mu (y)$

2 ：

• 若 $n=1$，易证。
• 否则设 $n={\prod }_{i=1}^q{p}_i^{x_i}$，根据 $\mu$ 的定义，我们只考虑含有每个质因数一次的 $k(k|n)$，算 $k$ 的贡献。即： $$\sum _{d|n}\mu (d)=\mu (1)+\sum _{i=1}^q{C}_q^i(-1{)}^i{1}^{q-i}=C_q^0(-1{)}^0{1}^q+\sum _{i=1}^q{C}_q^i(-1{)}^i{1}^{q-i}=\sum _{i=0}^q{C}_q^i(-1{)}^i{1}^{q-i}=(-1+1{)}^q=0$$

线性筛

其实根据定义就可以筛了。
参考代码：

void ola(int x){
	pr[1]=mu[1]=1;
	for(int i=2;i<=x;i++){
		if(!pr[i])
			prv[++cnt]=i,mu[i]=-1;
		for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
			int u=i*prv[j];
			pr[u]=1;
			if(i%prv[j]==0){
				mu[u]=0;
				break;
			}
			else{
				mu[u]=-mu[i];
			}
		}
	}	
}

狄利克雷卷积

狄利克雷卷积可以理解为一种定义新运算，符号是 $\ast$
考虑两个函数 $f(x),g(x)$，则 $(f\ast g)(x)={\sum }_{d|x}f(d)g(\frac{x}{d})$
这就是狄利克雷卷积，其中狄利克雷卷积包含以下性质：

1. $a\ast b=b\ast a$
2. $a\ast b\ast c=a\ast (b\ast c)$

整除分块

例：
求 ${\sum }_{k=1}^n\lfloor \frac{n}{k}\rfloor$
直接上代码了，请自行理解一下啊。
Code:

for(int l=1,r;l<=n;l=r+1){
	r=(n/(n/l));
	ans+=(r-l+1)*(n/l);
}

$\text{时间复杂度}O(\sqrt{n})$

莫比乌斯反演

考虑两个函数 $f(x),g(x)$ 满足 $f(n)={\sum }_{d|n}g(d)$，则$g(n)={\sum }_{d|n}f(d)\mu (\frac{n}{d})$
正确性证明：
定义函数：$I(n)=1,ϵ(n)=[n=1],\ast$为狄利克雷卷积符号
$\therefore \mu \ast I=ϵ,h\ast ϵ=h$
$\because f=g\ast I$
$\therefore f\ast \mu =g\ast I\ast \mu \to f\ast \mu =g\ast (I\ast \mu )\to g\ast ϵ=f\ast \mu \to g=f\ast \mu$
得证

例题

题目
题意：给定 $n,m,k$，求 ${\sum }_{i=1}^n{\sum }_{j=1}^m[\gcd (i,j)=k]$
令 $n=\lfloor \frac{n}{k}\rfloor ,m=\lfloor \frac{m}{k}\rfloor$
则 $ans={\sum }_{i=1}^n{\sum }_{j=1}^m[\gcd (i,j)=1]$

直接莫反：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)=1]=\sum _{i=1}^n\sum _{j=1}^m\sum _{d|\gcd (i,j)}\mu (d)=\sum _{i=1}^n\sum _{j=1}^m\sum _{d=1}^{\min (n,m)}\mu (d)[d|i][d|j] \\ =\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{i=1}^n[d|i]\sum _{j=1}^m[d|j]=\sum _{d=1}^{\min (n,m)}\mu (d)\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor \end{gathered}$$

那么我们就可以 $O(n)$ 求答案了，但还不够，考虑用整除分块以及前缀和优化，最终复杂度 $O(\sqrt{n})$
Code:

#include<bits/stdc++.h>
#define int long long
#define IOS ios::sync_with_stdio(false),cin.tie(NULL),cout.tie(NULL)
using namespace std;
const int N=2e6+1,M=2e6+1;
int t,n,m,k,ans=0;
struct fy{
	int prv[N],cnt,mu[N],sum[N];
	bool pr[N];
	void ola(int x){
		pr[1]=mu[1]=1;
		for(int i=2;i<=x;i++){
			if(!pr[i])
				prv[++cnt]=i,mu[i]=-1;
			for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
				int u=i*prv[j];
				pr[u]=1;
				if(i%prv[j]==0){
					mu[u]=0;
					break;
				}
				else{
					mu[u]=-mu[i];
				}
			}
		}
		for(int i=1;i<=x;i++)
			sum[i]=sum[i-1]+mu[i];
			
	}
}A;
signed main(){
	IOS;
	A.ola(N-1);
	cin>>t;
	while(t--){
		cin>>n>>m>>k;
		n/=k,m/=k;
		ans=0;
		for(int l=1,r;l<=min(n,m);l=r+1){
			r=min(n/(n/l),m/(m/l));//
			ans+=(A.sum[r]-A.sum[l-1])*(n/l)*(m/l);
		}
		cout<<ans<<"\n";
	}
	return 0;
}

注意一些细节。