容斥定理模板 HDU 4135 Co-prime

传送门

:解决[x,y]区间与n互质的数的个数->总个数-不互质(能被n素因子整除的数) 

#include<bits/stdc++.h>
#define il inline
#define pb push_back
#define fi first
#define se second
#define ms(_data,v) memset(_data,v,sizeof(_data))
#define sc(n) scanf("%d",&n)
#define SC(n,m) scanf("%d %d",&n,&m)
#define SZ(a) int((a).size())
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define drep(i,a,b)	for(int i=a;i>=b;--i)
using namespace std;
typedef long long ll;
const ll inf=0x3f3f3f3f;
const double PI=acos(-1.0);
const double eps=1e-9;
//解决[x,y]区间与n互质的数的个数->总个数-不互质(能被n素因子整除的数) 
ll a[1000],num,x,y;//a存质数
il void init(ll n) { //求一个数的质因子
	num=0;
	for(int 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;
}
ll que[10005];
il ll solve() { //容斥模板 
	ll k,t=0,sum=0;
	que[t++]=-1; //que[]存组合，并且奇加偶减 
	for(int i=0; i<num; i++) {
		k=t;
		for(int j=0; j<k; j++)	que[t++]=que[j]*a[i]*(-1);
	}
	ll nx=0,ny=0;
	for(int i=1; i<t; i++){
		nx+=(x-1)/que[i];
		ny+=(y)/que[i];
	}
	return (y-x+1)-(ny-nx);
}
int main() {
	ll T,n,sum;
	while(scanf("%lld",&T)!=EOF) {
		for(int i=1; i<=T; i++) {
			scanf("%lld%lld%lld",&x,&y,&n);
			init(n);
			printf("Case #%lld: %lld\n",i,solve());
		}
	}
	return 0;
}