莫比乌斯反演基础

一、先补充几个常用和式规则

$$\begin{gathered} \sum _{m|n}{a}_m=\sum _{m|n}{a}_{n/m} \\ \sum _{m|n}{a}_m=\sum _k\sum _{m>0}{a}_m[n=mk] \\ \sum _{m|n}\sum _{k|m}{a}_{k,m}=\sum _{k|n}\sum _{l|(n/k)}{a}_{k,kl} \end{gathered}$$

二、积性函数的概念

​ 如果 $a,b$ 互素，$f(ab)$ = $f(a)$ $\ast$ $f(b)$ ，就称 $f$ 为积性函数

​ 如果 $p$ 为素数，并且 $f(p^n)=f^n(p)$ 则f为完全积性函数

三、莫比乌斯函数的概念

$$\sum _{d|m}\mu (d)=[m=1]$$

​ 反演原理：
$$\begin{gathered} g(m)=\sum _{d|m}f(d)\Leftarrow \Rightarrow f(m)=\sum _{d|m}\mu (d)g(\frac{m}{d}) \\ g(n)=\sum _{n|d}f(d)\Leftarrow \Rightarrow f(n)=\sum _{n|d}\mu (\frac{d}{n})g(d) \end{gathered}$$
​ 证明如下：
$$\begin{gathered} 如果g(m)=\sum _{d|m}f(d)，那么 \\ \sum _{d|m}\mu (d)g(\frac{m}{d})=\sum _{d|m}\mu (\frac{m}{d})g(d) \\ =\sum _{d|m}\mu (\frac{m}{d})\sum _{k|d}f(k) \\ =\sum _{k|m}\sum _{d|(m/k)}\mu (\frac{m}{kd})f(k) \\ =\sum _{k|m}\sum _{d|(m/k)}\mu (d)f(k) \\ =\sum _{k|m}[m/k=1]f(k)=f(m) \end{gathered}$$

第二个式子的证明和第一个基本一致

同时莫比乌斯函数还有另外一个性质
$$g(x)=\sum _{d\ge 1}f(x\ge d)\Leftarrow \Rightarrow f(x)=\sum _{d\ge 1}\mu (d)g(x/d)$$
这个反演法则对所有满足 ${\sum }_{k,d\ge 1}|f(x/kd)|<\infty$

证明如下：先假设 $g(x)={\sum }_{d\ge 1}f(x/d)$ ，那么
$$\begin{gathered} \sum _{d\ge 1}\mu (d)g(x/d)=\sum _{d\ge 1}\mu (d)\sum _{k\ge 1}f(x/kd) \\ =\sum _{m\ge 1}f(x/m)\sum _{d,k\ge 1}\mu (d)[m=kd] \\ =\sum _{m\ge 1}f(x/m)\sum _{d|m}\mu (d) \\ =\sum _{m\ge 1}f(x/m)[m=1] \\ =f(x) \end{gathered}$$

四、整除分块(数论分块)

$$对于\sum _{i=1}^n\lfloor n/i\rfloor$$

循环遍历求当然可以，但是如果 $n$ 很大，并且配合其它的式子求解的话，时间复杂度就会非常高

首先我们想一下，对于每一个 $i$ 都是 $n/i$的下取整，这个有什么性质呢

比如 $n/i$的下取整得到了一个 $k$，那么这个 $i$ 一定小于等于 $n/k$， 这也就是说从 $i$ 到 $n/k$得到的答案一定都是 $k$

代码如下：

for(int i = 1; i <= n; ) {
    int r = n / (n / i);
    sum2 += (n / i) * (r - i + 1);
    i = r + 1;
}

五、例题

1.$P3455[POI2007]ZAP-Queries$

题目描述

密码学家正在尝试破解一种叫 $BSA$ 的密码。

他发现，在破解一条消息的同时，他还需要回答这样一种问题：

给出 $a,b,d$，求满足 $1\le x\le a$，$1\le y\le b$且 $gcd(x,y)=d$ 的二元组的数量。

因为要解决的问题实在太多了，他便过来寻求你的帮助。

输入格式

输入第一行一个整数 nn，代表要回答的问题个数。

接下来 $n$ 行，每行三个整数 $a,b,d$。

输出格式

对于每组询问，输出一个整数代表答案。

输入 #1

2
4 5 2
6 4 3

输出 #1

3
2

说明

