hdoj 4282 A very hard mathematic problem

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4282

解题思路：二分

///////////////////////////////////////////////////////////////////////////
//problem_id: hdoj 4282
//user_id: SCNU20102200088
///////////////////////////////////////////////////////////////////////////

#include <algorithm>
#include <iostream>
#include <iterator>
#include <iomanip>
#include <cstring>
#include <cstdlib>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cmath>
#include <queue>
#include <stack>
#include <list>
#include <set>
#include <map>
using namespace std;

///////////////////////////////////////////////////////////////////////////
#pragma comment(linker,"/STACK:1024000000,1024000000")

#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
///////////////////////////////////////////////////////////////////////////

///////////////////////////////////////////////////////////////////////////
const double EPS=1e-8;
const double PI=acos(-1.0);

const int x4[]={-1,0,1,0};
const int y4[]={0,1,0,-1};
const int x8[]={-1,-1,0,1,1,1,0,-1};
const int y8[]={0,1,1,1,0,-1,-1,-1};
///////////////////////////////////////////////////////////////////////////

///////////////////////////////////////////////////////////////////////////
typedef long long LL;

typedef int T;
T max(T a,T b){ return a>b? a:b; }
T min(T a,T b){ return a<b? a:b; }
T gcd(T a,T b){ return b==0? a:gcd(b,a%b); }
T lcm(T a,T b){ return a/gcd(a,b)*b; }
///////////////////////////////////////////////////////////////////////////

///////////////////////////////////////////////////////////////////////////
//Add Code:
LL power(LL a,LL n){
    LL ret=1;
    while(n){
        if(n&1) ret*=a;
        a*=a;
        n>>=1;
    }
    return ret;
}

LL cal(LL x,LL y,LL z){
    return power(x,z)+power(y,z)+x*y*z;
}

LL erfen(LL x,LL z,LL k){
    LL l=x+1,r=(LL)pow(k,1.0/z);
    while(r>=l){
        LL m=(l+r)/2,temp=cal(x,m,z);
        if(temp==k) return 1;
        if(temp>k) r=m-1;
        else l=m+1;
    }
    return 0;
}
///////////////////////////////////////////////////////////////////////////

int main(){
    ///////////////////////////////////////////////////////////////////////
    //Add Code:
    LL k;
    while(scanf("%I64d",&k)!=EOF){
        if(k==0) break;
        LL ans=0;
        for(LL z=2;z<32;z++){
            for(LL x=1;power(x,z)<k/2;x++) ans+=erfen(x,z,k);
        }
        printf("%I64d
",ans);
    }
    ///////////////////////////////////////////////////////////////////////
    return 0;
}

///////////////////////////////////////////////////////////////////////////
/*
Testcase:
Input:
9
53
6
0
Output:
1
1
0
*/
///////////////////////////////////////////////////////////////////////////