HDU4135容斥原理

 #include <cstdio>
 #include <string.h>
 #include <cmath>

 using namespace std;

 #define MAXSIZE 40000
 #define LL long long
 int prim[MAXSIZE];
 LL A,B,n,odd,even;
 int k;//k记录素数的个数

 void findPrim(LL n)
 {
     k=;
     for(int i=; i*i<=n; i++) //对n进行素因子分解
     {//筛法求素数
         if(n%i==)
         {
             prim[k++]=i;
             while(n%i==)
                 n/=i;
         }
     }
     if(n>)
         prim[k++]=n;
 }

 void getCount(int N)
 {
     long long mul=;
     int count=;
     int m=;
     while(N){
         if(N&) {
             count++;
             mul*=prim[m];
         }
         m++;
         N>>=;
     }
     long long t=B/mul-(A-)/mul;
     if(count%==) even+=t;
     else odd+=t;
 }

 int main()
 {
     int T,log;
     scanf("%d",&T);
     for(int x=;x<=T;x++)
     {
         odd=,even=,log=;
         scanf("%I64d%I64d%I64d",&A,&B,&n);
         findPrim(n);

         /*for(int i=0;i<k;i++) printf("%d ",prim[i]);
         printf("\n");*/

         while(log<(<<k)){
             getCount(log);
             log++;
         }
         printf("Case #%d: ",x);
         printf("%I64d\n",B-A+-odd+even);
     }
     return ;
 }