hdoj-4315 (容斥原理）

问题描述
给定数字N，要求您计算A和B之间的整数（与N相对质数）。相对于N而言，如果
两个整数除1以外没有共同的正除数，则称这两个整数为互质或相对质数。同样，如果它们的最大公约数为1，则数字1对每个整数都是质数。
输入值
输入的第一行包含T（0 <T <= 100）个测试用例，接下来的T行中的每一行包含三个整数A，B，N，其中（1 <= A <= B <= 10 15）和（1 <= N <= 10 9）。
输出量
对于每个测试用例，请打印A到B之间（包括N在内）的整数数目。请遵循以下输出格式。
样本输入
2
1 10 2
3 15 5
样本输出
案例1：5
案例2：10
暗示
在第一个测试用例中，相对于2质数为[1,10]的五个整数为{1,3,5,7,9}。

题意：给定一段区间 a-b 一个数n 求在a-b区间内与n互质数的个数
互质性质很容易联想到欧拉定理，但是区间范围大，不宜遍历
我们采用的方式是：先求出数n的质因数，n的质因数的倍数显然不与n互质。
故我们只需要求出在区间a-b中与n不互质的数的数量x 最后b-a-x 就是我们要求的a-b区间内互质的数。其中，a-b中与n不互质的数可先求出 1-b的不互质数再求出1-（a-1）的不互质数相减得到。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
const int N =10010;
typedef long long ll;
using namespace std;

ll p[N],k;

void get(ll n) //得出n的质因数的个数
{
	k=0;
	for(int i=2;i<=n/i;i++)
	{
		if(n%i==0)
		{
			p[k++]=i;
			while(n%i==0)
			n/=i;
		}
	}
	if(n>1) p[k++]=n; 
}
ll no(ll m)
{
	ll top=0,sum=0;
	ll que[N];
	que[top++]=-1;//队列数组保存n所有质因子任意不相同组合的乘积 
	for(int i=0;i<k;i++)//共k个
	{
		ll t=top;
		for(int j=0;j<t;j++)
			que[top++]=que[j]*p[i]*(-1);  //奇加偶减 que[0]=-1，每次进入循环的第一个第一个pue[top]存的是一个单独质因子，接下来的pop[top]存的是质因子之间的任意倍数。加负号是为了在最后一步求和时满足且方便奇加偶减。
	}
	for(int i=1;i<top;i++)
		sum+=m/que[i];
		return sum;
}
int main()
{
	ll n,t,a,b,k=1;
	scanf("%lld",&t);
	while(t--)
	{
		scanf("%lld %lld %lld",&a,&b,&n);
		get(n);
		printf("Case #%lld: %lld\n",k++,b-no(b)-(a-1-no(a-1)));
	}
	return 0;
}