【bzoj 2820】YY的GCD（莫比乌斯反演）

传送门biu~
参考博客：莫比乌斯反演

#include<bits/stdc++.h>
#define N 10000000
using namespace std;
int T,n,m,prime[N+5],mu[N+5],f[N+5];
bool mark[N+5];
void GetPrime(){
    mu[1]=1;
    for(int i=2;i<=N;i++){
        if(!mark[i]) prime[++prime[0]]=i,mu[i]=-1,f[i]=1;
        for(int j=1;j<=prime[0] && i*prime[j]<=N;j++){
            mark[i*prime[j]]=1;
            if(i%prime[j]==0){
                mu[i*prime[j]]=0;
                f[i*prime[j]]=mu[i];
                break;
            }
            mu[i*prime[j]]=-mu[i];
            f[i*prime[j]]=mu[i]-f[i];
        }
    }
    for(int i=1;i<=N;i++) f[i]+=f[i-1];
}
long long calc(int n,int m){
    if(n>m) swap(n,m);
    long long res=0;
    for(int r,i=1;i<=n;i=r+1){
        r=min(n/(n/i),m/(m/i));
        res+=(long long)(f[r]-f[i-1])*(n/i)*(m/i);
    }
    return res;
}
int main() {
    GetPrime();
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        printf("%lld\n",calc(n,m));
    }
    return 0;
}