题意:
给出一个范围(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;
}