bzoj 2301

　　

先预处理莫比乌斯函数，再分块求

不会用公式编辑TAT，直接贴题解吧。。

从结论来看貌似也能用容斥原理？

#include<bits/stdc++.h>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define inf 1e9
#define ll long long
#define succ(x) (1<<x)
#define NM 50000+5
using namespace std;
int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int mu[NM],sum[NM],a,b,c,d,k,T,p[NM],tot;
bool check[NM];
int slove(int x,int y){
    int t,ans=0;
    if((x/=k)>(y/=k))swap(x,y);
    for(int i=1;i<=x;i=t+1){
        t=min(x/(x/i),y/(y/i));
        if(t>x)t=x;
        ans+=(sum[t]-sum[i-1])*(x/i)*(y/i);
    }
    return ans;
}
int main(){
    freopen("data.in","r",stdin);
    T=read();
    sum[1]=mu[1]=1;
    inc(i,2,NM-5){
        if(!check[i]){
            p[++tot]=i;
            mu[i]=-1;
        }
        sum[i]=sum[i-1]+mu[i];
        inc(j,1,tot){
            if(i*p[j]>NM-5)break;
            check[i*p[j]]=true;
            if(i%p[j])mu[i*p[j]]=-mu[i];
            else break;
        }
    }
//    inc(i,1,10)printf("%d ",mu[i]);printf("
");
    while(T--){
        a=read();b=read();c=read();d=read();k=read();
        printf("%d
",slove(b,d)-slove(a-1,d)-slove(b,c-1)+slove(a-1,c-1));
    }
    return 0;
}