浅谈莫比乌斯反演

浅谈莫比乌斯反演

 

简介（$\text{cheˇ daˋn}$）：

莫比乌斯反演是数论中的重要内容。对于一些函数 ，如果很难直接求出它的值，而容易求出其倍数和或约数和 ，那么可以通过莫比乌斯反演简化运算，进行求值。
 

前置芝士：

引理$1$：$\forall a,b,c\in Z,\lfloor \frac{a}{bc}\rfloor =\lfloor \frac{\lfloor \frac{a}{b}\rfloor }{c}\rfloor$

引理$2$： ∀ n ∈ N + , ∣ { ⌊ n d ⌋ ∣ d ∈ N + , d ≤ n } ∣ ≤ ⌊ 2 n ⌋ \forall n \in N_+, \vert \{ \lfloor \frac{n}{d} \rfloor \vert d\in N_+,d\le n \}\vert \le \lfloor2\sqrt{n}\rfloor ∀n∈N+​,∣{⌊dn​⌋∣d∈N+​,d≤n}∣≤⌊2n ​⌋

Tips: ∣ { } ∣ \vert\{\}\vert ∣{}∣表示集合的长度，即集合中元素的个数。
 

1.整除分块：

用来求$\sum \lfloor \frac{n}{i}\rfloor$，直接模拟时间复杂度为$O(n)$，但是用整除分块可以优化到$O(\sqrt{n})$。
对于$\forall i\le n$ ，我们需要找到一个$i\le j_{max}\le n$ ，使得$\lfloor \frac{n}{i}\rfloor =\lfloor \frac{n}{j_{max}}\rfloor$。

很容易发现：$r=\lfloor \frac{n}{\lfloor \frac{n}{i}\rfloor }\rfloor$。（证明过程详见Oi-wiki）
 

例题：

P2261 [CQOI2007]余数求和

题目描述：
给出正整数 $n$ 和 $k$，请计算$G(n,k)=\sum _{i=1}^{n}k\mathrm{m}\mathrm{o}\mathrm{d}i$
其中 $k\mathrm{m}\mathrm{o}\mathrm{d}i$ 表示 $k$ 除以 $i$ 的余数。

$\mathtt{S}\mathtt{o}\mathtt{l}\mathtt{u}\mathtt{t}\mathtt{i}\mathtt{o}\mathtt{n}:$
一道非常基础的模板题，我们先来化简。
$Ans=\sum _{i=1}^{n}k\mathrm{m}\mathrm{o}\mathrm{d}i$
   $=\sum _{i=1}^{n}k-\lfloor \frac{k}{i}\rfloor \times i$
   $=n\times k-\sum _{i=1}^n\lfloor \frac{k}{i}\rfloor \times i$
然后就直接用整除分块做就行了
$\mathtt{C}\mathtt{o}\mathtt{d}\mathtt{e}:$

#include <bits/stdc++.h>

using namespace std;

long long n, k, ans;

int main() {
	scanf("%lld%lld", &n, &k);
	ans = n * k;
	for (int l = 1, r; l <= n; l = r + 1) {
		if (k / l) r = min(k / (k / l), n);
		else r = n;
		ans -= (k / l) * (r - l + 1) * (l + r) / 2;
	}
	printf("%lld", ans);
	return 0;
}

 

2.狄利克雷卷积与数论函数

定义：

• 数论函数： 值域为整数的函数。
• 狄利克雷卷积： $(f\ast g)(n)=\sum _{d|n}f(d)g(\frac{n}{d})$或者$(f\ast g)(n)=\sum _{x\times y=n}f(x)g(y)$

狄利克雷卷积的性质：

• 交换律： $f\ast g=g\ast f$
• 结合律： $(f\ast g)\ast h=f\ast (g\ast h)$
• 分配律： $f\ast (g+h)=f\ast g+f\ast h$
• $f\ast \epsilon =f$
   

积性函数：

定义：

