【POJ】poj 1001（高精度）

 

Description

Problems involving the computation of exact values of very large magnitude and precision are common. For example, the computation of the national debt is a taxing experience for many computer systems. 

This problem requires that you write a program to compute the exact value of Rn where R is a real number ( 0.0 < R < 99.999 ) and n is an integer such that 0 < n <= 25.

Input

The input will consist of a set of pairs of values for R and n. The R value will occupy columns 1 through 6, and the n value will be in columns 8 and 9.

Output

The output will consist of one line for each line of input giving the exact value of R^n. Leading zeros should be suppressed in the output. Insignificant trailing zeros must not be printed. Don't print the decimal point if the result is an integer.

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

这是一个高精度；可以通过模拟手算的方法来得出真确的答案，如果数据量大一点的话，也可以通过快速幂的方法来计算，不过这对高进度写的要求比较高，不过细心一点没有问题： 详细注释见代码：

代码：

 

/**********************************************************
    @auther Thunder01
    poj 1001
    算法： 高进度，如果数据强一点的话就应该用快速幂；
    ps：先看MAIN函数
    http://paste.ubuntu.com/8299713/
*********************************************************/
#include <iostream>///c++ 的输入输出，貌似没用
#include <cstdio>///和stdio.h 一样
#include <cstring>///调用了memset这个函数
#include <string>///没用这个头文件
#include <cmath>///数学库
#define MAX 1000
using namespace std; /// c++ 的命名空间，c的不用写，可以无视
int ans[MAX];
int a[MAX],b[MAX];
int lena, lenb;
void product()
{
    /**
    int temp;
    memset(ans,0,sizeof(ans));
    for ( int i = 0; i < lena; i++ )
    {
        temp = 0;
        for ( int j = 0; j < lenb; j++ )
        {
            ans[i+j] += temp + (a[i]) * ( b[j]);
            temp = ans[i+j] / 10;
            ans[i+j] = ans[i+j] % 10;
        }
        ans[i+lenb] += temp;
    }
    if (ans[lena+lenb-1] != 0 ) lena = lena + lenb;
    else lena = lena + lenb -1;
    for ( int i = 0; i < lena; i++ ) a[i] = ans[i];
    */
    int i, j;
    memset(ans, 0, sizeof(ans));///暂时储存结果，然后赋值给a[],如果直接变化a[]的话就会产生不必要的错误
    for (i=1; i<=lena; i++) ///笔算乘法
    {
        for (j=1; j<=lenb; j++)
        {
            ans[i+j-1] += a[i] * b[j];/// 逐位乘法
            if (ans[i+j-1]>9) ///大于9的话就向下一位进，这个是笔算，自己模拟一边就可以
            {
                ans[i+j] += ans[i+j-1] / 10;
                ans[i+j-1] %= 10;
            }
        }
    }
    if (ans[lena+lenb-1]>9)
    {
        ans[lena+lenb] += ans[lena+lenb-1] / 10;
        ans[lena+lenb-1] %= 10;
    }
    lena = lena + lenb;///乘法之后长度相加，无视前导0，后面会去掉
    for (i=1; i<=lena; i++) a[i] = ans[i];///答案赋值给a[],b[]始终是题目给的数据
}
int main()
{
    ///freopen("in.txt","r",stdin);///文件读入，无视
    int n;
    char str[6];
    while ( scanf("%s %d",str,&n) != EOF )///首先是输入；
    {
        memset(ans,0,sizeof(ans));
        int point = -1;
        int j,i;
        for ( i = 5, j = 1; i >=0 ; i--)///将数据存进数组，然后判断是否有小数点；注意这里是倒着存的。
            ///比如1.2345存进去之后就是54321了
        {
            if (str[i]!='.') a[j] = b[j++] = str[i] - '0';/// 没有小数点的话就是将字符转化位数字存入数组，
            ///注意这里将小数点省去了
            else point = i;///存小数点的位置比如1.2345, 那么point就是 1 模拟一下就可以明白，
            ///那么这个数字的小数的个数是5-point个；这里不考虑末尾是否有很多0
        }
        if ( point == -1 ) lena = lenb = 6;///如果没有小数的话得到的就是一个六位数
        else lena = lenb = 5;///有的话就是五位数
        for ( int i = 1; i < n; i++) product();///进行n次乘法
        if ( point == -1 )///没有小数点直接输出
        {
            for ( int i = lena-1; i >= 0; i--)///注意这里是倒着输出的，刚才是倒着存入数组的
            {
                printf("%d",a[i]);
            }
            printf("\n");
        }
        else
        {
            point = 5 - point;///小数的位数
            point *= n;///n次乘法之后的小数个数
            int down , up ;
            for (i=1; i<=lena; i++)
            {
                if (a[i] != 0)
                {
                    down = i;///小数点后面的0，注意这里是倒着存的所以前面的是小数，后面的是整数部分
                    break;
                }
            }
            for (j=lena; j>=1; j--)
            {
                if (a[j]!=0)
                {
                    up = j;///整数前面的0；
                    break;
                }
            }
            i = up;
            if (up<point) i = point;///看小数点的位置，也就是说去除小数点后面的0的时候不要取点小数点前面的0；
            ///也就输说如果是12340.00 这个的答案应该是1240,不能把后面的0全部取点
            j = down;
            if (j>point) j = point + 1;///同上
            for (; i>=j; i--)
            {
                if (i==point) printf(".");
                printf("%d", a[i]);
            }
            printf("\n");
        }
    }
}