这肯定要用卷积的,我们来推一推公式。
我们先令 n=min(n,m) , m=max(m,n)
则
ans=∑p∑ni=1∑mj=1gcd(i,j)==p
ans=∑p∑[np]i=1∑[mp]j=1gcd(i,j)==1
ans=∑p∑[np]i=1∑[mp]j=1∑d|gcd(i,j)μ(d)
ans=∑p∑[np]d=1μ(d)[npd][mpd]
现在,我们就可以在素数个数乘以
n√
下完成一个询问了,但是,这样显然会超时。
继续化简式子,令
k=pd
则
ans=∑nk=1∑p|kμ(kp)[nk][mk]
再令
F(x)=∑p|xμ(xp)
,这个东西可以
nlogn算出来的
最后就可以 n√ 完成一次询问了,这道题就解决了。
#include<bits/stdc++.h>
#define N 10000000
using namespace std;
typedef long long ll;
int prime[N+5],notp[N+5],mu[N+5],cnt,T,n,m,lst;
ll F[N+5],ans;
inline char nc()
{
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
int x=0,b=1;
char c=nc();
for(;!(c<='9'&&c>='0');c=nc())if(c=='-')b=-1;
for(;c<='9'&&c>='0';c=nc())x=x*10+c-'0';
return x*b;
}
inline void init()
{
notp[1]=true,mu[1]=1;
for(int i=2;i<=N;i++)
{
if(!notp[i])prime[++cnt]=i,mu[i]=-1;
for(int j=1;j<=cnt&&i*prime[j]<=N;j++)
{
notp[i*prime[j]]=true;
if(i%prime[j]==0)
{
mu[i*prime[j]]=0;
break;
}
mu[i*prime[j]]=-mu[i];
}
}
for(int i=1;i<=cnt;i++)
{
int p=prime[i];
for(int j=1;j*p<=N;j++)F[j*p]+=mu[j];
}
for(int i=1;i<=N;i++)F[i]+=F[i-1];
}
int main()
{
freopen("in.txt","r",stdin);
init();
T=read();
while(T--)
{
n=read(),m=read();if(n>m)swap(n,m);ans=0;
for(int i=1;i<=n;i=lst+1)
{
lst=min(n/(n/i),m/(m/i));
ans+=(ll)(F[lst]-F[i-1])*(n/i)*(m/i);
}
printf("%lld\n",ans);
}
return 0;
}