求高精度幂（Poj 1001)

Description

对数值很大、精度很高的数进行高精度计算是一类十分常见的问题。比如，对国债进行计算就是属于这类问题。 

现在要你解决的问题是：对一个实数R( 0.0 < R < 99.999 )，要求写程序精确计算 R 的 n 次方(Rⁿ)，其中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

#include <iostream>
#include<cstring>
using namespace std;
int a[111111],b[111111],temp[111111],la,lb;
void multi()
{
    memset(temp,0,sizeof(temp));
    int i,j;
    for(i=0;i<la;i++)
    {
        for(j=0;j<lb;j++)
        {
            temp[i+j]+=a[i]*b[j];//各位相加
            if(temp[i+j]>9)//如果相加大于9，进位
            {
                temp[i+j+1]+=temp[i+j]/10;
                temp[i+j]%=10;
            }
        }
    }
    if(temp[la+la-2]>9)//最高位大于9进位
    {
        temp[la+lb-1]+=temp[la+lb-2]/10;
        temp[la+lb-2]%=10;
    }
    la=la+lb;
    for(i=0;i<la;i++)//将得到的数拷入数组a中
    {
        a[i]=temp[i];
    }
}
int main()
{
    char s[10];
    int pos=-1,i,j,n;
    while(cin>>s>>n)
    {
        for(i=5,j=0;i>=0;i--)//判断小数点的位数
        {
            if(s[i]!='.')
            {
                a[j]=s[i]-'0';
                b[j]=a[j];//将个位上的数存入数组a与数组b中
                j++;
            }
            else
            {
                pos = i;//记录小数点的位数
            }
        }
        la=lb=(pos==-1?6:5);
        for(i=1;i<n;i++)//第一次运行R^2；所以只循环n-1次
        {
            multi();
        }
        if(pos==-1)//如果没有小数
        {
            for(i=la-1;i>=0;i--)
            {
                cout<<a[i];
            }
            cout<<endl;
        }
        else
        {
            int r=0,l=0;
            pos=5-pos;
            pos=pos*n;//小数的位数
            for(i=0;i<la;i++)//后消0
            {
                if(a[i]!=0)
                {
                    l=i;
                    break;
                }
            }
            for(i=la-1;i>=0;i--)//前消0
            {
                if(a[i]!=0)
                {
                    r=i;
                    break;
                }
            }
            if(r<pos)//如果计算的数小于1
            {
                r=pos-1;
            }
            if(l>pos)//如果所得数大于1
            {
                l=pos;
            }
            for(i=r;i>=l;i--)
            {
                if(i==pos-1)
                {
                    cout<<".";
                }
                cout<<a[i];
            }
            cout<<endl;
        }
    }
    return 0;
}