HDU 4135 Co-prime

思路：直接用求(b,1)范围内互质的数，(a-1,1)范围内互质的数。再求反

就是敲一下容斥模板

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#include<stack>
#include<algorithm>
using namespace std;
#define clc(a,b) memset(a,b,sizeof(a))
#define inf 0x3f3f3f3f
const int N=10010;
#define LL long long
// inline int r(){
//     int x=0,f=1;char ch=getchar();
//     while(ch>'9'||ch<'0'){if(ch=='-') f=-1;ch=getchar();}
//     while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
//     return x*f;
// }

int T;
LL a,b,n;
int pri[110];
int k;
void fun(LL x){
    k=0;
    for(int i=2;i*i<=x;i++){
        if(x&&x%i==0){
            pri[k++]=i;
            while(x&&x%i==0){
                x/=i;
            }
        }
    }
    if(x>1) pri[k++]=x;
}
LL exclu(LL num,int m){
    LL ans=0,tmp,times;
    for(int i=1;i<(LL)(1<<m);i++){
        tmp=1,times=0;
        for(int j=0;j<m;j++){
            if(i&((LL)(1<<j)))
                times++,tmp*=pri[j];
        }
        if(times&1) ans+=num/tmp;
        else ans-=num/tmp;
    }
    return ans;
}
int main(){
    scanf("%d",&T); int cas=1;
    while(T--){
        scanf("%I64d%I64d%I64d",&a,&b,&n);
        clc(pri,0);
        fun(n);
        LL ans1=exclu(a-1,k);
        LL ans2=exclu(b,k);
        printf("Case #%d: %I64d
",cas++,b-ans2-(a-1-ans1));
    }
    return 0;
}