hdu4135Co-prime 容斥原理水题

//问一个区间[a,b]与n互素的数的个数
//利用容斥原理可知
//在[a,b] 区间内对n的素数因子
//ans = 被一个数整除的数的个数 - 被两个数的最小公倍数整除的数的个数 + 被三个数的。。。
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std ;
const int maxn = 100010 ;
typedef __int64 ll ;
ll p[maxn] ;int len ;
void get_prime(ll n)
{
    len = 0 ;
    for(ll i = 2;i*i <= n;i++)
    {
        if(n%i == 0)p[++len] = i ;
        while(n%i == 0)n/=i ;
    }
    if(n>1)p[++len] = n;
}
ll dfs(int pos , ll n)
{
    ll ans = 0 ;
    for(int i = pos ;i <= len ;i++)
        ans += n/p[i] - dfs(i+1 , n/p[i]) ;
    return ans ;
}
int main()
{
    ll a , b ,n ;
    int T ;
    int cas =  0 ;
    scanf("%d" ,&T) ;
    while(T--)
    {
        scanf("%I64d%I64d%I64d" , &a , &b , &n);
        get_prime(n) ;
        ll ans = (b - dfs(1 , b)) - (a  - 1 - dfs(1 , a-1)) ;
        printf("Case #%d: " ,++cas) ;
        printf("%I64d " , ans) ;
    }
    return  0 ;
}