对于全部的测试点，保证 $1\le n\le 5\ast 10^4$, $1\le d\le a$, $b\le 5\ast 10^4$。

解题思路：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==d] \\ =\sum _{i=1}^n\sum _{j=1}^m[gcd(i/d,j/d)==1] \\ =\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }[gcd(i,j)==1] \\ =\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }\sum _{k|gcd(i,j)}\mu (k) \\ =\sum _{k=1}^{min(\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )}\mu (d)\ast \lfloor \frac{n}{kd}\rfloor \ast \lfloor \frac{m}{kd}\rfloor \end{gathered}$$
然后利用数论分块的知识就可以AC了

链接：https://www.luogu.com.cn/problem/P3455

AC代码

#include<cstdio>
#include<algorithm>

using namespace std;

const int N = 1e5 + 10;

bool vis[N];
int prime[N], mu[N], cnt, sum[N];

void getMu() {
	mu[1] = 1;
	for(int i = 2; i < N - 1; i++) {
		if(!vis[i]) {
			prime[cnt++] = i;
			mu[i] = -1;
		}
		for(int j = 0; j < cnt && i * prime[j] < N - 1; j++) {
			vis[i * prime[j]] = 1;
			mu[i * prime[j]] = -mu[i];
			if(i % prime[j] == 0) {
				mu[i * prime[j]] = 0;
				break;
			}
		}
	}
	for(int i = 1; i < N - 1; i++) {
		sum[i] = sum[i - 1] + mu[i];
	}
}

void solve(int a, int b, int d) {
	int ans = 0;
	a /= d; b/= d;
	int ed = min(a, b);
	for(int i = 1; i <= ed; ) {
		int r1 = a / (a / i);
		int r2 = b / (b / i);
		int r = min(r1, r2);
		ans += (sum[r] - sum[i - 1]) * (a / i) * (b / i);
		i = r + 1;
	}
	printf("%d\n", ans);
}

int main() {
	getMu();
	int t, a, b, d;
	scanf("%d", &t);
	while(t--) {
		scanf("%d%d%d", &a, &b, &d);
		solve(a, b, d);
	}
	return 0;
}

2.$YY的GCD$

题目描述

神犇 $YY$ 虐完数论后给某某出了一题

给定 $N,M$，求 $1\le x\le N,1\le y\le M$ 且 $gcd(x,y)$ 为质数的 $(x,y)$ 有多少对。

输入格式

第一行一个整数 $T$ 表述数据组数。

接下来 $T$ 行，每行两个正整数 $N,M$。

输出格式

$T$ 行，每行一个整数表示第 $i$ 组数据的结果。

输入输出样例

输入#1

2
10 10
100 100

输出#1

30
2791

说明/提示

$T=10^4$, $N,M\le 10^7$

解题思路：
$$\begin{gathered} \sum _{p\in prime}\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==p] \\ =\sum _{p\in prime}\sum _{i=1}^n\sum _{j=1}^m[gcd(i/p,j/p)==1] \\ =\sum _{p\in prime}\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{p}\rfloor }[gcd(i,j)==1] \\ =\sum _{p\in prime}\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{p}\rfloor }\sum _{d|gcd(i,j)}\mu (d) \\ =\sum _{p\in prime}\sum _{d=1}^{min(\lfloor \frac{n}{p}\rfloor ,\lfloor \frac{m}{p}\rfloor )}\mu (d)\ast \lfloor \frac{n}{pd}\rfloor \ast \lfloor \frac{m}{pd}\rfloor \end{gathered}$$
如果只是化到这里，就算用数论分块也是过不了，也会t掉

仔细看下就会发现，可以令 $T=pd$,原式就可以变为
$$\sum _{T=1}^{min(n,m)}\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor \sum _{p\in prime|T}\mu (\frac{T}{p})$$
可以看出来后面的那个和式是可以预处理的

链接：https://www.luogu.com.cn/problem/P3455

AC代码：

#include<cstdio>
#include<algorithm>

using namespace std;

const int N = 1e7 + 10;
const int maxx = 1e7 + 1;

bool vis[N];
int mu[N], prime[N], cnt, sum[N];

