HDOJ4135 Co-prime （容斥原理）

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

求A到B区间与N互质的数的个数，因为A，B不一定等于N，所以无法使用欧拉函数。正解是先求出区间内与N不互质的数的个数m，再用总数减去m即为答案。
N可分解为若干质因子，质因子的倍数与N不互质。所以我们先要求出N的质因子。但质因子的倍数必有重复，比如2和3的公共倍数有6，所以我们需要用容斥定理排重。
我们要用容斥定理分别求出【1，A】和【1，B】区间内与N不互质的数，相减得出【A,B】与N不互质的数，最后再 求出与N互质的数。

代码如下：

#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;

const int maxn = 1e3+10;
const int maxm = 1e4+10;
__int64 a[maxn],num;

//求质因子
void prime_factor(__int64 n){
    __int64 i;
    num=0;
    for (i=2;i*i<=n;i++){
        if (n%i==0) {
            a[num++]=i;
            while (n%i==0) n/=i;
        }
    }
    if (n>1) a[num++]=n;
}

__int64 iep(__int64 n){//使用队列数组实现
    __int64 que[maxm],i,j,k,t=0,sum=0;
    que[t++]=-1;
    for (i=0;i<num;i++){
        k=t;
        for (j=0;j<k;j++) que[t++]=que[j]*a[i]*(-1);
    }
    for (i=1;i<t;i++)  sum+=n/que[i];
    return sum;
}

int main(){
    int t,i;
    __int64 x,y,n,total;
    cin>>t;
    for (i=1;i<=t;i++){
        scanf("%I64d%I64d%I64d",&x,&y,&n);
        prime_factor(n);
        total=y-iep(y)-(x-1-iep(x-1));
        printf("Case #%d: ",i);
        printf("%I64d\n",total);
    }
    return 0;
}