莫比乌斯反演--HDU模板题

题意：http://acm.hdu.edu.cn/showproblem.php?pid=1695

直接上莫比乌斯模板。

 #include <bits/stdc++.h>
 using namespace std;
 const long long maxn=;

 long long mu[maxn],pri[maxn],v[maxn],sum[maxn],cnt=;

 void init()
 {
           mu[]=;
           for(long long i=;i<maxn;i++)
           {
                     if(!v[i])
                     {
                               pri[++cnt]=i;
                               mu[i]=-;
                     }
                     for(long long j=;i*pri[j]<maxn;j++)
                     {
                               v[i*pri[j]]=;
                               if(i%pri[j]==)
                               {
                                         mu[pri[j]*i]=;
                                         break;
                               }
                               mu[i*pri[j]]=-mu[i];
                     }
           }
           sum[]=;
           for(long long i=;i<maxn;i++)
                     sum[i]=sum[i-]+mu[i];
 }
 int main()
 {
           init();
           long long a,b,c,d,k;
           long long t;
           scanf("%lld",&t);
           for(long long ca=;ca<=t;ca++)
           {
                     scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);
                     if(k==)
                     {
                               printf("Case %lld: %lld\n",ca,);
                               continue;
                     }
                     b/=k,d/=k;
                     if(b>d)
                               swap(b,d);
                     long long ans1=,ans2=;
                     for(long long i=,r;i<=b;i=r+)
                     {
                               r=min(b/(b/i),d/(d/i));
                               ans1+=(sum[r]-sum[i-])*(b/i)*(d/i);
                               ans2+=(sum[r]-sum[i-])*(b/i)*(b/i);
                     }
                     printf("Case %lld: %lld\n",ca,ans1-ans2/);
           }
           return ;
 }