不妨设
n≥m
。
==∑i=1n∑j=1mf(gcd(i,j))∑d=1nf(d)∑x=1⌊nd⌋μ(x)⌊ndx⌋⌊mdx⌋∑D=1n∑d|Df(d)μ(Dd)⌊ndx⌋⌊mdx⌋
记 g=f∗μ ,只要我们线性筛出 g ,就可以
设 g(n)=∑d|nf(d)μ(nd) , n=∏ki=1pqii , d=∏ki=1pixi 。显然 xi=qi 或者 xi=qi−1 ,否则 μ(nd) 一定为 0 。如果
#include<cstdio>
#include<algorithm>
using namespace std;
#define LL long long
const int maxn=10000010;
int mnp[maxn],mnq[maxn],g[maxn],sg[maxn],prm[maxn],tot;
void init()
{
for (int i=2;i<=maxn;i++)
{
if (!mnp[i])
{
mnp[i]=i;
mnq[i]=1;
g[i]=1;
prm[++tot]=i;
}
for (int j=1;j<=tot&&(LL)i*prm[j]<=maxn;j++)
if (i%prm[j])
{
mnp[i*prm[j]]=prm[j];
mnq[i*prm[j]]=1;
if (mnq[i]==1) g[i*prm[j]]=-g[i];
}
else
{
mnp[i*prm[j]]=mnp[i]*prm[j];
mnq[i*prm[j]]=mnq[i]+1;
if (i/mnp[i]==1) g[i*prm[j]]=1;
else if (mnq[i/mnp[i]]==mnq[i*prm[j]]&&mnp[i/mnp[i]]!=mnp[i*prm[j]]) g[i*prm[j]]=-g[i/mnp[i]];
break;
}
sg[i]=sg[i-1]+g[i];
}
}
LL solve(int n,int m)
{
LL ret=0;
for (int i=1,j;i<=n&&i<=m;i=j+1)
{
j=min(n/(n/i),m/(m/i));
ret+=(LL)(sg[j]-sg[i-1])*(n/i)*(m/i);
}
return ret;
}
int main()
{
init();
int T,n,m;
scanf("%d",&T);
while (T--)
{
scanf("%d%d",&n,&m);
printf("%lld\n",solve(n,m));
}
}