bzoj 2301: [HAOI2011]Problem b

#include<cstdio>
#include<iostream>
#include<cstring>
#define M 50009
#define ll long long
using namespace std;
int n,a,b,c,d,k,mo[M+5],zhan[M],cnt,sum[M];
bool mark[M+5];
void mobiwus()
{
    mo[1]=1;
    for(int i=2;i<=M;i++)
      {
        if(!mark[i])
          {
            mo[i]=-1;
            zhan[++cnt]=i;
          }
        for(int j=1;zhan[j]*i<=M&&j<=cnt;j++)
          {
            mark[zhan[j]*i]=1;
            if(i%zhan[j])
              mo[i*zhan[j]]=mo[i]*mo[zhan[j]];
            else
              mo[i*zhan[j]]=0;
          }
      }
}
ll su(int n,int m)
{
    n/=k;
    m/=k;
    if(n>m)
      swap(n,m);
    int ans=0,last;
    for(int i=1;i<=n;i=last+1)
      {
        last=min(n/(n/i),m/(m/i));
        ans+=(ll)(n/i)*(m/i)*(sum[last]-sum[i-1]);
      }
      return ans;
}
int main()
{
    mobiwus();
    for(int i=1;i<=M;i++)
      sum[i]=sum[i-1]+mo[i];
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
      {
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        printf("%lld
",su(b,d)-su(a-1,d)-su(b,c-1)+su(a-1,c-1));
      }
    return 0;
}

莫比乌斯反演