hdu 2204 Eddy's爱好 (数论，容斥原理)

题意：

给出一个范围(1,N)问存在多少这样的数M^K(k>1)

题解：

我这样分析：

1、对于某个M^K,我们不妨将K分解，假设分解成为a,b,那么M^(a*b)，那么会发现一重复的问题，就是(M^a)^b和(M^b)^a是重复计算了，那么要使得不重复计算K一定不能分解，那么K必须为素数。

2、现在分析某个M^K,对于这样的K有多少个M能够满足，满足的个数为：N^(1/K),好啦，分析到这里可以通过枚举素数计算个数了，因为10^18小于2^60，那么也就是说只要枚举60以内的素数。

3、上面的分析为基础，继续分析会发现这样的特例对于a^i,b^j,如果a=t^j,b=t^i，那么这样就重复计算了，所以就要用到容斥原理把重复的删。对于这样的a^i,b^j,重复的个数就是N^(1/i/j);

这里注意一个点，如果直接这样计算，我保证会出现负数等非答案的解，因为精度问题。有一个技巧可以解决这个问题。就是在pow后面加一个esp然后在强转成long long，这样就不会出错了。

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
void cmax(int& a,int b){ if(b>a)a=b; }
void cmin(int& a,int b){ if(b<a)a=b; }
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=10007;
const int maxn=110;
const double eps=1e-8;
int p[]={2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61};
double n;

ll dfs(int s,ll x){
    ll res=0;
    for(int i=s;i<19;i++){
        ll cnt=(ll)(pow(n,1.0/(x*p[i]))+eps);
        res+=cnt-dfs(i+1,x*p[i]);
    }
    return res;
}

int main(){
    //freopen("E:\\read.txt","r",stdin);
    while(scanf("%lf",&n)!=EOF){
        printf("%I64d\n",dfs(0,1));
    }
    return 0;
}