hdu4135Coprime（容斥原理）

http://acm.hdu.edu.cn/showproblem.php?pid=4135

交了十几次，。最后发现理解错题意了，区间包含右端点，一直按不包含做的。。

求出N的质因数 根据容斥原理解出区间数能被质因数整除的 再减掉

#include <iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
#define LL __int64
LL sum1,sum2,a,b,n;
int g,x[1001];
void dfs(int i,LL y,int num)
{
    int j;
    if(y!=x[i])
    y = y*x[i];
    if(num%2!=0)
    {
        sum1+=(a-1)/y;
        sum2+=b/y;
    }
    else
    {
        sum2-=b/y;
        sum1-=(a-1)/y;
    }
    for(j = i+1 ; j <= g ; j++)
    dfs(j,y,num+1);
}
int main()
{
    int t,i,j,k = 0;
    cin>>t;
    while(t--)
    {
        k++;
        scanf("%I64d%I64d%I64d",&a,&b,&n);
        sum1 = 0;
        sum2 = 0;
        g = 0;
        x[g] = 0;
        LL yy = n;
        for(i = 2 ; i <= sqrt(yy*1.0) ; i++)
        {
           if(n%i==0)
           {
               g++;
               x[g] = i;
               while(n%i==0)
               n = n/i;
           }
        }
        if(n!=1)
        {
            g++;
            x[g] =n;
        }
        for(i = 1; i <= g ; i++)
        dfs(i,x[i],1);
        printf("Case #%d: ",k);
        LL sum = b-sum2-(a-1-sum1);
        printf("%I64d\n",sum);
    }
    return 0;
}