P2257 YY的GCD （莫比乌斯反演）

题目链接：点击这里

题目大意：
$t$ 组数据，每组数据给出 $N,M$ ，求 $1\le x\le N$ ，$1\le y\le M$ 且 $\gcd (x,y)$ 为质数的二元组 $(x,y)$ 的数目

题目分析：
题目所求式子为：
$$\sum _{i=1}^N\sum _{j=1}^M[\gcd (i,j)=prime]$$
我们套路的设$F(n)$为$\gcd (i,j)\%n=0$ 的个数 ， $f(d)$ 为 $\gcd (i,j)=d$ 的个数，即：
$$F(n)=\sum _{n|d}f(d)=\lfloor \frac{N}{n}\rfloor \lfloor \frac{M}{n}\rfloor$$
$$f(d)=\sum _{i=1}^N\sum _{j=1}^M[\gcd (i,j)=d]$$
由莫比乌斯反演可得：
$$f(n)=\sum _{n|d}\mu (\lfloor \frac{d}{n}\rfloor )F(d)$$
于是题目所求可以变成：
$$\sum _{p\in prime}f(p)=\sum _{p\in prime}\sum _{p|d}\mu (\lfloor \frac{d}{p}\rfloor )F(d)$$
$$=\sum _{p\in prime}\sum _{d=1}^{min(\lfloor \frac{N}{n}\rfloor ,\lfloor \frac{M}{n}\rfloor )}\mu (d)F(dp)$$
$$=\sum _{p\in prime}\sum _{d=1}^{min(\lfloor \frac{N}{n}\rfloor ,\lfloor \frac{M}{n}\rfloor )}\mu (d)\lfloor \frac{N}{dp}\rfloor \lfloor \frac{M}{dp}\rfloor$$
我们令 $x=dp$ ，此时有：
$$\sum _{p\in prime}\sum _{x=1}^{min(N,M)}\mu (\frac{x}{p})\lfloor \frac{N}{x}\rfloor \lfloor \frac{M}{x}\rfloor$$
$$=\sum _{x=1}^{min(N,M)}\sum _{p|x,p\in prime}\mu (\frac{x}{p})\lfloor \frac{N}{x}\rfloor \lfloor \frac{M}{x}\rfloor$$
$$=\sum _{x=1}^{min(N,M)}\lfloor \frac{N}{x}\rfloor \lfloor \frac{M}{x}\rfloor \sum _{p|x,p\in prime}\mu (\frac{x}{p})$$

推到这里就可以放心大胆的搞了，用整除分块来快速处理 $\lfloor \frac{N}{n}\rfloor$ 和 $\lfloor \frac{M}{n}\rfloor$ ，然后我们发现 $\mu (\frac{x}{p})$ 是可以提前处理的，我们提前打好这东西的前缀和然后套在整除分块中 $O(1)$ 求出 ${\sum }_{p|x,p\in prime}\mu (\frac{x}{p})$ 即可，时间复杂度为 $O(t\sqrt{n})$

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 1e7+5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
int cnt,pri[maxn],mu[maxn],g[maxn];
ll sum[maxn];
bool vis[maxn];
void get_mu() //筛莫比乌斯函数，求需要的函数的前缀和
{
	mu[1] = vis[1] = 1;
	for(int i = 2;i < maxn;i++)
	{
		if(!vis[i])
		{
			mu[i] = -1;
			pri[++cnt] = i;
		}
		for(int j = 1;j <= cnt && i*pri[j] < maxn;j++)
		{
			vis[i*pri[j]] = true;
			if(i%pri[j] != 0) mu[i*pri[j]] -= mu[i];
			else break;
		}
	}
	for(int i = 1;i <= cnt;i++)
		for(int j = 1;j*pri[i] < maxn;j++)
			g[j*pri[i]] += mu[j];
	for(int i = 1;i < maxn;i++)
		sum[i] = sum[i-1]+g[i];
}
int main()
{
	int t = read();
	get_mu();
	while(t--)
	{
		int n = read(),m = read();
		ll ans = 0;
		for(int l = 1,r;l <= min(n,m);l = r+1) //整除分块
		{
			r = min(n/(n/l),m/(m/l));
			ans += 1ll*(n/l)*(m/l)*(sum[r]-sum[l-1]);
		}
		printf("%lld\n",ans);
	}
	return 0;
}