思路:枚举r+二分k
分析:
1 题目要求的是找到一组最小的r*k,如果r*k相同那么就找r最小的。
2 很明显k>=2,根据n <= 10^12,那么可以知道r的最大值r<50,所以只要枚举枚举r的值,然后二分k的大小找到所有的解,存入一个结构体里面,然后在对结构体排序,那么这样就可以得到最后的ans
3 注意题目说了中心点最多一个蜡烛,所以写二分的时候应该注意判断的条件;
4 还有可能计算得到结果超了long long直接变成负数所以应该对或则个进行判断;
5 还有就是在二分的时候右边界right不能为n,应该是pow(n , 1.0/i);
6 还有hdu上面这一题不能用long long ,用__int64;
7 注意这种涉及到很多次求次方的问题,不要调用库函数自己写一个Pow函数。
代码:
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
typedef __int64 int64;
#define MAXN 1000
int64 n;
struct mulity{
int64 r;
int64 k;
int64 sum;
}m[MAXN];
/*排序函数*/
bool cmp(mulity m1 , mulity m2){
if(m1.sum < m2.sum)
return true;
else if(m1.sum == m2.sum && m1.r < m2.r)
return true;
return false;
}
/*Pow函数*/
int64 Pow(int64 mid , int64 n){
int64 ans = 1;
for(int i = 1 ; i <= n ; i++)
ans *= mid;
return ans;
}
int main(){
int i , j , cnt;
int64 right , left , mid;
while(scanf("%I64d" , &n) != EOF){
cnt = 0;
for(i = 1 ; i < 50 ; i++){/*枚举r的范围*/
right = pow(n , 1.0/i);
left = 2;
while(left <= right){
mid = (left+right)/2;
int64 tmp = 0;
for(j = 1 ; j <= i ; j++){
int64 tmp_sum = Pow(mid , j);
if(tmp_sum <= 0){/*如果返回的tmp_sum <= 0那么也就是说明超过了long long*/
tmp = -1;
break;
}
tmp += tmp_sum;
}
if(tmp == n-1 || tmp == n){/*如果计算的结果等于n或n-1*/
m[cnt].r = i;
m[cnt].k = mid;
m[cnt++].sum = i*mid;
break;
}
else if(tmp <= 0 || tmp > n)/*这里记得判断tmp<0的判断*/
right = mid-1;
else
left = mid+1;
}
}
sort(m , m+cnt , cmp);/*排序*/
printf("%I64d %I64d\n" , m[0].r , m[0].k);
}
return 0;
}