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

Source

East Central North America 1988

Translator

北京大学程序设计实习，Xie Di

思路：首先以字符串接收，然后挨个输出到数组内，大数相乘，最后得到小数点的位置，然后输出。

#include<stdio.h>
#include<string.h>
int a[1111111],b[1111111],temp[1111111];
int la,lb;
void multiply() {
	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]=temp[i+j]%10;
			}
		}
	if(temp[la+lb-2]>9) {
		temp[la+lb-1]+=temp[la+lb-2]/10;
		temp[la+lb-2]=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 position=-1,j=0;
		for(int i=5; i>=0; i--) {
			if(s[i]!='.') {
				a[j]=s[i]-'0';
				b[j]=a[j];
				j++;

			} else {//获得小数点的位置
				position=i;
			}
		}
		la=lb=(position==-1?6:5);
		for(int i=1; i<n; i++)
			multiply();//进行乘法运算
		if(position==-1) { //没有小数点则直接输出
			for(int i=la-1; i>=0; i--) {
				printf("%d",a[i]);
			}
			printf("\n");
		} else {
			int l,r;
			position=5-position;
			position*=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<position) r=position-1;
			if(l>position) l=position;
			for(int i=r; i>=l; i--) {
				if(i==position-1) {
					printf(".");
				}
				printf("%d",a[i]);

			}
			printf("\n");
		}
	}
	return 0;
}