hdu4135（容斥原理）

题目连接：HDU - 4135

开始学容斥原理。

#include<cstdio>
#include<cstring>
#define LL long long
const int maxn=100;

LL prime[maxn];
LL a,b,n;
LL func(LL up,int m)
{
    LL ans=0,temp=1;
    int flag;
    for(LL i=1;i<(1ll<<m);i++)
    {
        temp=1,flag=0;
        for(int j=0;j<m;j++) if(i&(1ll<<j))
            ++flag,temp*=prime[j];
        if(flag&1)  //容斥原理，奇加偶减
            ans+=up/temp;
        else ans-=up/temp;
    }
    return ans;
}

int factor(LL n)  //对n进行素数分解
{
    int num=0;
    for(LL i=2;i*i<=n;i++) if(n&&n%i==0)
    {
        prime[num++]=i;
        while(n&&n%i==0) n/=i;
    }
    if(n>1)
        prime[num++]=n;
    return num;
}
int main()
{
    int t;
    scanf("%d",&t);
    for(int i=1;i<=t;i++)
    {
        int num=0;
        scanf("%lld%lld%lld",&a,&b,&n);
        num=factor(n);
        LL ans=(b-func(b,num))-(a-1-func(a-1,num));
        printf("Case #%d: %lld
",i,ans);
    }
}