Eddy's爱好 hdu2204

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

 

#include <iostream>
#include <string.h>
#include <stdio.h>
#include <math.h>
using namespace std;
int a[]= {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59};
int b[]= {6,10,14,15,21,22,26,33,34,35,38,39,46,51,55,57,58};
long long n,an;
void dfs()
{
    int i;
    for(i=0; i<17; i++)
    {
        long long h=pow(n,1.0/a[i]);
        if(pow(h+1,a[i])<=n)h++;
        if(h<2)break;
        an+=h-1;
    }
    for(i=0; i<17; i++)
    {
        long long h=pow(n,1.0/b[i]);
        if(pow(h+1,b[i])<=n)h++;
        if(h<2)break;
        an-=h-1;
    }
    if(n-pow(2.0,30)>=0) an++;
    if(n-pow(2.0,42)>=0) an++;
    if(n-pow(3.0,30)>=0) an++;
    an++;
}
int main()
{
    while(cin>>n)
    {
        an=0;
        dfs();
        cout<<an<<endl;
    }
}