莫比乌斯反演——[HAOI2011]Problem b

题目链接

题解：
列出式子是一个很明显的莫比乌斯反演，然后再加上容斥就可以很快过了。

#include <bits/stdc++.h>
#define int long long
using namespace std;
typedef long long ll;
const int N=5e5+5;
const int MOD=1e9+7;
int a,b,c,d,k;
int prime[N],vis[N],mob[N],cnt;
void init()
{
    mob[1]=1;
    for(int i=2;i<N;i++) {
        if(!vis[i]) prime[cnt++]=i,mob[i]=-1;
        for(int j=0;prime[j]*i<N&&j<cnt;j++){
            vis[prime[j]*i]=1;
            if(i%prime[j]==0) {
                mob[i*prime[j]]=0;
                break;
            }mob[i*prime[j]]=-mob[i];
        }
    }
    for(int i=1;i<N;i++) mob[i]+=mob[i-1];
}
ll work(int n,int m)
{
    ll res=0;
    n/=k,m/=k;
    for(int l=1,r,x;l<=min(n,m);l=r+1){
        r=min(n/(n/l),m/(m/l));
        res+=(mob[r]-mob[l-1])*(n/l)*(m/l);
    }
    return res;
}
signed main()
{
    init();
    int T; scanf("%d",&T);
    while(T--) {
        scanf("%lld%lld%lld%lld%lld", &a, &b, &c, &d, &k);
        printf("%lld\n",work(b,d)-work(a-1,d)-work(c-1,b)+work(a-1,c-1));
    }
}