关闭

hdu4430 Yukari's Birthday(二分)

110人阅读 评论(0) 收藏 举报
分类:

题意:
给定一个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;
}
0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:17486次
    • 积分:1209
    • 等级:
    • 排名:千里之外
    • 原创:108篇
    • 转载:4篇
    • 译文:0篇
    • 评论:5条
    最新评论