POJ解题报告——1001 Exponentiation
问题描述
求高精度幂,输入是一个实数 R (0.0 < R < 99.999) 和 指数 n (0 < n <= 25),
要求精确计算 R 的 n 次方(Rn)。
如 95.123 12
结果为548815620517731830194541.899025343415715973535967221869852721
详见 http://poj.org/problem?id=1001
解题思路
题目中实数的位数是给定的,如果把实数当作整数来计算,最后确定结果的位数,求高精度幂就相当于基本的大整数乘法。
在知道小数点有 m 位的情况下,小数部分一定有 m*n 位(输出时去掉最后的 0 )。
具体需要在读入的时候记录小数点的位置,并将实数 R 转换成 整数 r 。
循环乘数并进位。
最后依次输出整数部分,小数点和小数部分。
题目其实不难,但一直答案错误,上网查了下,发现一组测试数据是10.000 10,答案是10
题目要求整数不输出小数点,估计大部分错的人是没看到这句话。
代码如下
输出的时候,判断下最后一个非零数在哪,如果不在小数部分,小数点‘.’不要输出
代码以0为低位,size()为高位,计算过程不区分小数和整数部分。
因为一开始没想到,有小数点但结果是整数的情况,代码写的不够连贯。
#include <iostream>
#include <vector>
#include <string>
using namespace std;
int main()
{
string R;
int n,m;
while (cin >> R >> n) {
vector<int>a;
a.push_back(1);
int r = 0;
m = 0;
for (int i = 0; i < R.length(); ++i) {
if (R[i] == '.')
m = R.length() - i - 1;
else
r = r * 10 + R[i] - '0';
}
for (int i = 0; i < n; ++i) {
for (int j = 0; j < a.size(); ++j)
a[j] *= r;
int k = 0;
while (a[k] > 9 || (k + 1 < a.size())) {
if (k + 1 == a.size())
a.push_back(0);
a[k + 1] += a[k] / 10;
a[k] = a[k] % 10;
k++;
}
}
m = m * n;
while (a.size() < m)
a.push_back(0);
for (int i = a.size() - 1; i >= m; --i)
cout << a[i];
int set0 = 0;
while (a[set0] == 0) ++set0;
if (m > set0)cout << ".";
for (int i = m - 1; i >= set0; --i)
cout << a[i];
cout << endl;
}
return 0;
}
测试数据
为了方便只看到这篇文章的人测试。这里给出POJ全部的实例数据和蛮坑的那个测点
ps,第4行9和第7行1前面还有一个空格。(如果你的代码完全以字符串读入,小心别看错了。)
Sample Input
95.123 12
0.4321 20
5.1234 15
6.7592 9
98.999 10
1.0100 12
10.000 1
Sample Output
548815620517731830194541.899025343415715973535967221869852721
.00000005148554641076956121994511276767154838481760200726351203835429763013462401
43992025569.928573701266488041146654993318703707511666295476720493953024
29448126.764121021618164430206909037173276672
90429072743629540498.107596019456651774561044010001
1.126825030131969720661201
10