【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=0,mk=1;  char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')mk=-1;  ch=getchar();}
    while(ch>='0'&&ch<='9'){z=(z<<3)+(z<<1)+ch-'0';  ch=getchar();}
    return z*mk;
}
void wt(ll x){if(x<0)  putchar('-'),x=-x;
    int wtp=0;  char wtc[52];
    while(x)  wtc[++wtp]=(x%10)+'0',x/=10;
    while(wtp)  putchar(wtc[wtp--]);
}
int n,m;
int miu[11000000],prm[11000000],prp=0;  bool prg[11000000];
ll f[11000000],s[11000000];
void gtmiu(){
    memset(prg,0,sizeof(prg));
    miu[1]=1;
    for(int i=2;i<=10000000;++i){
        if(!prg[i])  prm[++prp]=i,miu[i]=-1;
        for(int j=1;j<=prp && prm[j]*i<=10000000;++j){
            prg[i*prm[j]]=true;
            if(!(i%prm[j])){  miu[i*prm[j]]=0;  break;}
            miu[i*prm[j]]=-miu[i];
        }
    }
    for(int i=1;i<=prp;++i)for(int j=1;j*prm[i]<=10000000;++j)
        f[j*prm[i]]+=miu[j];
    for(int i=1;i<=10000000;++i)  s[i]=s[i-1]+f[i];
}
int main(){
    //freopen("ddd.in","r",stdin);
    gtmiu();
int T;  cin>>T;  while(T--){
    n=rd(),m=rd();
    ll ans=0;
    if(n>m)  swap(n,m);
    for(int i=1,j;i<=n;i=j+1){
        j=min(n/(n/i),m/(m/i));
        ans+=(s[j]-s[i-1])*(n/i)*(m/i);
    }
    wt(ans),putchar('
');
    continue;
}
    return 0;
}