莫比乌斯反演初步2：BZOJ 2301: [HAOI2011]Problem b

题目描述：戳这里
题解：
这题看似和HDU1965很像，其实就是很像。（HDU1965题解）
反演的过程是一样的，只是这题的数据更加大一点。所以我们一定要用除法分块。
除法的原理可以看这个blog，讲得很清楚。
还有一点，这题的a和c不是1，所以要容斥一下：
答案是：mu(b,d)-mu(a-1,d)-mu(b,c-1)+mu(a-1,c-1)

代码如下：

#include<cstdio>
#include<string>
#include<cstring>
#define ll long long
using namespace std;
const int maxn=50005;
int T,m,a,b,c,d,k,su[maxn],f[maxn];
bool vis[maxn];
void make_p(){
    f[1]=vis[1]=1;
    for (int i=2;i<=50000;i++){
        if (!vis[i]) su[++m]=i,f[i]=-1;
        for (int j=1;j<=m&&i*su[j]<=50000;j++) {
            vis[i*su[j]]=1;
            if (i%su[j]==0) {f[i*su[j]]=0; break;}
            f[i*su[j]]=-f[i];
        }
    }
    for (int i=2;i<=50000;i++) f[i]+=f[i-1];
}
ll getmu(int x,int y){
    x/=k； y/=k;
    int minn=min(x,y);
    int r=0; ll ret=0;
    for (int i=1;i<=minn;i=r+1) {
        r=min(x/(x/i),y/(y/i));
        ret+=1ll*(f[r]-f[i-1])*(x/i)*(y/i);
    }
    return ret;
}
int main(){
    make_p();
    scanf("%d",&T);
    while (T--) {
        scanf("%d %d %d %d %d",&a,&b,&c,&d,&k);
        printf("%lld\n",getmu(b,d)-getmu(a-1,d)-getmu(b,c-1)+getmu(a-1,c-1));
    }
    return 0;
}