poj1001高精度

链接:http://poj.org/problem?id=1001

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

 

这道高精度计算题需要应用字符串，因为常规存储类型无法存储如此巨大的数据，只能将数字以字符串的方式存储与处理。

由于字符串无法直接相乘，需要将字符串数组中的每一位字符转化成整型（与’0’的差），再应用乘法法则（笔算乘法的方法）进行计算。

 

本题中需要注意的点：

1.先找到小数点的位置，纪录并去掉。

2.去掉无意义的0。

3.纪录小数点后有效位的位数。

4.每次计算中都要去除无意义的0。

5.计算最后结果中小数点的位置。

6.将小数点加入到计算结果中。

7.去掉最终结果中的无意义的0并输出。

 

 

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int pointpos;
int output[176],num[6],n1,ll;
void multiply()
{
	int i;
	for(i=0;i<ll;i++){
		output[i]*=n1;
	}
	    i=0;
	while(output[i]>9||i<ll){
		output[i+1]+=output[i]/10;
		output[i]=output[i]%10;
		i++;
	}
	ll=i+1;
}
int main()
{
	char c[8];
	int n,i,j,len,temp,first,last,flag;
	while(scanf("%s%d",c,&n)!=EOF){
	   pointpos=0;
	   temp=len=strlen(c);
	   flag=0;
	   for(i=0;i<len;i++){
   		if(c[i]=='.') flag=1;
   	}
	   for(i=0;i<len;i++){
        	if(c[i]!='0'){
               first=i;
               break;
	        }
        }
        if(flag){
  	   for(i=len-1;i>=first;i--){
  	   	if(c[i]=='.'){
  	   		last=i-1;
	   	  	break;
	   	  }
     	if(c[i]!='0'){
        last=i;
        break;
	    }
   	  }
        }else last=len-1;
   	  temp=0;
     	memset(output,0,sizeof(output));
  	   for(i=last;i>=first;i--){
           if(c[i]!='.'){
           	output[temp]=num[temp]=c[i]-'0';
           	temp++;
           } else{
           	pointpos=last-i;
           }
	    }
	    if(n<2){
	    	for(i=first;i<=last;i++)
    		printf("%c",c[i]);
    		printf("\n");
    		continue;
    	}
    	n1=0;
    	for(i=temp-1;i>=0;i--){
	    	n1=n1*10+num[i];
	    }
	    if(n1==0){
    	printf("0\n");
		continue;	
    	}
	    ll=temp;
	    for(i=1;i<n;i++){
    		multiply();
    	}
    	for(i=0;i<177;i++){
    		if(i==pointpos) {last=pointpos;break;}
	    	if(output[i]!=0){
	    	last=i;
			break;	
	    	}
	    }
    	for(i=175;i>=0;i--){
	    	if(output[i]!=0){
	    	first=i;
			break;	
	    	}
	    	if(i==pointpos) {first=pointpos-1;break;}
	    }
	    pointpos*=n;
	    if(pointpos-1<last){
    		for(i=first;i>=last;i--)
    		printf("%d",output[i]);
    	}else if(pointpos>first){
    		printf(".");
	    	for(i=0;i<pointpos-first-1;i++)
	    	printf("0");
	    	for(i=first;i>=last;i--)
    		printf("%d",output[i]);
	    }else{
    		for(i=first;i>=last;i--){
    			printf("%d",output[i]);
		    	if(pointpos==i)printf(".");
		    }
    	   }
    	   printf("\n");
 }		
	return 0;
}