莫比乌斯反演从零开始（数学不好的来看这个）

困扰了我半年的莫比乌斯，今天终于有了进展。本人是数论小白，数学基础不好，对数学不敏感，对数学公式推导转化头疼，所以才学的慢，下面我以一个数学不好的人的视角学习。

为了解决哪些问题？

$\sum \sum f(x)$这样的式子，比如${\sum }_{i=1}^n{\sum }_{j=1}^{n}gcd(i,j)$之类的

莫比乌斯函数

$$u(n)=\{\begin{array}{ll}1& n=1\\ (-1{)}^k& n无平方数因子，k为不同因子个数\\ 0& n有大于1的平方数因子\end{array}$$

重要的性质

$$\sum _{d|n}u(d)=[n=1]$$

莫比乌斯函数的线性筛
类似线性筛素数，线性筛出莫比乌斯函数的值及其前缀和。

const int N=1e6+10,M=1e6;
int cnt,vis[N],prime[N],Smu[N],mu[N];
void init_mu() { //线性筛莫比乌斯函数
	Smu[1]=mu[1]=1;
	cnt=0;
	for(int i=2; i<=M; i++) {
		if(vis[i]==0) {
			prime[++cnt]=i;
			mu[i]=-1;
		}
		for(int j=1; j<=cnt; j++) {
			if(i*prime[j]>M)break;
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)break;
			else mu[i*prime[j]]=-mu[i];
		}
	}
	for(int i=2;i<=M;i++)Smu[i]=Smu[i-1]+mu[i]; 
}

如何使用?

前置知识：整除分块（一学就懂）

套路1
$$[gcd(i,j)=1]=\sum _{d|gcd(i,j)}u(d)$$
$$[gcd(a_1,a_2,a_3\dots \dots a_n)=1]=\sum _{d|gcd(a_1,a_2,a_3\dots \dots a_n)}u(d)$$
套路2：枚举d，这样就不用判断 d|gcd(i,j) ，就可以将其消去
$$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|gcd(i,j)}1=\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=d]=\sum _{i=1}^{n/d}\sum _{j=1}^{m/d}=[\frac{n}{d}][\frac{m}{d}]$$
即求满足：gcd(i,j)等于d或d的倍数的对数。
证明：d|i 的 i 有$[\frac{n}{d}]$个，d|j 的 j 有$[\frac{m}{d}]$个。这些 i,j 两两组合的gcd均为d或d的倍数。

例题1.1

$$\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=1]（n<m）$$

1.由套路1变形得到
$$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|gcd(i,j)}u(d)$$
2.由套路2变形得到
$$\sum _{d=1}^{n}u(d)\ast [\frac{n}{d}][\frac{m}{d}]$$
3.u(d)的前缀和筛出来，最后整除分块即可

由于线性筛是O(n)的，所以单次查询情况下，时间复杂度并没有降下来，但在多次查询的情况下，是快了很多的。

代码如下

long long solve(long long n,long long m) {
	if(n==0||m==0)return 0;
	if(n>m)swap(n,m);
	long long ans=0;
	for(int l=1,r; l<=n; l=r+1) { //整除分块
		r=min(n/(n/l),m/(m/l));
		ans+=(Smu[r]-Smu[l-1])*(n/l)*(m/l);
	}
	return ans;
}

例题1.2

$$\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=k]（n<m）$$
变形得到
$$\sum _{i=1}^{n/k}\sum _{j=1}^{m/k}[gcd(i,j)=1]$$
就是例题1.1了。

long long solve(long long n,long long m,long long k) {
	n/=k,m/=k;
	if(n==0||m==0)return 0;
	if(n>m)swap(n,m);
	long long ans=0;
	for(int l=1,r; l<=n; l=r+1) { //整除分块
		r=min(n/(n/l),m/(m/l));
		ans+=(Smu[r]-Smu[l-1])*(n/l)*(m/l);
	}
	return ans;
}

例题1.3

$$\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=prime]（n<m）$$
1.枚举 p，变形得到
$$\sum _{p\in prime}\sum _{i=1}^{n/p}\sum _{j=1}^{m/p}[gcd(i,j)=1]$$
2.后面这个式子就是例题1嘛，变=>
$$\sum _{p\in prime}\sum _{d=1}^{n/p}u(d)\ast [\frac{n}{pd}]\ast [\frac{m}{pd}]$$
3.令k=pd，变=>
$$\sum _{p\in prime}\sum _{d=1}^{n/p}u(\frac{k}{p})\ast [\frac{n}{k}]\ast [\frac{m}{k}]$$
4.枚举k，变=>
$$\sum _{k=1}^n\sum _{p\in prime,p|k}u(\frac{k}{p})\ast [\frac{n}{k}]\ast [\frac{m}{k}]$$

5.令$f(x)={\sum }_{p\in prime,p|x}u(\frac{x}{p})$，变=>
$$\sum _{k=1}^{n}f(k)\ast [\frac{n}{k}]\ast [\frac{m}{k}]$$
6.由整除分块可知，若可以求出f(x)函数的前缀和，就可以轻而易举的求出这个式子了。怎么求？？？（等等，之后再写）

