hdu 4135 Co-prime

一个容斥定理的运用，代码实现感觉有点难理解，不过自己手动模拟一下程序就好了

#include <iostream>
#include <cmath>
#include <cstdio>
using namespace std;

typedef long long ll;

int fac[100000];//存的是n的素因子
ll Q[100000];
int num_fac;
void get_factor(ll n)
{
num_fac=0;
for(int i=2;i<=sqrt(n);i++)
{
if(n%i==0)
{
fac[num_fac++]=i;
while(n%i==0) n/=i;
}
}
if(n>1) fac[num_fac++]=n;
}
ll rc(ll n)
{
ll k,t,ans;
ans=t=0;
Q[t++]=-1;//
for(ll i=0;i<num_fac;i++)
{
k=t;
for(ll j=0;j<k;j++)
Q[t++]=(-1) * Q[j] * fac[i];
}
for(ll i=1;i<t;i++)
ans += n/Q[i];
return ans;
}
int main()
{
int t,cs=1;
scanf("%d",&t);
while(t--)
{
ll a,b,n;
scanf("%lld %lld %lld",&a,&b,&n);
get_factor(n);
ll ans;
ans= (b-rc(b)) - (a-1-rc(a-1)) ;//[a,b]
printf("Case #%d: %lld\n",cs++,ans);
}
return 0;
}
 
---------------------