BZOJ 1101 [POI2007]Zap 莫比乌斯反演

题目大意：求

看到gcd(i,j)==d 首先想到把i,j同时除d，变成统计互质的数的个数。
令题目转化为求

想统计众多式子中等于1的个数，可以想到莫比乌斯函数的一个性质
（这里的d不同于题目中的d，题目中的d已经没用了）

令上式中的，则代入原式有

将条件转换一下，得到

将式子由枚举i,j转换成枚举d，得到

而(a’/d)与(b’/d)的取值在一定区间内不会发生变化，分块预处理，算法时间复杂度为O(n*sqrt(n))，其中n=min(a’,b’)

#include<cstdio>
#include<algorithm>
#define N 50000
using namespace std;
int mu[N+5],n,ans,prime[N],tot,sum[N];
bool np[N+5];
void init(){
    mu[1]=1;
    for(int i=2;i<=N;i++) {
        if(!np[i]) prime[++tot]=i, mu[i]=-1;
        for(int j=1;prime[j]*i<=N;j++){
            np[i*prime[j]]=true;
            if(!(i%prime[j])) break;
            mu[i*prime[j]]=mu[i]*mu[prime[j]];
        }
    }
    return ;
}
int main() {
    init();
    for(int i=1;i<=N;i++)
        sum[i]=sum[i-1]+mu[i];
    int T;
    scanf("%d",&T);
    while(T--) {
        int a,b,d;
        scanf("%d%d%d",&a,&b,&d);
        if(a>b) swap(a,b);
        a/=d, b/=d, ans=0;
        int l=1,r;
        while(l<=a) {
            int r=min(a/(a/l), b/(b/l));
            ans+=(a/l)*(b/l)*(sum[r]-sum[l-1]);
            l=r+1;
        }
        printf("%d\n",ans);
    }
    return 0;
}