• 对于$gcd(a,b)=1$，有$\{\begin{array}{l}f(1)=1\\ f(ab)=f(a)\times f(b)\end{array}$，则$f(x)$被称为积性函数。
   
• 对于$\forall a,b\in C$，有$\{\begin{array}{l}f(1)=1\\ f(ab)=f(a)\times f(b)\end{array}$，则$f(x)$被称为完全积性函数。
   

结论：

• 两个积性函数的狄利克雷卷积还是积性函数。
• 积性函数的逆还是积性函数。

常见的积性函数：

• 恒等函数：$\mathrm{I}(n)\equiv 1$（无论$n$是啥，恒等于1，废柴一个 ）
• 元函数：$\epsilon (n)=[n=1]$
• 单位函数：$id(n)=n$
• 莫比乌斯函数：$\mu (x)=\{\begin{array}{ll}0& \exists d>1且d^2|x\\ (-1{)}^k& x=\prod _{i=1}^k{p}_i(p_i\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e})\end{array}$
  莫比乌斯函数的性质：$\mu \ast \mathrm{I}=\sum _{d|n}\mu (d)=[n=1]=\epsilon$
  莫比乌斯函数的模板如下：

  void Mobius() {
      Mu[1] = 1;
      for (int i = 2; i <= Maxn; i++) {
   	    if (!vis[i]) Mu[i] = -1, p[++cnt] = i;
   	    for (int j = 1; j <= cnt && (int)i * p[j] <= Maxn; j++) {
   		    vis[i * p[j]] = true;
   		    if (i % p[j] == 0) break;
   		    Mu[i * p[j]] = -Mu[i];
   	    }
   	Mu[i] += Mu[i - 1];
      }
  }	
  

 

然后进入正题：

3.莫比乌斯反演

设两个数论函数$F(x),f(x)$，则有：
$$\{\begin{array}{ll}F(n)=\sum _{d|n}f(d)& ①\\ f(n)=\sum _{d|n}\mu (d)F(\frac{n}{d})& ②\\ f(n)=\sum _{n|d}\mu (\frac{d}{n})F(d)& ③\end{array}$$
证：
由①可得：$$F\ast \mu =f\ast \mathrm{I}\ast \mu$$
由莫比乌斯函数的性质可得：$$F\ast \mu =f\ast \epsilon =f$$
所以，$$f=\mu \ast F=\sum _{d|n}\mu (d)F(\frac{n}{d})$$
 

例题：

P3455 [POI2007]ZAP-Queries

题目描述：

密码学家正在尝试破解一种叫 $BSA$ 的密码。
他发现，在破解一条消息的同时，他还需要回答这样一种问题：
给出 $a,b,d$求满足 $1\le x\le a，1\le y\le b$，且 $\gcd (x,y)=d$ 的二元组 $(x,y)$ 的数量。
因为要解决的问题实在太多了，他便过来寻求你的帮助。

 
输入格式：

输入第一行一个整数 $n$，代表要回答的问题个数。
接下来 $n$ 行，每行三个整数 $a,b,d$。

 
输出格式：

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

 

 
$\mathtt{S}\mathtt{o}\mathtt{l}\mathtt{u}\mathtt{t}\mathtt{i}\mathtt{o}\mathtt{n}:$
首先我们先把题意转换成人能看懂 的样子。
题目让我们求的东西：

$$\sum _{x=1}^a\sum _{y=1}^b[\gcd (x,y)=d]$$

我们先把$d$提出来：

$$\sum _{x=1}^{\frac{a}{d}}\sum _{y=1}^{\frac{b}{d}}[\gcd (x,y)=1]$$

很容易发现：

$$[\gcd (x,y)=1]=\sum _{k|\gcd (x,y)}\mu (k)$$

然后，开始化简：

$$\sum _{x=1}^{\frac{a}{d}}\sum _{y=1}^{\frac{b}{d}}\sum _{k|\gcd (x,y)}\mu (k)$$

拆掉$\gcd$，枚举k

$$\sum _{x=1}^{\frac{a}{d}}\sum _{y=1}^{\frac{b}{d}}\sum _{k=1}^{\min (\frac{a}{d},\frac{b}{d})}\mu (k)[k|x][k|y]$$

$$\sum _{k=1}^{\min (\frac{a}{d},\frac{b}{d})}\mu (k)\sum _{x=1}^{\frac{a}{d}}[k|x]\sum _{y=1}^{\frac{b}{d}}[k|y]$$

$$\sum _{k=1}^{\min (\frac{a}{d},\frac{b}{d})}\mu (k)\lfloor \frac{a}{dk}\rfloor \lfloor \frac{b}{dk}\rfloor$$

很容易看出，最终结果就用莫比乌斯函数和数论分块计算就行了。
 
$\mathtt{C}\mathtt{o}\mathtt{d}\mathtt{e}:$

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

using namespace std; 

int N, M, D, cnt;
int Mu[50005], vis[50005], p[50005];

void Mobius() {
	Mu[1] = 1;
	for (int i = 2; i <= 50005; i++) {
		if (!vis[i]) Mu[i] = -1, p[++cnt] = i;
		for (int j = 1; j <= cnt && (int)i * p[j] <= 50005; j++) {
			vis[i * p[j]] = true;
			if (i % p[j] == 0) break;
			Mu[i * p[j]] = -Mu[i];
		}
		Mu[i] += Mu[i - 1];
	}
}

void Solve() {
	int ans = 0;
	scanf("%lld%lld%lld", &N, &M, &D);
	if (N > M) swap(N, M);
	N /= D, M /= D;
	for (int l = 1, r; l <= N; l = r + 1) {
		r = min((M / (M / l)), (N / (N / l)));
		ans += (int)(Mu[r] - Mu[l - 1]) * (N / l) * (M / l);
	}
	printf("%lld\n", ans);
}

signed main() {
	Mobius();
	int T = 0;
	scanf("%lld", &T);
	while (T--) Solve();
	return 0;
}

 

P2257 YY的GCD

题目描述：

神犇 YY 虐完数论后给傻× kAc 出了一题
给定 $N,M$，求 $1\le x\le N$，$1\le y\le M$ 且 $\gcd (x,y)$ 为质数的 $(x,y)$ 有多少对。

 
输入格式：

第一行一个整数 $T$ 表述数据组数。
接下来 $T$ 行，每行两个正整数，$N,M$。

 
输出格式：

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

 

说明/提示：
$T=10^4，N,M\le 10^7$
 
$\mathtt{S}\mathtt{o}\mathtt{l}\mathtt{u}\mathtt{t}\mathtt{i}\mathtt{o}\mathtt{n}:$
这道题是前面那道题的升级版。
老规矩，我们先把题目化成人能看懂的样子：

$$\sum _{d\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}}\sum _{i=1}^x\sum _{j=1}^y[\gcd (x,y)=d]$$

根据套路，我们先把$d$提出来：

$$\sum _{d\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}}\sum _{i=1}^{\frac{x}{d}}\sum _{j=1}^{\frac{y}{d}}[\gcd (x,y)=1]$$

然后转换成莫比乌斯函数的形式：

$$\sum _{d\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}}\sum _{i=1}^{\frac{x}{d}}\sum _{j=1}^{\frac{y}{d}}\sum _{k|\gcd (x,y)}\mu (k)$$

