POJ 1001 求高精度幂（高精度）

Description

对数值很大、精度很高的数进行高精度计算是一类十分常见的问题。比如，对国债进行计算就是属于这类问题。现在要你解决的问题是：对一个实数R，要求写程序精确计算R的n次方R^n，其中n 是整数并且0.0 < R < 99.999,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

Solution

高精度乘法，注意前置零和后置零以及整数时去掉小数点，其中有一个显然的结论，小数点原来位置是pos，那么结果小数点的位置应为pos*n

Code

#include<stdio.h>
#include<string.h>
int a[1111111],b[1111111],temp[1111111],la,lb;
void multi()
{
memset(temp,0,sizeof(temp));
for(int i=0;i<la;i++)
for(int j=0;j<lb;j++)
{
temp[i+j]+=a[i]*b[j];
if(temp[i+j]>9)
{
temp[i+j+1]+=temp[i+j]/10;
temp[i+j]%=10;
}
}
if(temp[la+lb-2]>9)//进位
{
temp[la+lb-1]+=temp[la+lb-2]/10;
temp[la+lb-2]%=10;
}
la=la+lb;
for(int i=0;i<la;i++)
a[i]=temp[i];
}
int main()
{
char s[10];
int n;
while(~scanf("%s%d",s,&n))
{
int pos=-1;//小数点位置
for(int i=5,j=0;i>=0;i--)
if(s[i]!='.')
a[j]=s[i]-'0',b[j]=a[j],j++;
else
pos=i;
la=lb=(pos==-1?6:5);
for(int i=1;i<n;i++)
multi();
if(pos==-1)//没有小数点直接输出结果即可
{
for(int i=la-1;i>=0;i--)
printf("%d",a[i]);
printf("\n");
}
else//有小数点
{
int l,r;
pos=5-pos;
pos*=n;//结果中小数点位置
for(int i=0;i<la;i++)//后置零
if(a[i]!=0)
{
l=i;
break;
}
for(int i=la-1;i>=0;i--)//前置零
if(a[i]!=0)
{
r=i;
break;
}
if(r<pos)r=pos-1;
if(l>pos)l=pos;
for(int i=r;i>=l;i--)
{
if(i==pos-1)printf(".");
printf("%d",a[i]);
}
printf("\n");
}
}
return 0;
}