void getMu() {
	mu[1] = 1;
	for(int i = 2; i < maxx; i++) {
		if(!vis[i]) {
			prime[cnt++] = i;
			mu[i] = -1;
		}
		for(int j = 0; j < cnt && i * prime[j] < maxx; j++) {
			vis[i * prime[j]] = 1;
			mu[i * prime[j]] = -mu[i];
			if(i % prime[j] == 0) {
				mu[i * prime[j]] = 0;
				break;
			}
		} 
	}
	
	for(int i = 0; i < cnt; i++) {
		for(int j = 1; j * prime[i] < maxx; j++) {
			sum[prime[i] * j] += mu[j];
		}
	}
	
	for(int i = 1; i < maxx; i++) {
		sum[i] += sum[i - 1];
	}
	
}

void solve(int n, int m) {
	long long ans = 0;
	for(int i = 1; i <= min(n, m); ) {
		int r = min(n / (n / i), m / (m / i));
		ans += 	1ll * (sum[r] - sum[i - 1]) * (n / i) * (m / i);
		i = r + 1;
	} 
	printf("%lld\n", ans);
}

int main() {
	getMu();
	int n, m;
	int t;
	scanf("%d", &t);
	while(t--) {
		scanf("%d%d", &n, &m);
		solve(n, m);
	}
	return 0;
} 

3.最小公倍数之和

题目描述

对于 $A_1,A_2,...,A_N$,求

${\sum }_{i=1}^N{\sum }_{j=1}^{N}lcm(A_i,A_j)$的值。

$lcm(a,b)$ 表示 $a$ 和 $b$ 的最小公倍数

输入格式

第一行，一个整数 $N$。

第二行， $N$ 个整数 $A_1,A_2,...,A_N$。

输出格式

一个整数，表示所求的值。

输入输出样例

输入#1

2
2 3

输出#1

17

说明/提示

对于 $30\%$ 的数据, $1\le N\le 1000$; $1\le A_i\le 50000$;

对于另外 $30\%$ 的数据, $1\le N\le 50000$; $1\le A_i\le 1000$;

对于 $100\%$ 的数据, $1\le N\le 50000$; $1\le A_i\le 50000$。

解题思路：
$$\begin{gathered} \sum _{i=1}^N\sum _{j=1}^{N}lcm(A_i,A_j) \\ =\sum _{i=1}^N\sum _{j=1}^N\frac{A_i\ast A_j}{gcd(A_i,A_j)} \\ =\sum _{k=1}^N\sum _{i=1}^N\sum _{j=1}^N\frac{A_i\ast A_j}{k}[gcd(A_i,A_j)==k] \end{gathered}$$
写到这发现如果艾弗森表达式写成 $[gcd(A_i/k,A_j/k)==1]$，前面的 $i,j$范围没办法缩，因为 $A_i\ne i$,

所以我们转化下思路，怎么让 $i$ 表示 $A_i$呢？如果我们记录 $A_i$ 的最大值 $n$，然后再用一个数组，比如说

$C_i$来记录每一个 $A_i$ 出现的次数，然后就可以写成：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{n}lcm(i,j)\ast c_i\ast c_j \\ =\sum _{i=1}^n\sum _{j=1}^n\frac{i\ast j}{gcd(i,j)}\ast c_i\ast c_j \\ =\sum _{k=1}^n\sum _{i=1}^n\sum _{j=1}^n\frac{i\ast j}{k}[gcd(i,j)==k]\ast c_i\ast c_j \\ =\sum _{k=1}^n\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{k}\rfloor }i\ast j\ast k[gcd(i,j)==1]\ast c_{ik}\ast c_{jk} \\ =\sum _{k=1}^n\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{k}\rfloor }i\ast j\ast k\sum _{d|gcd(i,j)}\mu (d)\ast c_{ik}\ast c_{jk} \\ =\sum _{k=1}^n\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{i=1}^{\lfloor \frac{n}{kd}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{kd}\rfloor }i\ast j\ast k\ast d^2\ast c_{ikd}\ast c_{jkd}\ast \mu (d) \\ =\sum _{T=1}^n\sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{T}\rfloor }i\ast j\ast T\ast c_{iT}\ast c_{jT}\sum _{d|T}\mu (d)\ast d \\ =\sum _{T=1}^{n}T\ast (\sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }i\ast c_{iT}{)}^2\sum _{d|T}\mu (d)\ast d \end{gathered}$$
链接：https://www.luogu.com.cn/problem/P3911

AC代码：

#include<cstdio>
#include<algorithm>

using namespace std;