把$\gcd$提出来，枚举$k$：

$$\sum _{d\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}}\sum _{i=1}^{\frac{x}{d}}\sum _{j=1}^{\frac{y}{d}}\sum _{k=1}^{\min (\frac{x}{d},\frac{y}{d})}\mu (k)[k|x][k|y]$$

各回各家：

$$\sum _{d\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}}\sum _{k=1}^{\min (\frac{x}{d},\frac{y}{d})}\mu (k)\sum _{i=1}^{\frac{x}{d}}[k|x]\sum _{j=1}^{\frac{y}{d}}[k|y]$$

后面的2个$\Sigma$可以直接算：

$$\sum _{d\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}}\sum _{k=1}^{\min (\frac{x}{d},\frac{y}{d})}\mu (k)\lfloor \frac{x}{dk}\rfloor \lfloor \frac{y}{dk}\rfloor$$

枚举$dk$为$T$：

$$\sum _{T=1}^{\min (x,y)}\lfloor \frac{x}{T}\rfloor \lfloor \frac{y}{T}\rfloor \sum _{t|T,t\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}}\mu (\frac{T}{t})$$

令

$$f(x)=\sum _{p|x,p\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}}\mu (\frac{x}{p})$$

设$x=i\times y$，$y$是$x$的最小质因子。
①当$x\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}$时，很容易得到：$f(x)=\mu (\frac{x}{x})=\mu (1)=1$
②$i\mathrm{m}\mathrm{o}\mathrm{d}y=0$时，即$x$存在多个最小质因子：

