容斥A - Eddy's爱好 HDU - 2204

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

通过观察，可以发现若一个数可以表示成x^(k*t)，则可以表示成(x^k)^t。因此指数必然为素数。

然后对n开(x^k)^t次幂，求出这样得出来的数字是n以内有几个数能被某个数的(x^k)^t次幂表示。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<math.h>
#include<set>
#include<vector>
#include<sstream>
#include<queue>
#define ll long long
#define ull unsigned long long
#define PI 3.1415926535897932384626
#define inf 0x3f3f3f3f
using namespace std;
int prime[]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61};
ull ans,n;
void dfs(int v,int pos,ull num)
{
    ull temp;
    for(int i=v;i<18;i++)
    {
        temp=prime[i]*num;
        //cout<<"temp"<<temp<<endl;
        if(pos&1)
            ans+=(ull)(pow(n,1.0/temp));
        else
            ans-=(ull)(pow(n,1.0/temp));
        dfs(i+1,pos+1,temp);
    }
}
int main()
{
    while(scanf("%llu",&n)!=EOF)
    {
        ans=0;
        dfs(0,1,1);
        printf("%llu\n",ans);
    }
}