const int N = 1e5 + 10;
bool vis[N];
int prime[N], cnt, mu[N], sum[N], x, maxx = 0, s[N];

void getMu() {
	mu[1] = 1;
	for(int i = 2; i < N; i++) {
		if(!vis[i]) {
			prime[cnt++] = i;
			mu[i] = -1;
		}
		for(int j = 0; j < cnt && i * prime[j] < N; j++) {
			vis[i * prime[j]] = 1;
			mu[i * prime[j]] = -mu[i];
			if(i % prime[j] == 0) {
				mu[i * prime[j]] = 0;
				break;
			}
		}
	}
	
	for(int i = 1; i < N; i++) {
		for(int j = 1; j * i < N; j++) {
			sum[i * j] += mu[i] * i;
		}
	}
	
}

void solve() {
	long long ans = 0;
	for(int T = 1; T <= maxx; T++) {
		long long sumt = 0;
		for(int i = 1; i <= maxx / T; i++) {
			sumt += 1ll * i * s[i * T];
		}
		ans += 1ll * T * sumt * sumt * sum[T];
	}
	printf("%lld\n", ans);
}
.
int main() {
	getMu();
	int n;
	scanf("%d", &n);
	for(int i = 0; i < n; i++) {
		scanf("%d", &x);
		maxx = max(maxx, x);
		s[x]++;
	}
	solve();
	return 0;
}

4.能量采集

题目描述

栋栋有一块长方形的地，他在地上种了一种能量植物，这种植物可以采集太阳光的能量。在这些植物采集能量后，栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起。

栋栋的植物种得非常整齐，一共有 $n$ 列，每列有 $m$ 棵，植物的横竖间距都一样，因此对于每一棵植物，栋栋可以用一个坐标 $(x,y)$ 来表示，其中 xx 的范围是 $1$ 至 $n$ ，$y$ 的范围是 $1$ 至 $m$ ，表示是在第 $x$ 列的第 $y$ 棵。

由于能量汇集机器较大，不便移动，栋栋将它放在了一个角上，坐标正好是 $(0,0)$。

能量汇集机器在汇集的过程中有一定的能量损失。如果一棵植物与能量汇集机器连接而成的线段上有 $k$ 棵植物，则能量的损失为 $2k+1$。例如，当能量汇集机器收集坐标为 $(2,4)$ 的植物时，由于连接线段上存在一棵植物 $(1,2)$ ，会产生 $3$ 的能量损失。注意，如果一棵植物与能量汇集机器连接的线段上没有植物，则能量损失为 $1$ 。现在要计算总的能量损失。

下面给出了一个能量采集的例子，其中 $n=5$ ，$m=4$，一共有 $20$ 棵植物，在每棵植物上标明了能量汇集机器收集它的能量时产生的能量损失。

在这个例子中，总共产生了 $36$ 的能量损失。

输入格式

一行两个整数 $n,m$。

输出格式

仅包含一个整数，表示总共产生的能量损失。

输入输出样例

输入#1

5 4

输出 #1

3 6

输入#2

3 4

输出#2

20

说明/提示

$1\le n,m\le 10^5$

解题思路：

不难发现从 $(0,0)$ 到 $(x,y)$ 的一条线上不算上这两个点一共有 $gcd(x,y)-1$ 个点

所以可以列个式子：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^m(2\ast gcd(i,j)-1) \\ =2\ast (\sum _{i=1}^n\sum _{j=1}^{m}gcd(i,j))-n\ast m \end{gathered}$$
所以这题就相当于求解这个和式
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{m}gcd(i,j) \\ =\sum _{k=1}^{min(n,m)}\sum _{i=1}^n\sum _{j=1}^{m}k[gcd(i,j)==k] \\ =\sum _{k=1}^{min(n,m)}\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }k[gcd(i,j)==1] \\ =\sum _{k=1}^{min(n,m)}\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }k\sum _{d|gcd(i,j)}\mu (d) \\ =\sum _{k=1}^{min(n,m)}\sum _{d=1}^{min(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor )}k\ast \mu (d)\ast \lfloor \frac{n}{kd}\rfloor \ast \lfloor \frac{m}{kd}\rfloor \\ =\sum _{T=1}^{min(n,m)}\lfloor \frac{n}{T}\rfloor \ast \lfloor \frac{m}{T}\rfloor \ast \sum _{k|T}\mu (\frac{T}{k})\ast k \end{gathered}$$
链接：https://www.luogu.com.cn/problem/P1447

