bzoj 1101 [POI2007]Zap——反演

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1101

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=5e4+;
int T,n,m,D,u[N],pri[N];bool vis[N];
int rdn()
{
  int ret=;bool fx=;char ch=getchar();
  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return fx?ret:-ret;
}
void init()
{
  int lm=5e4,cnt=;u[]=;
  for(int i=;i<=lm;i++)
    {
      if(!vis[i])pri[++cnt]=i,u[i]=-;
      for(int j=;j<=cnt&&(ll)i*pri[j]<=lm;j++)
    {
      vis[i*pri[j]]=;
      if(i%pri[j]==){u[i*pri[j]]=;break;}
      u[i*pri[j]]=-u[i];
    }
      u[i]+=u[i-];
    }
}
int main()
{
  T=rdn();init();
  while(T--)
    {
      n=rdn();m=rdn();D=rdn();
      n/=D;m/=D;if(n>m)swap(n,m);
      ll ans=;
      for(int i=,d0,d1,j;i<=n;i=j+)
    {
      d0=n/i;d1=m/i; j=min(n/d0,m/d1);
      ans+=(ll)(u[j]-u[i-])*d0*d1;
    }
      printf("%lld\n",ans);
    }
  return ;
}