Eddy's爱好

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 1266    Accepted Submission(s): 541


Problem Description
Ignatius 喜欢收集蝴蝶标本和邮票,但是Eddy的爱好很特别,他对数字比较感兴趣,他曾经一度沉迷于素数,而现在他对于一些新的特殊数比较有兴趣。
这些特殊数是这样的:这些数都能表示成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
 
 1 #include <iostream>
 2 #include <string.h>
 3 #include <stdio.h>
 4 #include <math.h>
 5 using namespace std;
 6 int a[]= {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59};
 7 int b[]= {6,10,14,15,21,22,26,33,34,35,38,39,46,51,55,57,58};
 8 long long n,an;
 9 void dfs()
10 {
11     int i;
12     for(i=0; i<17; i++)
13     {
14         long long h=pow(n,1.0/a[i]);
15         if(pow(h+1,a[i])<=n)h++;
16         if(h<2)break;
17         an+=h-1;
18     }
19     for(i=0; i<17; i++)
20     {
21         long long h=pow(n,1.0/b[i]);
22         if(pow(h+1,b[i])<=n)h++;
23         if(h<2)break;
24         an-=h-1;
25     }
26     if(n-pow(2.0,30)>=0) an++;
27     if(n-pow(2.0,42)>=0) an++;
28     if(n-pow(3.0,30)>=0) an++;
29     an++;
30 }
31 int main()
32 {
33     while(cin>>n)
34     {
35         an=0;
36         dfs();
37         cout<<an<<endl;
38     }
39 }
View Code

 

posted on 2014-03-28 20:45  ERKE  阅读(162)  评论(0编辑  收藏  举报