标签：

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

Co-prime

2
1 10 2
3 15 5

Case #1: 5
Case #2: 10

Hint
In the first test case, the five integers in range [1,10] which are relatively prime to 2 are {1,3,5,7,9}. 
 

#include <cstdio>
#include <vector>
using namespace std;
vector<int>fac;
void get_fac(__int64 n){
    for(__int64 i=2;i*i<=n;i++){
        if(n%i==0){
            fac.push_back((int)i);
            while(n%i==0)n/=i;
        }
    }
    if(n!=1)fac.push_back((int)n);
}
__int64 solve(__int64 a){
    //cnt表示有几个集合，lcm是集合的元素个数
    __int64 sum=0;
    int cnt=fac.size();
    for(int num=1;num<(1<<cnt);num++){//cnt个集合，容斥定理中就有2的cnt次方-1项,num是二进制压缩//3为11，代表A1交A2
        __int64 mult=1,ones=0;
        for(int i=0;i<cnt;i++){//i表示第几个集合，例如i=2，表示为100，表示为第三个集合
            if(num&(1<<i)){//当i为0的时候，判断第一个集合是否进行交集操作
            ones++;//当ones为偶数的时候，那么就是减，当noes为奇数的时候为加
            mult*=fac[i];
            }
        }
    if(ones%2) sum+=a/mult;
    else sum-=a/mult;
    }
    return a-sum;
}
int main()
{
    //freopen("F://in.txt","r",stdin);
    int t;
    scanf("%d",&t);
    for(int Case=1;Case<=t;Case++){
        __int64 a,b,n;
        fac.clear();
        scanf("%I64d%I64d%I64d",&a,&b,&n);
        get_fac(n);
        __int64 ans=solve(b)-solve(a-1);
        printf("Case #%d: %I64d\n",Case,ans);
    }
    return 0;
}

Co-prime

标签：

原文地址：http://blog.csdn.net/burning_newbie/article/details/48001177