baoj 2301: [HAOI2011]Problem b （莫比乌斯反演）

题目大意：

对于给出的n个询问，每次求有多少个数对(x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y) = k，gcd(x,y)函数为x和y的最大公约数。

思路：

先看懂题目设ans[n][m]为1 < i < n,1 < j < m，gcd(i,j)=k的个数，那么答案就是ans[b][d]-ans[a-1][d]-ans[b][c-1]+ans[a-1][c-1]。
我们先设f[i]为gcd(x,y)=i的个数，这东西好像不太好求，所以我们引入莫比乌斯反演。多设一个F[i]为i|gcd（x，y）。这个F显然很好求F[i]=(x/i)*(y/i)除法为下取整，然后就是反演！！一般形式为
第一个$F[i]=\sum_{d|i}f[d]=>f[n]=\sum_{d|n}u[d]*F[n/d]$
第二个$F[i]=\sum_{i|d}f[d]=>f[n]=\sum_{n|d}u[d/n]*F[d]$
下取整可以根号做，很ok。

程序：

#include<cstdio>
#include<algorithm>
#include<iostream>
#define N 500005
#define LL long long

LL ans,t,a,b,c,d,k,cnt;
LL sum[N],mu[N],p[N],np[N];

void getmu(LL n){
    mu[1]=1;
    for (LL i=2;i<=n;i++){
        if (!np[i]){ 
            mu[i]=-1;
            p[++cnt]=i;
        }
        for (LL j=1;j<=cnt&&i*p[j]<=n;j++){
            np[i*p[j]]=1;
            if (i%p[j]==0){
                mu[i*p[j]]=0;
                break;
            }
            mu[i*p[j]]=-mu[i];
        }
    }
    sum[1]=1;
    for (LL i=2;i<=n;i++) sum[i]=sum[i-1]+mu[i];
}

LL cal(LL n,LL m,LL k){
    if (n>m) std::swap(n,m);
    LL ans=0,j=0;
    for (LL i=1;i<=n/k;i=j+1){
        j=std::min(n/(n/i/k)/k,m/(m/i/k)/k);
        ans+=(n/i/k)*(m/i/k)*(sum[j]-sum[i-1]);
    }
    return ans;
}

int main(){
    getmu(500000);
    scanf("%lld",&t);
    while (t){
        scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);
        ans=cal(b,d,k)-cal(a-1,d,k)-cal(c-1,b,k)+cal(a-1,c-1,k);
        printf("%lld\n",ans);
        t--;
    }
}