【BZOJ2820】ygy的gcd

不知道为什么不想写总结，只是（因为是用别人的权限号交的所以）屯一个代码

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 using namespace std;
 #define ll long long
 const ll inf=(ll)2e9;
 ll rd(){ll z=,mk=;  char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')mk=-;  ch=getchar();}
     while(ch>=''&&ch<=''){z=(z<<)+(z<<)+ch-'';  ch=getchar();}
     return z*mk;
 }
 void wt(ll x){if(x<)  putchar('-'),x=-x;
     int wtp=;  char wtc[];
     while(x)  wtc[++wtp]=(x%)+'',x/=;
     while(wtp)  putchar(wtc[wtp--]);
 }
 int n,m;
 int miu[],prm[],prp=;  bool prg[];
 ll f[],s[];
 void gtmiu(){
     memset(prg,,sizeof(prg));
     miu[]=;
     for(int i=;i<=;++i){
         if(!prg[i])  prm[++prp]=i,miu[i]=-;
         for(int j=;j<=prp && prm[j]*i<=;++j){
             prg[i*prm[j]]=true;
             if(!(i%prm[j])){  miu[i*prm[j]]=;  break;}
             miu[i*prm[j]]=-miu[i];
         }
     }
     for(int i=;i<=prp;++i)for(int j=;j*prm[i]<=;++j)
         f[j*prm[i]]+=miu[j];
     for(int i=;i<=;++i)  s[i]=s[i-]+f[i];
 }
 int main(){
     //freopen("ddd.in","r",stdin);
     gtmiu();
 int T;  cin>>T;  while(T--){
     n=rd(),m=rd();
     ll ans=;
     if(n>m)  swap(n,m);
     for(int i=,j;i<=n;i=j+){
         j=min(n/(n/i),m/(m/i));
         ans+=(s[j]-s[i-])*(n/i)*(m/i);
     }
     wt(ans),putchar('\n');
     continue;
 }
     return ;
 }