AC代码：

#include<cstdio>
#include<algorithm>

using namespace std;

const int N = 1e6 + 10;
bool vis[N];
int prime[N], cnt, mu[N], sum[N];

void getMu() {
	mu[1] = 1;
	for(int i = 2; i < N; i++) {
		if(!vis[i]) {
			prime[cnt++] = i;
			mu[i] = -1;
		}
		for(int j = 0; j < cnt && i * prime[j] < N; j++) {
			vis[i * prime[j]] = 1;
			mu[i * prime[j]] = -mu[i];
			if(i % prime[j] == 0) {
				mu[i * prime[j]] = 0;
				break;
			}
		}
	}
	
	for(int i = 1; i < N; i++) {
		for(int j = 1; j * i < N; j++) {
			sum[i * j] += mu[j] * i;
		}
	}
	
	for(int i = 1; i < N; i++) {
		sum[i] += sum[i - 1];
	}
}

void solve(int n, int m) {
	long long ans = 0;
	for(int i = 1; i <= min(n, m); ) {
		int r = min(n / (n / i), m / (m / i));
		ans += 1ll * (sum[r] - sum[i - 1]) * (n / i) * (m / i);
		i = r + 1;
	}
	printf("%lld\n", 2 * ans - 1ll * n * m);
}

int main() {
	getMu();
	int n, m;
	scanf("%d%d", &n, &m);
	solve(n, m);
	return 0;
} 

这样处理后面那个和式是可以过的，但是这题还可以再优化一下

首先 ${\sum }_{d|m}\phi (d)=m\Leftarrow \Rightarrow \phi (m)={\sum }_{d|m}\mu (d)\frac{m}{d}$

然后发现后面那个要预处理的和式转换一下就等于 $\phi (T)$ ，所以只要在素数筛里求 $\phi (T)$就可以

再补充一下欧拉函数的知识($p$ 代表素数)

如果 $n=p_1^{k_1}\ast p_2^{k_2}\ast ...\ast p_m^{k_m}$ ,那么 $\phi (n)=\phi (p_1^{k_1})\ast \phi (p_2^{k_2})\ast ...\ast \phi (p_m^{k_m})$, 因为欧拉函数是积性函数

还有 $\phi (p^k)=p^k-p^{k-1}$

所以
$$\begin{gathered} \phi (n)=p_1^{k_1}\ast p_2^{k_2}\ast ...\ast p_m^{k_m}\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast ...\ast (1-\frac{1}{p_m}) \\ =n\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast ...\ast (1-\frac{1}{p_m}) \end{gathered}$$
所以如果 $p$ 是素数，并且 $p$ 是 $x$ 的因子
$$\phi (p\ast x)=p\ast x\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast ...\ast (1-\frac{1}{p_m})=\phi (x)\ast p$$
如果 $p$ 不是 $x$ 的因子, 那么 $p$ 和 $x$ 互素，则
$$\phi (p\ast x)=\phi (p)\ast \phi (x)$$
AC代码：

#include<cstdio>
#include<algorithm>

using namespace std;
typedef long long ll;
const int N = 1e6 + 10;
bool vis[N];
int prime[N], cnt, mu[N];
ll sum[N], phi[N];

void getMu() {
	phi[1] = sum[1] = 1;
	for(int i = 2; i < N; i++) {
		if(!vis[i]) {
			prime[cnt++] = i;
			phi[i] = i - 1;
		}
		for(int j = 0; j < cnt && i * prime[j] < N; j++) {
			vis[i * prime[j]] = 1;
			if(i % prime[j] == 0) {
				phi[prime[j] * i] = phi[i] * prime[j];
				break;
			}
			phi[prime[j] * i] = phi[prime[j]] * phi[i];
		}
		sum[i] = sum[i - 1] + phi[i];
	}
}

void solve(int n, int m) {
	ll ans = 0;
	for(int i = 1; i <= min(n, m); ) {
		int r = min(n / (n / i), m / (m / i));
		ans += 1ll * (sum[r] - sum[i - 1]) * (n / i) * (m / i);
		i = r + 1;
	}
	printf("%lld\n", 2 * ans - 1ll * n * m);
}

int main() {
	getMu();
	int n, m;
	scanf("%d%d", &n, &m);
	solve(n, m);
	return 0;
}