poj1001 Exponentiation --高精度

题意: 

        题意很简单，就是让求一个数R的N次方。

        其中，有多组输入，对应多组输出，数据范围（ 0.0 < R < 99.999 ，  0 < n <= 25 ）。

思路：

       （1） 从数据规模来看，求R的N次方，R为浮点型，long long 最大可表示20位左右，没有任何疑问，要用高精度来求（仔细想想，为什么？）。

       （2） 本题要用高精度的话，就要处理小数点，可以先把R转化为整型，最后再转化回来。

代码:

         本来自己写了个代码本打算发上来，不过感觉写的实在太差了(刚开始各种WA，后来在Discuss搜数据才A。。。) ，先发一个网上找的感觉比较清晰的。。。

//Memory Time 
//1232K  0MS 


#include<iostream>
#include<string>
using namespace std;


const int size=1000;  //大数位数


void mult(char* A,char* B,char* ans)
{
	int i,j,k;


	int fract;   //总小数个数
	int dot=-1;  //小数点位置
	for(k=0;A[k]!='\0';k++)
		if(A[k]=='.')
			dot=k;
	int lena=k;//数据的长度
	if(dot==-1)
		fract=0;
	else
		fract=lena-dot-1;//总小数个数


	dot=-1;
	for(k=0;B[k]!='\0';k++)
		if(B[k]=='.')
			dot=k;
	int lenb=k;
	if(dot==-1)
		fract+=0;
	else
		fract+=(lenb-dot-1);  //总小数个数


	int a[size+1]={0};
	int b[size+1]={0};
	int pa=0,pb=0;


	/*倒序*/


	for(i=lena-1;i>=0;i--)
	{
		if(A[i]=='.')
			continue;
		a[pa++]=A[i]-'0';//转移到int型数组 a [  ]  里
	}
	for(j=lenb-1;j>=0;j--)
	{
		if(B[j]=='.')     //暂时删除小数点
			continue;
		b[pb++]=B[j]-'0';//转移到int 型数组b [  ]  里
	}


	int c[2*size+1]={0};
	int lenc;
	for(pb=0;pb<lenb;pb++)//lenb >=lena ！！     for    for  高精度核心算法
	{
		int w=0;  //低位到高位的进位
		for(pa=0;pa<=lena;pa++)  //       = 为了处理最后的进位
		{
			int temp=a[pa]*b[pb]+w;
			w=temp/10;
			temp=(c[pa+pb]+=temp%10);
			c[lenc=pa+pb]=temp%10;//注意： lenc 
			w+=temp/10;
		}
	}


	/*倒序，得到没有小数点的ans*/


	for(pa=0,pb=lenc;pb>=0;pb--)
		ans[pa++]=c[pb]+'0';//   char类型   ans [  ]
	  
	   ans[pa]='\0';
	    lena=pa;


	 /*插入小数点*/


	bool flag=true; //标记是否需要删除小数末尾的0
	if(fract==0)   //小数位数为0，无需插入小数点
		flag=false;
	else if(fract<lena) //小数位数小于ans长度，在ans内部插入小数点
	{
		ans[lena+1]='\0';
		for(i=0,pa=lena;pa>0;pa--,i++)
		{
			if(i==fract)
			{
				ans[pa]='.';
				break;
			}
			else
				ans[pa]=ans[pa-1];
		}
				
	}
	else //小数位数大于等于ans长度，在ans前面恰当位置插入小数点
	{
		char temp[size+1];
		strcpy(temp,ans);
		ans[0]='0';
		ans[1]='.';
		for(int i=0;i<fract-lena;i++)  //补充0
			ans[i+2]='0';
		for(j=i,pa=0;pa<lena;pa++)
			ans[j++]=temp[pa];
		ans[j]='\0';
	}


	/*删除  ans 小数  末尾  的0*/


	if(flag)
	{
		lena=strlen(ans);
		pa=lena-1;
		while(ans[pa]=='0')
			ans[pa--]='\0';
		if(ans[pa]=='.')   //小数全为0
			ans[pa--]='\0';
	}


	/*删除ans整数开头的0，但至少保留1个0*/


	pa=0;
	while(ans[pa]=='0')  //寻找ans开头第一个不为0的位置
		pa++;


	if(ans[pa]=='\0')  //没有小数
	{
		ans[0]='0';
		ans[1]='\0';
	}
	else  //有小数
	{
		for(i=0;ans[pa]!='\0';i++,pa++)
			ans[i]=ans[pa];
		ans[i]='\0';
	}
	return;
}


char a[size+1];
char ans[size*size+1];//百万


int main(  )
{
	int b;
	while(cin>>a>>b)
	{
		memset(ans,'\0',sizeof(ans));
		ans[0]='1';//初始化
		ans[3]='\0';


		for(int i=1;i<=b;i++)//b次
			mult(a,ans,ans);


		cout<<ans<<endl;
	}
	return 0;
}

    说实话，示例程序写的不算很简洁高效，但能0ms过，着实令人惊奇。思路也很简单，就是把R的n次方 一次次模拟为两个数相乘。 这也提醒我们，只要满足时间和空间要求，算法越简单越好。                         切记！！！！