[bzoj2820]YY的GCD

#include<bits/stdc++.h>
using namespace std;
#define N 10000001
int t,n,m,k,mu[N],vis[N],p[N];
long long ans,f[N];
void linear(){
    mu[1]=1;
    for(int i=2;i<N;i++){
        if (!vis[i]){
            p[++p[0]]=i;
            mu[i]=-1;
            f[i]=1;
        }
        for(int j=1;i*p[j]<N;j++){
            vis[i*p[j]]=1;
            if (i%p[j]){
                mu[i*p[j]]=-mu[i];
                f[i*p[j]]=mu[i]-f[i];
            }
            else{
                mu[i*p[j]]=0;
                f[i*p[j]]=mu[i];
                break;
            }
        }
    }
    for(int i=1;i<N;i++)f[i]+=f[i-1];
}
int main(){
    scanf("%d",&t);
    linear();
    while (t--){
        scanf("%d%d",&n,&m);
        if (n>m)swap(n,m);
        ans=0;
        for(int i=1,j;i<=n;i=j+1){
            j=min(n/(n/i),m/(m/i));
            ans+=(f[j]-f[i-1])*(n/i)*(m/i);
        }
        printf("%lld\n",ans);
    }
}