hdu 2204 Eddy's爱好

// 一个整数N，1<=N<=1000000000000000000(10^18).
// 输出在在1到N之间形式如M^K的数的总数
// 容斥原理
// 枚举k=集合{2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59}中选取1或多个数乘积的值 
而且大于三个质数的乘积就大于60 而2^60>10^18 所以最多只取3个素数

#include <iostream>
#include <algorithm>
#include <queue>
#include <vector>
#include <math.h>
#include <stdio.h>
#include <string.h>
using namespace std;
#define maxm 100010
#define maxn 1000110
int prim[]={,,,,,,,,,,,,,,,,};// 17个
int main()
{
    int i,j,ct,mid;
    int rc[];
    __int64 n,m,ans,tp;
    while(scanf("%I64d",&n)!=EOF){
        ans=;
        // if(n<4) { printf("%d\n",1);continue;}
        int x=(<<)-;
        for(i=;i<=x;i++){
            mid=i;
            ct=j=;
            while(mid){
                if(mid&) rc[ct++]=j;
                j++;
                if(ct>) break;
                mid=mid>>;
            }
            if(ct>) continue;
            tp=;
            for(j=;j<ct;j++)
                tp=tp*prim[rc[j]];
            if(tp>=) continue;
            m=pow(n+1.0,1.0/tp);
            m++;
            while(pow(m,0.0+tp)>n) m--;
            if(ct%) ans+=(m-);
            else ans-=(m-);
        }
         printf("%I64d\n",ans+);
    }
    return ;
}