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=0;
    for(int i=2; i*i<=n; i++) //对n进行素因子分解
    {//筛法求素数
        if(n%i==0)
        {
            prim[k++]=i;
            while(n%i==0)
                n/=i;
        }
    }
    if(n>1)
        prim[k++]=n;
}

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

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

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

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