hdu 4235 容斥原理模板题

#include<cstdio>
#include<cstdlib>
#include<cstring>

using namespace std;

typedef long long ll;
const int maxn = 1000000+10;
ll que[maxn],fac[maxn],num;

void divid(ll n)
{
    num=0;
    for(ll i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            while(n%i==0) n/=i;
            fac[num++]=i;
        }
    }
    if(n!=1) fac[num++]=n;
}
ll solve(ll n)
{
    int k,t=0;
    ll ans=0;
    que[t++]=-1;
    for(ll i=0;i<num;i++)
    {
        k=t;
        for(ll j=0;j<k;j++)
        {
            que[t++]=-1*que[j]*fac[i];
        }
    }
    for(ll i=1;i<t;i++) ans+=n/que[i];
    return ans;
}


int main()
{
    int T;
    scanf("%d",&T);
    int cas=0;
    while(T--)
    {
        ll a,b,n;
        scanf("%I64d%I64d%I64d",&a,&b,&n);
        divid(n);
        ll res=b-solve(b)-(a-1-solve(a-1));
        printf("Case #%d: %I64d\n",++cas,res);
    }
    return 0;
}