• 当$i$没有多个相同的质因子时，有且只有$p=y$时$\mu (\frac{x}{p})=\mu (i)\ne 0$。因为当$p\ne y$时，$\exists y>1且y^2|x$，
• 根据莫比乌斯函数的性质，很容易知道：这时的函数值为$0$，所以$f(x)=\mu (i)$。
• 当$i$有多个相同的质因子时，对于任意$p$，都有$\mu (\frac{x}{p})=\mu (i)=0$。
• 所以当$imody=0$时，$f(x)=\mu (i)$。

③$i\mathrm{m}\mathrm{o}\mathrm{d}y\ne 0$时，即$x$只存在一个最小质因子。

• $f(i)=\sum _{p\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e},p|i}\mu (\frac{i}{p}),f(x)=\sum _{p\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e},p|x}\mu (\frac{i\times y}{p})$
• 根据$\mu$线性筛的过程，很容易知道$\mu (\frac{i\times y}{p})=-\mu (\frac{i}{p})$，所以我们可以从$f(x)$中找到任意一个属于$f(i)$的项。
• 因此$f(x)$比$f(i)$多了一项$\mu (\frac{i\times y}{y})=\mu (i)$。
• 所以当$imody\ne 0$时，$f(x)=-f(i)+\mu (i)$。

综上所述，我们得出：
$$f(x)=\{\begin{array}{ll}\mu (1)=1,& case1:x\in \mathtt{p}\mathtt{r}\mathtt{i}\mathtt{m}\mathtt{e}\\ \mu (i),& case2:i\mathrm{m}\mathrm{o}\mathrm{d}y=0.\\ -f(i)+\mu (i),& case3:i\mathrm{m}\mathrm{o}\mathrm{d}y\ne 0.\end{array}$$
所以，$$Ans=\sum _{T=1}^{\min (x,y)}\lfloor \frac{x}{T}\rfloor \lfloor \frac{y}{T}\rfloor f(T)$$
我们只需要计算$f(T)$的前缀和与数论分块就行了。
$Code:$

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

using namespace std;

const int Maxn = 1e7 + 5, lnMaxn = Maxn / 16;
int N, M, T, cnt;
int f[Maxn], Mu[Maxn], p[lnMaxn];
bool vis[Maxn];

void Mobius() {
	Mu[1] = vis[1] = 1;
	for (int i = 2; i < Maxn; i++) {
		if (!vis[i]) Mu[i] = -1, f[i] = 1, p[++cnt] = i;   //case 1
		for (int j = 1; j <= cnt && i * p[j] < Maxn; j++) {
			vis[i * p[j]] = true;
			if (i % p[j] == 0) {
				f[i * p[j]] = Mu[i];  //case 2
				break;
			}
			else {
				f[i * p[j]] = -f[i] + Mu[i];  //case 3
			}
			Mu[i * p[j]] = -Mu[i];
		}
		f[i] += f[i - 1]; 
	}
}

void Solve() {
	int ans = 0;
	scanf("%llu%llu", &N, &M);
	if (N > M) swap(N, M);
	for (int l = 1, r; l <= N; l = r + 1) {
		r = min(N / (N / l), (M / (M / l)));
		ans += (f[r] - f[l - 1]) * (N / l) * (M / l);
	}
	printf("%llu\n", ans);
}

signed main() {
	scanf("%llu", &T);
	Mobius();
	while (T--) Solve();
	return 0;
} 

 

P3327 [SDOI2015]约数个数和

题目描述：

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

输入格式：

输入文件包含多组测试数据。
第一行，一个整数 $T$，表示测试数据的组数。
接下来的 $T$ 行，每行两个整数 $n,m$。

输出格式：

T 行，每行一个整数，表示你所求的答案。

说明/提示：
对于$100\%$的数据，$1\le T,n,m\le 50000$
 
$\mathtt{S}\mathtt{o}\mathtt{l}\mathtt{u}\mathtt{t}\mathtt{i}\mathtt{o}\mathtt{n}:$
首先我们得知道一个公式：$$d(ij)=\sum _{x|i}\sum _{y|j}[\gcd (x,y)=1]$$

然后进行化简：

$$\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}[\gcd (x,y)=1]$$

