hdu4135 容斥原理

题意：

给出A，B，N，求出[A,B]范围内与N互素的数字的个数。

即B范围内的不与N互素的数减去A范围内不与N互素的数，把[1，A]，[1，B]中不与N互素的数分别求出来，再减掉就是和N互素的数了。

那么首先将N分解质因数，因子和因子的倍数可以被除尽，一定不与N互素，把这些数都记下来，筛法求素数即可~

然后把这些数和倍数加起来？果断会有重复，所以这就到容斥原理了，其实这个画个韦恩图就好~

奇数加偶数减，用队列，位运算，dfs都可以实现，不过那个位运算的搞了好久现在暂时还是没懂。。

下面的代码是队列。

#include<stdio.h>
#include<string.h>
using namespace std;
long long a[1000],num;
void prime(long long n)
{
long long i;
num=0;
for(i=2;i*i<=n;i++)
{
if(n%i==0)
{
a[num++]=i;
while(n%i==0)
n=n/i;
}
}
if(n>1)
a[num++]=n;
}
long long get(long long m)
{
long long que[10000],k,t=0,sum=0;
que[t++]=-1;
for(int i=0;i<num;i++)
{
k=t;
for(int j=0;j<k;j++)
que[t++]=que[j]*a[i]*(-1);
}
for(int i=1;i<t;i++)
sum=sum+m/que[i];
return sum;
}

int main()
{
long long T,x,y,n,cnt;
while(scanf("%lld",&T)!=EOF)
{
for(int i=1;i<=T;i++)
{
scanf("%lld%lld%lld",&x,&y,&n);
prime(n);
cnt=y-get(y)-(x-1-get(x-1));
printf("Case #%d: ",i);
printf("%lld\n",cnt);
}
}
return 0;
}