思路:枚举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; }