【数学】莫比乌斯反演（以 P2522 和 P3327 为例）

一、莫比乌斯函数
$$\mu (n)=\{\begin{array}{ll}1& n=1\\ 0& n\text{含有平方因子}\\ (-1{)}^k& k\text{为}n\text{的本质不同质因子个数}\end{array}$$
性质：
$$\sum _{d|n}\mu (d)=\{\begin{array}{ll}1& n=1\\ 0& n\ne 1\end{array}$$
证明：

$\forall n\in \mathbb{Z}$，通过唯一分解定理得：$n=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}\dots p_k^{{\alpha }_k}$

$\forall d|n$，通过唯一分解定理得：$d=p_1^{{\beta }_1}{p}_2^{{\beta }_2}\dots p_k^{{\beta }_k}$，${\beta }_i\in [0,{\alpha }_i]$

若 $\exists i,{\beta }_i\ge 2$，那么由莫比乌斯函数的定义得，这个 $d$ 对于和没有贡献。

反之，$\forall i,{\beta }_i\le 1$，记 $m$ 为 ${\beta }_i=1$ 的 $i$ 的个数，那么 ${\sum }_{d|n}\mu (d)={\sum }_{d|n}(-1{)}^m={\sum }_{i=0}^k{C}_k^i(-1{)}^i$

解释：从 $k$ 个数中选出 $i$ 个 ${\beta }_j=1,j\in [1,k]$，方案数为 $C_k^i$，贡献为 $(-1{)}^i$。

观察发现为二项式定理，故：${\sum }_{i=0}^k{C}_k^i(-1{)}^i=(1+(-1){)}^k=0$

所以，只有 $n=1$ 的时候，${\sum }_{d|n}\mu (d)=1$，证毕。

二、莫比乌斯反演

若有 $F(n)={\sum }_{d|n}f(d)$，那么有 $f(n)={\sum }_{d|n}\mu (d)F(\frac{n}{d})$

若有 $F(n)={\sum }_{n|d}f(d)$，那么有 $f(n)={\sum }_{n|d}\mu (\frac{d}{n})F(d)$

其中实际应用通常使用②式，且 $F$ 好求，而 $f$ 不好求，但是满足条件，所以可以通过 $F$ 来求解 $f$。

一式证明：
$$\begin{array}{rlrl}& & & \sum _{d|n}\mu (d)\sum _{i|\frac{n}{d}}f(i)\\ & =& & \sum _{d|n}f(i)\sum _{d|\frac{n}{i}}\mu (d)\to \sum _{d|\frac{n}{i}}\mu (d)=1当且仅当i为n\\ & =& & f(n)\end{array}$$
二式证明：
$$\begin{array}{rl}\sum _{n|d}\mu (\frac{d}{n})F(d)=& \sum _{n|d}\mu (\frac{n}{d})\sum _{d|i}f(i)\\ =& \sum _{n|i}f(i)\sum _{d^{\prime }|\frac{i}{n}}\mu (d^{\prime })\to d^{\prime }=\frac{n}{d}\\ =& f(n)\end{array}$$
三、例题精讲

P2522 [HAOI2011] Problem b

对于给出的 $n$ 个询问，每次求有多少个数对 $(x,y)$，满足 $a\le x\le b$，$c\le y\le d$，且 $\gcd (x,y)=k$，$\gcd (x,y)$ 函数为 $x$ 和 $y$ 的最大公约数。

$1\le n,k\le 5\times 10^4$，$1\le a\le b\le 5\times 10^4$，$1\le c\le d\le 5\times 10^4$。

---

将数对 $(x,y)$ 映射成平面直角坐标系中的一个点，那么对于每一次查询相当于在左下角为 $(a,c)$，右上角为 $(b,d)$ 的矩形内查询有多少个点 $(x,y)$ 满足 $\gcd (x,y)=k$。

那么，通过差分可以将问题简化为求左下角 $(0,0)$，右上角为 $(a,b)$ 的矩形内查询有多少个点 $(x,y)$ 满足 $\gcd (x,y)=k$。

考虑如何设计 $F(n)$ 与 $f(n)$，由于为了 $F$ 好求，所以设 $F(n)={\sum }_{x,y}[n|\gcd (x,y)]$，由于要满足约数和的关系，所以设 $f(n)={\sum }_{x,y}[\gcd (x,y)=n]$。

那么，观察可以发现，$F(n)={\sum }_{n|d}f(d)$。所以，可以套用莫比乌斯反演得：$f(n)={\sum }_{n|d}\mu (\frac{n}{d})F(d)$

那么，考虑如何求解 $F(n)$：
$$\begin{array}{rl}F(n)=& \sum _{x=1}^a\sum _{y=1}^b[n|\gcd (x,y)]\\ =& \sum _{x=1}^a\sum _{y=1}^b[n|x,n|y]\\ =& \lfloor \frac{a}{n}\rfloor \lfloor \frac{b}{n}\rfloor \end{array}$$
则：$f(n)={\sum }_{n|d}\mu (\frac{n}{d})\lfloor \frac{a}{d}\rfloor \lfloor \frac{b}{d}\rfloor$

继续化简方便求解（尽量化成区间的和）：
$$\begin{array}{rl}f(n)=& \sum _{n|d}\mu (\frac{n}{d})\lfloor \frac{a}{d}\rfloor \lfloor \frac{b}{d}\rfloor \\ =& \sum _{d^{\prime }\in \mathbb{Z}}\mu (d^{\prime })\lfloor \frac{a}{d^{\prime }n}\rfloor \lfloor \frac{b}{d^{\prime }n}\rfloor \to d^{\prime }=\frac{n}{d}\\ =& \sum _{d^{\prime }\in \mathbb{Z}}\mu (d^{\prime })\lfloor \frac{a^{\prime }}{d^{\prime }}\rfloor \lfloor \frac{b^{\prime }}{d^{\prime }}\rfloor \to a^{\prime }=\lfloor \frac{a}{n}\rfloor ,b^{\prime }=\lfloor \frac{b}{n}\rfloor \end{array}$$
疑难解答

为什么 $\lfloor \frac{a}{d^{\prime }n}\rfloor =\lfloor \frac{a^{\prime }}{d^{\prime }}\rfloor =\lfloor \frac{\lfloor \frac{a}{n}\rfloor }{d^{\prime }}\rfloor$？

• $\forall a,b,c\in \mathbb{Z},\lfloor \frac{a}{bc}\rfloor =\lfloor \frac{\lfloor \frac{a}{b}\rfloor }{c}\rfloor$

  证明：$\frac{a}{b}=\lfloor \frac{a}{b}\rfloor +r$，$r\in [0,1)$

  $\lfloor \frac{a}{bc}\rfloor =\lfloor \frac{a}{b}\times \frac{1}{c}\rfloor =\lfloor \frac{1}{c}(\lfloor \frac{a}{b}\rfloor +r)\rfloor =\lfloor \frac{\lfloor \frac{a}{b}\rfloor }{c}+\frac{r}{c}\rfloor =\lfloor \frac{\lfloor \frac{a}{b}\rfloor }{c}\rfloor$

最后，通过整除分块计算最终的式子即可，$\mu$ 函数预处理前缀和。查询的答案就是 $f(k)$

---

#include <bits/stdc++.h>
#define int long long

using namespace std;

typedef pair<int, int> PII;
typedef long long LL;

const int N = 5e4 + 10;

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

void prework() { //筛法求莫比乌斯函数
	mu[1] = 1;
	for (int i = 2; i < N; i ++) {
		if (!vis[i]) prime[ ++ cnt] = i, mu[i] = -1;
		for (int j = 1; prime[j] * i < N; j ++) {
			vis[prime[j] * i] = 1;
			if (i % prime[j] == 0) break;
			mu[prime[j] * i] = -mu[i];
		}
	}
	for (int i = 1; i < N; i ++) //前缀和
		mu[i] += mu[i - 1];
}

int f(int a, int b, int k) {
	a = a / k, b = b / k;
	int res = 0;
	for (int l = 1, r; l <= min(a, b); l = r + 1) { //整除分块
		r = min(a / (a / l), b / (b / l));
		res += (mu[r] - mu[l - 1]) * (a / l) * (b / l);
	}
	return res;
}

void solve() {
	int a, b, c, d, k;
	cin >> a >> b >> c >> d >> k;

	cout << f(b, d, k) - f(a - 1, d, k) - f(b, c - 1, k) + f(a - 1, c - 1, k) << endl; //差分
}

signed main() {
	cin.tie(0);
	cout.tie(0);
	ios::sync_with_stdio(0);

	prework();

	int dt;
	cin >> dt;

	while (dt --)
		solve();

	return 0;
}

---

P3327 [SDOI2015] 约数个数和

设 $d(x)$ 为 $x$ 的约数个数，给定 $n,m$，求
$$\sum _{i=1}^n\sum _{j=1}^{m}d(ij)$$

$1\le T,n,m\le 50000$

---

$$\mathrm{d}(ij)=\sum _{x|i}\sum _{y|i}[\gcd (x,y)=1]$$

注：$\mathrm{d}(i)$ 表示 $i$ 的约数个数

记 $\mathrm{res}$ 为答案，则 $\mathrm{res}={\sum }_{i=1}^n{\sum }_{j=1}^m{\sum }_{x|i}{\sum }_{y|j}[\gcd (x,y)=1]$

观察发现后半部分与上题惊人的类似，所以不难想到可以设：
$$\begin{array}{rl}& F(n)=\sum _{i=1}^N\sum _{j=1}^M\sum _{x|i}\sum _{y|j}[n|\gcd (x,y)]\\ & f(n)=\sum _{i=1}^N\sum _{j=1}^M\sum _{x|i}\sum _{y|j}[\gcd (x,y)=n]\end{array}$$
则有：$F(n)={\sum }_{n|d}f(d)$

故：$f(n)={\sum }_{n|d}\mu (\frac{n}{d})F(d)$

继续化简：
$$\begin{array}{rl}f(n)=& \sum _{n|d}\mu (\frac{n}{d})F(d)\\ =& \sum _{n|d}\mu (\frac{n}{d})\sum _{i=1}^N\sum _{j=1}^M\sum _{x|i}\sum _{y|j}[n|\gcd (x,y)]\\ =& \sum _{n|d}\mu (\frac{n}{d})\sum _{i=1}^N\sum _{j=1}^M\sum _{x|i}\sum _{y|j}[n|x,n|y]\\ =& \sum _{n|d}\mu (\frac{n}{d})\sum _{x=1}^N\sum _{y=1}^M\sum _{x|i}\sum _{y|j}[n|x,n|y]\\ =& \sum _{n|d}\mu (\frac{n}{d})\sum _{x=1}^N\sum _{y=1}^M\lfloor \frac{N}{x}\rfloor \lfloor \frac{M}{y}\rfloor [n|x,n|y]\to 观察到x,y必须都是n的倍数，所以引出下面的x^{\prime },y^{\prime }\\ =& \sum _{n|d}\mu (\frac{n}{d})\sum _{x^{\prime }=1}^{\frac{N}{n}}\sum _{y^{\prime }=1}^{\frac{M}{n}}\lfloor \frac{N}{x^{\prime }n}\rfloor \lfloor \frac{M}{y^{\prime }n}\rfloor \to x^{\prime }=\lfloor \frac{x}{n}\rfloor ,y^{\prime }=\lfloor \frac{y}{n}\rfloor \\ =& \sum _{n|d}\mu (\frac{n}{d})\sum _{x^{\prime }=1}^{N^{\prime }}\sum _{y^{\prime }=1}^{M^{\prime }}\lfloor \frac{N^{\prime }}{x^{\prime }}\rfloor \lfloor \frac{M^{\prime }}{y^{\prime }}\rfloor \to N^{\prime }=\lfloor \frac{N}{n}\rfloor ,M^{\prime }=\lfloor \frac{M}{n}\rfloor \end{array}$$
经过一番推导，可以发现后面的式子可以用整除分块做了，即：
$$\begin{array}{rl}& \sum _{x^{\prime }=1}^{N^{\prime }}\sum _{y^{\prime }=1}^{M^{\prime }}\lfloor \frac{N^{\prime }}{x^{\prime }}\rfloor \lfloor \frac{M^{\prime }}{y^{\prime }}\rfloor \\ =& \sum _{x^{\prime }=1}^{N^{\prime }}\lfloor \frac{N^{\prime }}{x^{\prime }}\rfloor \sum _{y^{\prime }=1}^{M^{\prime }}\lfloor \frac{M^{\prime }}{y^{\prime }}\rfloor \end{array}$$
那么，令 $g(n)={\sum }_{i=1}^n\lfloor \frac{n}{i}\rfloor$，则 ${\sum }_{x^{\prime }=1}^{N^{\prime }}\lfloor \frac{N^{\prime }}{x^{\prime }}\rfloor {\sum }_{y^{\prime }=1}^{M^{\prime }}\lfloor \frac{M^{\prime }}{y^{\prime }}\rfloor =g(N^{\prime })g(M^{\prime })$

代入会原式得：
$$\begin{array}{rl}f(n)=& \sum _{n|d}\mu (\frac{n}{d})F(d)\\ =& \sum _{n|d}\mu (\frac{n}{d})g(\lfloor \frac{N}{d}\rfloor )g(\lfloor \frac{M}{d}\rfloor )\\ =& \sum _{d^{\prime }=1}\mu (d^{\prime })g(\lfloor \frac{N}{n{d}^{\prime }}\rfloor )g(\lfloor \frac{M}{n{d}^{\prime }}\rfloor )\to d^{\prime }=\frac{d}{n}\\ =& \sum _{d^{\prime }=1}\mu (d^{\prime })g(\lfloor \frac{N^{\prime }}{d^{\prime }}\rfloor )g(\lfloor \frac{M^{\prime }}{d^{\prime }}\rfloor )\to N^{\prime }=\lfloor \frac{N}{n}\rfloor ,M^{\prime }=\lfloor \frac{M}{n}\rfloor \end{array}$$
这时候，发现可以在此进行整除分块，这样求一个 $f(n)$ 的时间复杂度为 $\sqrt{n}$，则有 $q$ 次查询，总时间复杂度为 $O(q\sqrt{n})$

最终的答案为：
$$f(1)=\sum _{d=1}\mu (d^{\prime })g(\lfloor \frac{N}{d}\rfloor )g(\lfloor \frac{M}{d}\rfloor )$$

---

#include <bits/stdc++.h>

using namespace std;

typedef pair<int, int> PII;
typedef long long LL;

const int N = 5e4 + 10;

int n, m;
int mu[N], prime[N], cnt, vis[N];
int g[N];

void prework() {
	mu[1] = 1;
	for (int i = 2; i < N; i ++) {
		if (!vis[i]) prime[ ++ cnt] = i, mu[i] = -1;
		for (int j = 1; prime[j] * i < N; j ++) {
			vis[i * prime[j]] = 1;
			if (i % prime[j] == 0) break;
			mu[i * prime[j]] = -mu[i];
		}
	}
	for (int i = 1; i < N; i ++)
		mu[i] += mu[i - 1];
	
	for (int i = 1; i < N; i ++) {
		for (int l = 1, r; l <= i; l = r + 1) {
			r = min(i, i / (i / l));
			g[i] += (r - l + 1) * (i / l);
		}
	}
}

void solve() {
	cin >> n >> m;

	LL res = 0;
	for (int l = 1, r; l <= min(n, m); l = r + 1) {
		r = min(min(n / (n / l), m / (m / l)), min(n, m));
		res += (LL)(mu[r] - mu[l - 1]) * g[n / l] * g[m / l];
	}

	cout << res << '\n';
}

signed main() {
	cin.tie(0);
	cout.tie(0);
	ios::sync_with_stdio(0);

	prework();
	int dt;
	cin >> dt;

	while (dt --)
		solve();

	return 0;
}