例题1.4

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

1.枚举 x，gcd(i,j)=x
$$\sum _{x=1}^{n}x\times \sum _{i=1}^n\sum _{j=1}^n[gcd(i,j)=x]$$
$$=\sum _{x=1}^{n}x\times \sum _{i=1}^{n/x}\sum _{j=1}^{n/x}[gcd(i,j)=1]$$

2.根据套路1变形
$$=\sum _{x=1}^{n}x\times \sum _{i=1}^{n/x}\sum _{j=1}^{n/x}\sum _{d|gcd(i,j)}u(d)$$

3.再根据套路2，枚举d
$$=\sum _{x=1}^{n}x\times \sum _{d=1}^{n/x}[\frac{n}{dx}][\frac{n}{dx}]u(d)$$

令$f(i)={\sum }_{d=1}^i[\frac{i}{d}][\frac{i}{d}]u(d)$，这个式子整除分块 $O(n^{1/2})$ 求
$$原式=\sum _{x=1}^{n}x\times f([\frac{n}{x}])$$
这个式子也整除分块 $O(n^{1/2})$ 求

总时间复杂度： $O(n)$

long long f(long long n) {
	long long ans=0;
	for(int l=1,r;l<=n;l=r+1){
		r=n/(n/l);
		ans+=1ll*(n/l)*(n/l)*(Smu[r]-Smu[l-1]);
	}
	return ans;
}
int main() {
	init_mu();
	long long n,ans=0;
    cin>>n;
    for(int l=1,r;l<=n;l=r+1){
    	r=n/(n/l);
    	ans+=1ll*(r-l+1)*(l+r)/2*f(n/l);
	}
	cout<<ans;
}

例题2.1

$$\sum _{i=1}^n\sum _{j=1}^{m}f(ij)，n<=m$$
【注意】其中 $f(x)$ 表示 x 的因子个数。

1.将 d(i,j) 展开

$$f(ij)=\sum _{x|i}\sum _{y|j}[gcd(x,y)=1]$$
2.原式转化，变=>
$$\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}[gcd(x,y)=1]$$

3.更换枚举项，枚举x,y，变=>
$$\sum _{x=1}^n\sum _{y=1}^m[\frac{n}{x}][\frac{m}{y}][gcd(x,y)=1]$$
$$=\sum _{i=1}^n\sum _{j=1}^m[\frac{n}{i}][\frac{m}{j}][gcd(i,j)=1]$$

4.套路1变形，变=>
$$\sum _{i=1}^n\sum _{j=1}^m[\frac{n}{i}][\frac{m}{j}]\sum _{d|gcd(i,j)}\mu (d)$$

5.与套路2的思路一致，枚举d，这样就不用考虑 d|gcd(i,j)，变=>
$$\sum _{d=1}^n\sum _{x=1}^{n/d}\sum _{y=1}^{m/d}[\frac{n}{dx}][\frac{m}{dy}]\mu (d)$$
$$=\sum _{d=1}^n\mu (d)\sum _{x=1}^{n/d}[\frac{\frac{n}{d}}{x}]\sum _{y=1}^{m/d}[\frac{\frac{m}{d}}{y}]$$

6.整除分块
优化：后面两坨的形式一样，可以使用整除分块预处理出函数 $g(x)={\sum }_{i=1}^x[\frac{x}{i}]$的所有值，预处理时间复杂度$O(n^{3/2})$。然后每次O(n)求，时间复杂度$O(Tn+n^{3/2})$。
继续优化：预处理完第一步之后，求得式子就是${\sum }_{d=1}^{n}u(d)\times g(n/d)\times g(m/d)$ ，这步还可以再整除分块。时间复杂度$O(T{n}^{1/2}+n^{3/2})$。

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+10,M=5e4;
long long mu[N],prime[N],cnt,vis[N],Smu[N],g[N];
void init_Smu() {
	Smu[1]=mu[1]=1;
	for(int i=2; i<=M; i++) {
		if(vis[i]==0) {
			prime[++cnt]=i;
			mu[i]=-1;
		}
		for(int j=1; j<=cnt; j++) {
			if(i*prime[j]>M)break;
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)break;
			else mu[i*prime[j]]=-mu[i];
		}
	}
	for(int i=2; i<=M; i++)Smu[i]=Smu[i-1]+mu[i];
	for(int i=1; i<=M; i++) {
		for(int l=1,r=2; l<=i; l=r+1) {
			r=i/(i/l);
			g[i]+=(r-l+1)*(i/l);
		}
	}
}
int main() {
	long long T,n,m;
	cin>>T;
	init_Smu();
	while(T--) {
		cin>>n>>m;
		if(n>m)swap(n,m);
		long long ans=0;
		for(int l=1,r=2;l<=n;l=r+1){
			r=min(n/(n/l),m/(m/l));
			ans+=(Smu[r]-Smu[l-1])*g[n/l]*g[m/l];
		}
		cout<<ans<<endl;
	}
	return 0;
}

例题2.2

例题链接没A掉

