BZOJ 2820 YY的GCD 莫比乌斯反演

题意：求有多少个数对(x,y)，使得x<=m,y<=n且GCD(x,y)为质数

解法：具体去见ACdream的博客里面讲的还是很详细的 这里

///BZOJ 2820

#include <bits/stdc++.h>
using namespace std;
const int maxn = 10000000;
typedef long long LL;
bool mark[10000005];
int pri[1000005],mu[10000005];
LL sum[10000005];
int tot;
void GetMobius(){
    mu[1]=1;
    tot=0;
    for(int i=2; i<maxn; i++){
        if(!mark[i]) pri[++tot]=i, mu[i]=-1;
        for(int j=1; j<=tot&&i*pri[j]<maxn; j++){
            mark[i*pri[j]]=1;
            if(i%pri[j]==0){
                mu[i*pri[j]]=0;
                break;
            }
            else mu[i*pri[j]]=-mu[i];
        }
    }
    for(int i=1; i<=tot; i++){
        for(int j=pri[i]; j<maxn; j+=pri[i]){
            sum[j]+=mu[j/pri[i]];
        }
    }
    for(int i=1; i<maxn; i++) sum[i]+=sum[i-1];
}

int main(){
    GetMobius();
    int T;
    int n, m;
    scanf("%d", &T);
    while(T--){
        scanf("%d%d", &n,&m);
        LL ans = 0LL;
        int x=min(n,m), pos;
        for(int d=1; d<=x; d=pos+1){
            pos=min(n/(n/d), m/(m/d));
            ans+=(sum[pos]-sum[d-1])*(n/d)*(m/d);
        }
        printf("%lld\n", ans);
    }
    return 0;
}