BZOJ 2301: [HAOI2011]Problem b (莫比乌斯反演+分块+容斥)

因为这里的a和c不一定为1,所以我们要使用容斥
并且因为询问组数较多，我们每次不能线性的使用莫比乌斯反演，我们要对区间进行分块
分块其实很简单,具体看代码吧

我的代码：

#include<cstdio>
#include<algorithm>
using namespace std;

typedef long long ll;
const int maxn=50000+100;

int mu[maxn],sum[maxn];
bool isprime[maxn];
int prime[maxn],tot;

void init(){

    for(int i=0;i<maxn;i++) isprime[i]=true;
    isprime[0]=isprime[1]=false;
    tot=0;
    sum[0]=0;
    sum[1]=mu[1]=1;
    for(int i=2;i<maxn;i++){

        if(isprime[i]) prime[tot++]=i,mu[i]=-1;
        for(int j=0;j<tot && i*prime[j]<maxn;j++){

            isprime[i*prime[j]]=false;
            if(i%prime[j]) mu[i*prime[j]]=-mu[i];
            else break;
        }
    }
    for(int i=2;i<maxn;i++) sum[i]=sum[i-1]+mu[i];
}

ll Mobius(int a,int b){

    if(a>b) swap(a,b);
    int pos;
    ll res=0;
    for(int i=1;i<=a;i=pos+1){

        pos=min(a/(a/i),b/(b/i));
        res+=(ll)(sum[pos]-sum[i-1])*(a/i)*(b/i);
    }
    return res;
}

int main(){

    int n;
    scanf("%d",&n);
    init();
    while(n--){

        int a,b,c,d,k;
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        a--,c--; 
        a/=k;b/=k;
        c/=k;d/=k;
        printf("%lld\n",Mobius(b,d)-Mobius(a,d)-Mobius(b,c)+Mobius(a,c));       
    }

}