HDU4430 Yukari's Birthday(枚举+二分)

Yukari's Birthday  HDU4430

 

就是枚举+二分：

注意处理怎样判断溢出。。。（因为题目只要10^12）

先前还以为要用到快速幂和等比数列的快速求和（但肯定会超__int64）

而且这样判断会超时的。。。

还有题目中的And it's optional to place at most one candle at the center of the cake. （中间的蜡烛可有可无）

还有观察数据就知道：因为n最大10^12，r最多枚举到40，然后二分k的结果，看是否有符合条件的k存在。

#include <stdio.h>
#include <iostream>
#include <algorithm>
using namespace std;
const __int64 N=1000010;
__int64 fun1(__int64 k,__int64 r,__int64 n)
{
    __int64 s=0,t=k;
    for(int i=1;i<=r;i++)
    {
        s+=t;
        if(s > n)//判断在求和的过程中是否已经比n大（可省时很多且不会超__int64）
            return N*N;
        t*=k;
    }
    return s;
}
__int64 fun(__int64 r,__int64 n)
{
    if(r == 1)
        return n;
    __int64 mid,p,lef=2,rig=N;
    while(lef <= rig)
    {
        mid=(lef+rig)/2;
        p=fun1(mid,r,n);
        if(p == n)
            return mid;
        if(n > p)
            lef=mid+1;
        else
            rig=mid-1;
    }
    return -1;
}
int main()
{
    __int64 n,p,r,rr,rp;
    while(scanf("%I64d",&n)!=EOF)
    {
        rr=N;rp=N;
        for(r=1;r<=40;r++)//只要到40就行了
        {
            p=fun(r,n);//中间放蜡烛
            if(p!=-1)
                if(r*p < rr*rp || (r*p == rr*rp && r < rr))
                {
                    rr=r;
                    rp=p;
                }
            p=fun(r,n-1);//中间不放蜡烛
            if(p!=-1)
                if(r*p < rr*rp || (r*p == rr*rp && r < rr))
                {
                    rr=r;
                    rp=p;
                }
        }
        printf("%I64d %I64d
",rr,rp);
    }
    return 0;
}