HDU 2204 Eddy's爱好

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2204

题意：求n范围内，可以表示成x^k(k>1)的数的个数。

思路：枚举指数k，那么n范围内有pow(n,1.0/k)-1个答案（1^k不算），但是也存在重复情况 x^14 = ( x^2 )^7 = ( x^7 )^2 ，分别用2和7计算时，有些数
会被计算两次，所以用容斥定理解决。2^60>10^18,我们先求60之内的素数，然后dfs容斥计算，如果指数乘起来达到60及其以上就退出。

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <utility>
using namespace std;

#define rep(i,j,k) for (int i=j;i<=k;i++)
#define Rrep(i,j,k) for (int i=j;i>=k;i--)

#define Clean(x,y) memset(x,y,sizeof(x))
#define LL long long
#define ULL unsigned long long
#define inf 0x7fffffff
#define mod %100000007

LL n;
LL ans;

int prime[60];
int sum = 0;
bool check(int x)
{
    for(int i = 2; i * i <= x; i++ )
        if ( x % i == 0 ) return false;
    return true;
}

void init()
{
    sum = 0;
    rep(i,2,60)
        if ( check(i) )
        {
            prime[++sum] = i;
        }
}


void dfs( int now , int cnt , int x )
{
    if ( now > sum  || x >= 60 ) return;
    LL temp=(LL)pow(n,1.0/x);

    if ( temp > 1 )
    {
        if ( cnt & 1 ) ans += temp - 1;
        else ans -= temp - 1;
    }
    rep( i , now+1 , sum)
        dfs( i , cnt + 1 , x*prime[i]  );
}



int main()
{
    init();
    while( cin>>n )
    {
        ans = 1;
        rep(i,1,sum)
                dfs( i , 1 , prime[i] );
        cout<<ans<<endl;
    }
    return 0;
}