HDU 4135

状态压缩

其实刚开始没看懂为什么要用位运算

后来看了别人的注释

逐渐明白

我也加上注释吧

#include <iostream>
using namespace std;

long long Prime[50];//存放N的质因子


long long getNonCoPrime(long long num,int m)//计算非互质数的多少
{
    long long ans=0,flag,temp;
    for(long long i =1;i<(long long)(1<<m);i++)//用二进制串来表示每个质因子是否被取
    {
        temp = 1;
        flag = 0;
        for(int j = 0;j<m;j++)//判断每位是否取到
        {
            if(i&((long long)(1<<j)))//如果取了
            {
                flag++;
                temp*=Prime[j];
            }
        }
        if(flag&1)//容斥原理 ，奇加偶减 
        {
                ans+=num/temp;
        }else{
            ans-=num/temp;
        }
        
    }
    return ans;
}


int main()
{
    int T;
    cin>>T;
    int t =T;
    
    while(T--)
    {
        long long A,B;
        int N;
        cin>>A>>B>>N;
        int j =0;
        for(int i = 2;i*i<=N;i++)//获得质因子
        {
            if(N&&N%i==0)
            {
                Prime[j]=i;
                j++;
                while(N&&N%i==0)//避免重复
                {
                    N/=i;
                }
            }
        }
        if(N>1)
            Prime[j++]=N;
        long long result = B-A-(getNonCoPrime(B,j) -getNonCoPrime(A-1,j))+1;
        cout<<"Case #"<<t-T<<": "<<result<<endl;
    }
    return 0;
}