luogu P2257 YY的GCD

背景：

以下图片均来自我的$PDF$文件，谢绝转载。

题目传送门：

https://www.luogu.org/problemnew/show/P2257

题意：

思路：

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define _ (int)1e7+10
using namespace std;
	int mu[_],prime[_],g[_];
	bool bz[_];
	int t=0;
	LL sum[_];
void init(int ma)
{
	mu[1]=1;
	bz[0]=bz[1]=true;
	for(int i=2;i<=ma;i++)
	{
		if(!bz[i]) prime[++t]=i,mu[i]=-1,g[i]=1;
		for(int j=1;j<=t&&(LL)i*prime[j]<=ma;j++)
		{
			bz[i*prime[j]]=true;
			if(!(i%prime[j]))
			{
				mu[i*prime[j]]=0;
				g[i*prime[j]]=mu[i];
				break;
			}
			mu[i*prime[j]]=-mu[i];
			g[i*prime[j]]=-g[i]+mu[i];
		}
		sum[i]=sum[i-1]+g[i];
	}
}
int main()
{
	int T,n,m;
	init(_);
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d %d",&n,&m);
		int l=1,r;
		LL ans=0;
		while(l<=min(n,m))
		{
			r=min(n/(n/l),m/(m/l));
			ans+=(LL)(n/l)*(LL)(m/l)*(LL)(sum[r]-sum[l-1]);
			l=r+1;
		}
		printf("%lld\n",ans);
	}
}