题目描述： T组询问，一张 n×m 的数表，其中 aij 为能同时整除 i 和 j 的所有自然数之和。给定 b，计算数表中不大于 b 的数之和。T,n<=1e5。
问题分析：
$$\sum _{i=1}^n\sum _{j=1}^{m}f(gcd(i,j))[f(gcd(i,j))<=b]，n<=m$$
【注意】其中 $f(x)$ 表示 x 的因子的和。

1.先不考虑后面那个限制，枚举 $gcd(i,j)=x$ ，并用套路1变形
$$\sum _{x=1}^n\sum _{i=1}^n\sum _{j=1}^{m}f(x)[gcd(i,j)=x]$$
=$$\sum _{x=1}^n\sum _{i=1}^{n/x}\sum _{j=1}^{m/x}f(x)[gcd(i,j)=1]$$
$$=\sum _{x=1}^{n}f(x)\sum _{i=1}^{n/x}\sum _{j=1}^{m/x}\sum _{d|gcd(i,j)}u(d)$$
2.根据套路2，枚举 d 变形
$$\sum _{x=1}^{n}f(x)\sum _{d=1}^{n/x}[\frac{n}{dx}][\frac{m}{dx}]u(d)$$
3.观察后面那个式子，令 $g(i)={\sum }_{d=1}^i[\frac{i}{d}][\frac{m/x}{d}]u(d)$ 是可以整除分块O(sqrt(n))求出来的
$$\sum _{x=1}^{n}f(x)g([\frac{n}{x}])$$
4.发现又可以整除分块。
方法1：O(n) 求出 f(x) 的时候，把 f(x)>b 的项当成0 ，并计算前缀和。TLE
方法2：对查询的 b 升序排序，用树状数组维护前缀和，最初 f(x) 为 0，随便 b 的增大，单点修改

【错误】

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10,M=1e5+5;
const long long mod=(1ll<<31);

struct quest {
	long long n,m,b,id;
} q[N];
struct node {
	long long v,id;
} f[N];
int cmp1(node x,node y) {
	return x.v<y.v;
}
int cmp(quest x,quest y) {
	return x.b<y.b;
}
long long ans[N],c[N];

long long cnt,vis[N],prime[N],Smu[N],mu[N];
void init_Smu() { //线性筛莫比乌斯函数
	Smu[1]=mu[1]=1;
	for(int i=2; i<=M; i++) {
		if(vis[i]==0) {
			prime[++cnt]=i;
			mu[i]=-1;
		}
		for(int j=1; j<=cnt; j++) {
			if(i*prime[j]>M)break;
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)break;
			else mu[i*prime[j]]=-mu[i];
		}
	}
	for(int i=2; i<=M; i++)Smu[i]=Smu[i-1]+mu[i];
}
void init_S() {
	for(int i=1; i<=M; i++) {
		for(int j=i; j<=M; j+=i) {
			f[j].v+=i;
		}
		f[i].id=i;
	}
	sort(f+1,f+M+1,cmp1);
}

int lowbit(int x) {
	return x&(-x);
}
void update(long long k,long long x) {
	for(int i=k; i<=M; i+=lowbit(i))c[i]+=x;
}
long long query(long long l,long long r) {
	long long ans=0;
	for(int i=r; i>0; i-=lowbit(i))ans+=c[i];
	for(int i=l-1; i>0; i-=lowbit(i))ans-=c[i];
	return ans;
}
long long g(long long n,long long m) {
	long long ans=0;
	for(int l=1,r; l<=n; l=r+1) {
		r=min(n/(n/l),m/(m/l));
		ans+=(Smu[r]-Smu[l-1])*(n/l)*(m/l);
		ans=ans%mod;
	}
	return ans;
}

int main() {
	freopen("1538.txt","r",stdin);
	init_S();
	init_Smu();
	long long T;
	cin>>T;
	for(int i=1; i<=T; i++) {
		scanf("%lld%lld%lld",&q[i].n,&q[i].m,&q[i].b);
		if(q[i].n>q[i].m)swap(q[i].n,q[i].m);
		q[i].id=i;
	}
	sort(q+1,q+T,cmp);

    int now=0;
	for(int i=1; i<=T; i++) {
		while(f[now+1].v<=q[i].b&&now<M) {
			now++;
			update(f[now].id,f[now].v);
		}
		long long ret=0;
		for(int l=1,r; l<=q[i].n; l=r+1) {
			r=q[i].n/(q[i].n/l);
			ret+=query(l,r)*g(q[i].n/l,q[i].m/l);
			ret=ret%mod;
		}
		ans[q[i].id]=ret;
	}
	for(int i=1; i<=T; i++)printf("%lld\n",ans[i]);
	return 0;
}

例题3

$$\sum _{i=a}^b\sum _{j=c}^d[gcd(i,j)=k]$$

1.容斥一下就是例题1.2了

$$ans=\sum _{i=1}^b\sum _{j=1}^d[gcd(i,j)=k]-\sum _{i=1}^{a-1}\sum _{j=1}^d[gcd(i,j)=k]-\sum _{i=1}^b\sum _{j=1}^{c-1}[gcd(i,j)==k]+\sum _{i=1}^{a-1}\sum _{j=1}^{c-1}[gcd(i,j)=k]$$