POJ2109-Power of Cryptography

全解题报告索引目录 -> 【北大ACM – POJ试题分类】

转载请注明出处：http://exp-blog.com/2018/06/15/pid-409/

-------------------------------------------------------------------------

1. /*
2. Author: Exp
3. Date: 2017-12-07
4. Code: POJ 2109
5. Problem: Power of Cryptography
6. URL: http://poj.org/problem?id=2109
7. */
9. /*
10. 【题意分析】
11. 有指数函数 k^n = p ，
12. 其中k、n、p均为整数且 1<=k<=10^9 , 1<=n<= 200 , 1<=p<10^101
13. 给定 n 和 p ，求底数 k
15. 【解题思路】
16. 考虑到数值存储问题和精度问题，这题最直观的思路应该是使用 高精度算法 求解。
17. 而事实上，这题也可用公式法求解，但需要一些技巧。
20. 开方公式：k = n-sqrt(p)
21. 但C++的数学函数库并没有提供k次方的开方函数，此时需要转换一下公式：
22. k = p^(1/n)
24. 对p开k次方等价于求p的1/k次方，此时我们就可以用pow函数求解了：
25. k = pow(p, 1.0/n)
28. 其实严格来说，如果这题没有限制 底数k 是整数，就不可能通过公式投机取巧。
29. 简单来说，如果要使用公式法，那么题目中所有运算都只能基于double类型进行（int会溢出）
31. double的取值范围为10^(-307)~10^308，但小数精度只有前16位（可自行搜索double的精度丢失问题）.
32. 也是就说，当我们用double存储p的时候, 它就已经开始出现误差, 其误差范围在10^(-15)的数量级左右.
34. 此时套用公式对p开n次方根，须知开方运算是不会扩大误差范围的，
35. 所以 n-sqrt(p) 的小数位误差范围依旧在10^(-15)的数量级以内，
36. 又因为 k = n-sqrt(p) ，亦即计算所得的 n 的小数位误差范围也在10^(-15)的数量级以内,
37. 显然这个误差级数仅会对n的小数部分存在影响，四舍五入后对整数部分是无影响的.
38. 而题目已经限定了，n、k、p均是整数，因此使用公式法可以直接得到准确结果.
40. 假若题目不存在整数限制，当n极大时，k会极小（无限迫近1，对小数精度极高），
41. 此时公式法则会因为精度问题而失效.
43. */
45. #include <math.h>
46. #include <iostream>
47. using namespace std;
49. int main(void) {
50. double n , p;
51. while( cin >> n >> p) {
52. double tmp = pow(p, 1 / n); // p开n次方
53. int k = floor(tmp + 0.5); // 四舍五入（+0.5后向下取整）
54. cout << k << endl;
55. }
56. return 0;
57. }

不要惊讶！程序就是这么短，这就是“技巧”与“算法”的差别

用double避开高精度算法，可以说是这题最大的BUG 

有兴趣的同学也可以去挑战 二分+高精