hdu 4135 Co-prime 复习容斥

题意，思路：http://www.cnblogs.com/jiangjing/archive/2013/06/03/3115470.html

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

上面大神链接是用队列写的，我是位运算写的

#include <cstdio>
#include <vector>

using namespace std;
int num[100000],cnt;

long long solve(long long x)
{
    long long ans=0;//ans要设成long long
    for(int i=1; i<(1<<cnt); i++) //这里i一定要从1开始，不然会有错，cnt一般都很小，也就是10
    {
        int cnt1=0,val=1;
        for(int j=0; j<cnt; j++)
        {
            if(i&(1<<j))
            {
                cnt1++;
                val*=num[j];
            }
        }
        if(cnt1&1) ans+=x/val;// a/b表示1到a里面有几个数是可以被b整除的
        else ans-=x/val;
    }
    return x-ans;
}

int main()
{
    long long a,b,n;
    int t;
    scanf("%d",&t);
    for(int kase=1; kase<=t; kase++)
    {
        scanf("%I64d%I64d%I64d",&a,&b,&n);
        cnt=0;
        for(int i=2; i*i<=n; i++) //i即使是int也是可以的
            if(n%i==0)
            {
                num[cnt++]=i;//存的是n的所有的素因子
                while(n%i==0) n/=i;//这里直接用n除也可以
            }
        if(n>1) num[cnt++]=n;//这句一定不要漏了
        long long ans1=solve(a-1);
        long long ans2=solve(b);
        printf("Case #%d: %I64d\n",kase,ans2-ans1);
    }
    return 0;
}