POJ1001

来源:互联网 发布:监控软件流量 编辑:程序博客网 时间:2024/06/02 01:21

求高精度幂
Time Limit: 500MS
Memory Limit: 10000K
Total Submissions: 167412
Accepted: 40583
Description
对数值很大、精度很高的数进行高精度计算是一类十分常见的问题。比如,对国债进行计算就是属于这类问题。

现在要你解决的问题是:对一个实数R( 0.0 < R < 99.999 ),要求写程序精确计算 R 的 n 次方(Rn),其中n 是整数并且 0 < n <= 25。
Input
T输入包括多组 R 和 n。 R 的值占第 1 到第 6 列,n 的值占第 8 和第 9 列。
Output
对于每组输入,要求输出一行,该行包含精确的 R 的 n 次方。输出需要去掉前导的 0 后不要的 0 。如果输出是整数,不要输出小数点。
Sample Input
95.123 12
0.4321 20
5.1234 15
6.7592 9
98.999 10
1.0100 12
Sample Output
548815620517731830194541.899025343415715973535967221869852721
.00000005148554641076956121994511276767154838481760200726351203835429763013462401
43992025569.928573701266488041146654993318703707511666295476720493953024
29448126.764121021618164430206909037173276672
90429072743629540498.107596019456651774561044010001
1.126825030131969720661201
Source
East Central North America 1988
Translator
北京大学程序设计实习,Xie Di

代码很简单就不解释了

#include <iostream>#include <cstring>using namespace std;#define size 10000int main(){    char c[size];    int n;    while (cin >> c >> n) {        int i, j, k;        int arr[size] = {0};        int len = strlen(c);        int mul;        int cry;        int dp;        int tmp;        // 处理无效0 ascii('0') == 48        for (i = len - 1; c[i] == 48; i--) {            len--;        }        c[len] = 0;        // 计算乘数 小数点位置        dp = 0;        mul = 0;        for (i = 0; i < len; i++) {            if (c[i] == '.') {                dp = len - i - 1;            }            else {                mul = mul * 10 + c[i] - 48;            }        }        dp = dp*n;        // 初始化乘积数组,逆序放入        tmp = mul;        k = 0;        while (tmp) {            arr[k++] = tmp % 10;            tmp /= 10;        }        // 连乘        for (i = 1; i < n; i++) {            // 初始化进位            cry = 0;            for (j = 0; j < k; j++) {                arr[j] = arr[j] * mul + cry;                cry = arr[j] / 10;                arr[j] %= 10;            }            while (cry) {                arr[k++] = cry%10;                cry /= 10;            }        }        if (k < dp)            k = dp;        while (k>0) {            if (k == dp)                cout << '.';            cout << arr[--k];        }        cout << endl;    }    return 0;}
0 0