洛谷P3327：[SDOI2015]约数个数和（莫比乌斯反演 | 莫比乌斯函数）

---

首先 $d(i\ast j)={\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]$$
先改变一下枚举项，枚举因子 x,y：$$\sum _{x=1}^N\sum _{y=1}^M\sum _{i=1}^{\lfloor \frac{N}{x}\rfloor }\sum _{j=1}^{\lfloor \frac{M}{y}\rfloor }[gcd(x,y)=1]$$
然后用${\sum }_{d|gcd(i,j)}\mu (d)$替换掉$[gcd(i,j)=1]$：$$\sum _{x=1}^N\sum _{y=1}^M\sum _{i=1}^{\lfloor \frac{N}{x}\rfloor }\sum _{j=1}^{\lfloor \frac{M}{y}\rfloor }\sum _{d|gcd(x,y)}\mu (d)$$
再次改变枚举项，枚举 d，将 d 提取到前面来：$$\sum _{d=1}^N\mu (d)\sum _{x=1}^{\lfloor \frac{N}{d}\rfloor }\sum _{y=1}^{\lfloor \frac{M}{d}\rfloor }\sum _{i=1}^{\lfloor \frac{N}{x\ast d}\rfloor }\sum _{j=1}^{\lfloor \frac{M}{y\ast d}\rfloor }$$
把后面两个求和去掉：$$\sum _{d=1}^N\mu (d)\sum _{x=1}^{\lfloor \frac{N}{d}\rfloor }\lfloor \frac{N}{x\ast d}\rfloor \sum _{y=1}^{\lfloor \frac{M}{d}\rfloor }\lfloor \frac{M}{y\ast d}\rfloor$$
到这一步就可以做了，后面两个可以数论分块，$O(n\ast \sqrt{n})$预处理一下${\sum }_{i=1}^p\lfloor \frac{p}{i}\rfloor$就可以O（1）算出后面一坨，再预处理一下莫比乌斯函数前缀和，就可以$O(\sqrt{n})$处理单组，多组复杂度为：$O(T\ast \sqrt{n})$

---

另外一种推法：我们要求的是$$\sum _{i=1}^N\sum _{j=1}^M\sum _{x|i}\sum _{y|j}[gcd(x,y)=1]$$
改变枚举项，枚举 x,y：$$\sum _{x=1}^N\sum _{y=1}^M\sum _{i=1}^{\lfloor \frac{N}{x}\rfloor }\sum _{j=1}^{\lfloor \frac{M}{y}\rfloor }[gcd(x,y)=1]$$
接下来是不同的地方：设$$f(d)=\sum _{x=1}^N\sum _{y=1}^M\sum _{i=1}^{\lfloor \frac{N}{x}\rfloor }\sum _{j=1}^{\lfloor \frac{M}{y}\rfloor }[gcd(x,y)=d]$$
$$F(p)=\sum _{p|d}f(d)=\sum _{x=1}^N\sum _{y=1}^M\sum _{i=1}^{\lfloor \frac{N}{x}\rfloor }\sum _{j=1}^{\lfloor \frac{M}{y}\rfloor }[p|gcd(x,y)]$$
把 p 提出来：枚举 x,y是p的倍数，i,j 是 $p\ast x,p\ast y$的倍数（因为i,j是 x,y的倍数，x,y现在是p的倍数，所以i,j是 $p\ast x,p\ast y$的倍数）：$$F(p)=\sum _{p|d}f(d)=\sum _{x=1}^{\lfloor \frac{N}{p}\rfloor }\sum _{y=1}^{\lfloor \frac{M}{p}\rfloor }\sum _{i=1}^{\lfloor \frac{N}{x\ast p}\rfloor }\sum _{j=1}^{\lfloor \frac{M}{y\ast p}\rfloor }1=\sum _{x=1}^{\lfloor \frac{N}{p}\rfloor }\sum _{y=1}^{\lfloor \frac{M}{p}\rfloor }\lfloor \frac{N}{x\ast p}\rfloor \ast \lfloor \frac{M}{y\ast p}\rfloor$$

反演一下：$$f(d)=\sum _{d|p}\mu (\frac{p}{d})\ast F(p)$$(这里p 最多枚举到 n)
我们要求的是： $$f(1)=\sum _{p=1}^n\mu (p)\ast F(p)=\sum _{p=1}^n\mu (p)\ast \sum _{x=1}^{\lfloor \frac{N}{p}\rfloor }\sum _{y=1}^{\lfloor \frac{M}{p}\rfloor }\lfloor \frac{N}{x\ast p}\rfloor \ast \lfloor \frac{M}{y\ast p}\rfloor$$
这里已经上面一样，可以分块做了

---

代码：

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e7 + 10;
int t,n,m;
typedef long long ll;
bool ispri[maxn];
int pri[maxn],mu[maxn];
ll sum[maxn];
void sieve(int n) {
	ispri[0] = ispri[1] = true;
	pri[0] = 0;
	mu[1] = 1;
	for(int i = 2; i <= n; i++) {
		if(!ispri[i]) pri[++pri[0]] = i,mu[i] = -1;
		for(int j = 1; j <= pri[0] && i * pri[j] <= n; j++) {
			ispri[i * pri[j]] = true;
			if(i % pri[j] == 0) break;
			mu[i * pri[j]] = -mu[i];
		}
	}
	for(int i = 1; i <= n; i++) 
		mu[i] += mu[i - 1]; 
}
void prework(int n) {
	int i,j;
	for(i = 1; i <= n; i = j + 1) {
		j = n / (n / i);
		sum[n] += 1ll * (j - i + 1) * (n / i);
	}
}
int main() {
	scanf("%d",&t);
	sieve(100000);
	for(int i = 1; i <= 50000; i++)
		prework(i);
	while(t--) {
		scanf("%d%d",&n,&m);
		if(n > m) swap(n,m);
		int i,j;
		ll res = 0;
		for(i = 1; i <= n; i = j + 1) {
			j = min(n / (n / i),m / (m / i));
			int p1 = n / i,p2 = m / i;
			res += 1ll * (mu[j] - mu[i - 1]) * (sum[p1] * sum[p2]);
		}
		printf("%lld\n",res);
	}
	return 0;
}