hdu 4135 Co-prime 容斥原理

简单的容斥原理！！！

代码如下：

#include<iostream>
#include<stdio.h>
#include<cstring>
#include<algorithm>
#define I(x) scanf("%d",&x)
#define ll __int64
#define MAX 500000
using namespace std;
int prime[MAX],cnt,e[MAX],num;
ll a,b,res1,res2;
bool f[MAX];
void init()
{
    cnt=0;
    memset(f,0,sizeof(f));
    for(int i=2;i<MAX;i++){
        if(!f[i]) prime[cnt++]=i;
        for(int j=0;j<cnt&&i*prime[j]<MAX;j++){
            f[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
void fac(ll n)
{
    num=0;
    for(int i=0;i<cnt&&prime[i]*prime[i]<=n;i++){
        if(n%prime[i]==0){
            e[num++]=prime[i];
            n/=prime[i];
            while(n%prime[i]==0) n/=prime[i];
        }
    }
    if(n>1) e[num++]=n;
}
void dfs(ll sum,int d,int c,ll &ans,ll n)
{
    if(sum>n) return;
    ans-=(ll)(n/sum*c);
    for(int i=d+1;i<num;i++)
        dfs((ll)sum*e[i],i,-c,ans,n);
}
int main(){
    int t,ca=0;
    ll n;
    init();
    I(t);
    while(t--){
        scanf("%I64d%I64d%I64d",&a,&b,&n);
        fac(n);
        res1=a-1;
        res2=b;
        for(int i=0;i<num;i++){
            dfs(e[i],i,1,res1,a-1);
            dfs(e[i],i,1,res2,b);
        }
        printf("Case #%d: %I64d
",++ca,res2-res1);
    }
    return 0;
}