Eddy's爱好
Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 2162 Accepted Submission(s): 939
Problem Description
Ignatius 喜欢收集蝴蝶标本和邮票,但是Eddy的爱好很特别,他对数字比较感兴趣,他曾经一度沉迷于素数,而现在他对于一些新的特殊数比较有兴趣。
这些特殊数是这样的:这些数都能表示成M^K,M和K是正整数且K>1。
正当他再度沉迷的时候,他发现不知道什么时候才能知道这样的数字的数量,因此他又求助于你这位聪明的程序员,请你帮他用程序解决这个问题。
为了简化,问题是这样的:给你一个正整数N,确定在1到N之间有多少个可以表示成M^K(K>1)的数。
这些特殊数是这样的:这些数都能表示成M^K,M和K是正整数且K>1。
正当他再度沉迷的时候,他发现不知道什么时候才能知道这样的数字的数量,因此他又求助于你这位聪明的程序员,请你帮他用程序解决这个问题。
为了简化,问题是这样的:给你一个正整数N,确定在1到N之间有多少个可以表示成M^K(K>1)的数。
Input
本题有多组测试数据,每组包含一个整数N,1<=N<=1000000000000000000(10^18).
Output
对于每组输入,请输出在在1到N之间形式如M^K的数的总数。
每组输出占一行。
每组输出占一行。
Sample Input
10 36 1000000000000000000
Sample Output
4 9 1001003332
Author
Eddy
Recommend
lcy
容斥原理求解,x^p<=n,开根号求大于1的x的个数,对于指数使用容斥原理,注意精度问题。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<queue>
#include<vector>
#include<cmath>
#include<string>
using namespace std;
typedef long long LL;
const int MAXN=1000008;
int prime[MAXN];
bool vis[MAXN];
LL ans,n;
int tot;
int getPrime(int n){
tot=0;
memset(vis,0,sizeof(vis));
for(int i=2;i<=n;i++){
if(!vis[i]){
prime[tot++]=i;
}
for(int j=0;j<tot&&i*prime[j]<=n;j++){
vis[i*prime[j]]=true;
if(i%prime[j]==0){
break;
}
}
}
return tot;
}
void dfs(int dep,int mul,int num){
for(int i=num;i<tot&&prime[i]*mul<=(int)(log(n)/log(2.0));i++){
int tp=mul*prime[i];
if(dep&1){
ans-=(LL)(pow(n,1/(double)tp)+0.000001)-1;//1不算,注意精度
}else{
ans+=(LL)(pow(n,1/(double)tp)+0.000001)-1;
}
dfs(dep+1,tp,i+1);
}
}
int main(){
getPrime(1000000);
while(~scanf("%I64d",&n)){
ans=1;//1满足条件
dfs(0,1,0);
printf("%I64d\n",ans);
}
return 0;
}