POJ解题报告——1001 Exponentiation

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