hdu4430 Yukari's Birthday(二分)

题意：
给定一个n，计算满足x^1+x^2+……+x^r=n且使x*r最小（如果有多个，取r最小的那个）的x，r值。
由于2^40大于1e12，所以我们可以枚举r从1到40，但这样可能会超时，所以可以根据3^26大于2e12，我们可以枚举r从1到26，并单独计算x=2的情况。
对于给定的每个r，x^1+x^2+……+x^r这个式子就关于x单调递增，所以就可以用二分来做。

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <stdlib.h>
#include <stack>
#include <vector>
#include <math.h>
#include <string.h>
#include <queue>
#define msc(X) memset(X,-1,sizeof(X))
#define ms(X) memset(X,0,sizeof(X))
typedef long long LL;
using namespace std;
LL res,p2;
int p1;
void cal(LL n)
{
    if(((n+2)&(n+1))==0){
            int k=1;
            for(int i=2;i<=40;i++)
                if(!(n&(1ll<<i)))
                    {
                        k=i-1;
                        break;
                    }
            if(res>=(k<<1)) {
                res=(k<<1);
                p1=k;
                p2=2;
            }
        }
    for(int i=1;i<27;i++)
    {
        LL l=(LL)pow(2.0*n,1.0/(i+1.0)),
        r=(LL)pow((double)n,1.0/i)+5ll;
        while(l<r){
            LL m=(l+r)>>1;
            LL tmp=m;
            for(int j=1;j<i;j++)
                tmp=tmp*m+m;
            if(tmp<n) l=m+1;
            else if(tmp>n) r=m;
            else{
                if(res>m*i){
                    res=m*i;
                    p1=i;
                    p2=m;
                }
                break;
            }
        }
    }
}
int main(int argc, char const *argv[])
{
    LL n;
    while(scanf("%I64d",&n)!=EOF){
        res=1ll<<40;
        cal(n);
        cal(n-1);
        printf("%d %I64d\n",p1,p2 );
    }
    return 0;
}