枚举$\gcd (x,y)$的因数：

$$\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}\sum _{d|\gcd (x,y)}\mu (d)$$
$$\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}\sum _{d=1}^{\min (n,m)}\mu (d)[d|\gcd (x,y)]$$

把莫比乌斯函数提到前面：

$$\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}[d|\gcd (x,y)]$$

把枚举$i,j$的因子，改成枚举$x,y$：

$$\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{x=1}^n\sum _{y=1}^m[d|\gcd (x,y)]\lfloor \frac{n}{x}\rfloor \lfloor \frac{n}{y}\rfloor$$

把$d$提出来，这样$\gcd (x,y)$就可以消去了：

$$\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{x=1}^{\frac{n}{d}}\sum _{y=1}^{\frac{m}{d}}\lfloor \frac{n}{dx}\rfloor \lfloor \frac{m}{dy}\rfloor$$

然后整理一下：

$$\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{x=1}^{\lfloor \frac{n}{d}\rfloor }\lfloor \frac{\frac{n}{d}}{x}\rfloor \sum _{y=1}^{\frac{m}{d}}\lfloor \frac{\lfloor \frac{m}{d}\rfloor }{y}\rfloor$$

令：

$$f(x)=\sum _{i=1}^x\lfloor \frac{x}{i}\rfloor$$

很容易发现：

$$\begin{gathered} f(\lfloor \frac{n}{d}\rfloor )=\sum _{x=1}^{\lfloor \frac{n}{d}\rfloor }\lfloor \frac{\lfloor \frac{n}{d}\rfloor }{x}\rfloor \\ f(\lfloor \frac{m}{d}\rfloor )=\sum _{x=1}^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{\lfloor \frac{m}{d}\rfloor }{x}\rfloor \end{gathered}$$

所以：

$$Ans=\sum _{d=1}^{\min (n,m)}\mu (d)\ast f(\lfloor \frac{n}{d}\rfloor )\ast f(\lfloor \frac{m}{d}\rfloor )$$

我们只需要处理出$\mu (d)$的前缀和与$f(x)$就行了。
 
$\mathtt{C}\mathtt{o}\mathtt{d}\mathtt{e}:$

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

using namespace std;

int T, N, M;
int p[500100], Mu[500100], f[500100];
bool vis[500100];

void Mobius() {
	Mu[1] = vis[1] = 1;
	for (int i = 2; i <= 5e4; i++) {  //处理μ(d)的前缀和
		if (!vis[i]) Mu[i] = -1, p[++p[0]] = i;
		for (int j = 1; j <= p[0] && i * p[j] <= 5e4; j++) {
			vis[i * p[j]] = true;
			if (i % p[j] == 0) break;
			Mu[i * p[j]] = -Mu[i];
		}
		Mu[i] += Mu[i - 1];
	}
	for (int i = 1; i <= 5e4; i++)  //处理f(x)
		for (int l = 1, r; l <= i; l = r + 1)  {
			r = i / (i / l);
			f[i] += (r - l + 1) * (i / l);
		}
}

void Solve() {
	int ans = 0;
	scanf("%lld%lld", &N, &M);
	if (N > M) swap(N, M);
	for (int l = 1, r; l <= N; l = r + 1) {
		r = min(N / (N / l), M / (M / l));
		ans += (Mu[r] - Mu[l - 1]) * f[N / l] * f[M / l];
	}
	printf("%lld\n", ans);
}

signed main() {
	scanf("%lld", &T);
	Mobius();
	while (T--) Solve();
	return 0;
}