hdu 2204 http://acm.hdu.edu.cn/showproblem.php?pid=2204

//直接模拟容斥原理的过程解决:(纠结了一天 忘了减去1) 把1去掉 最后再加上。
//1:  10^18<2^60;
//2:  为了避免运算重叠可以只利用质数{2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59};
//3:  利用容斥原理要注意(-1)^n求交
#include<stdio.h>
#include<math.h>
__int64 ans;
__int64 test;
int num[19]={2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59};
void inclu(int judge,int cnt,int pox )
{
    __int64 tem;
    for(int i=pox;i<19;i++)
    {
        tem=(__int64)(pow(test,1.0/(cnt*num[i])));//利用容斥原理
        tem-=1; // 1 的情况。
        if(tem<1)
          break;
        judge%2?ans+=tem:ans-=tem; //求交 一个一个求 
        inclu(judge+1,cnt*num[i],i+1);
    }
}
int main()
{
    while(scanf("%I64d",&test)!=EOF)
    {
        ans=0;
        inclu(1,1,0);
        printf("%I64d\n",ans+1);
    }
  return 0;
}

你可能感兴趣的:(hdu 2204 http://acm.hdu.edu.cn/showproblem.php?pid=2204)