POJ 1001 Exponentiation

POJ 1001 题目描述:http://poj.org/problem?id=1001

 题目整体是一个求实数的n次幂的问题。可能会遇到2个问题，第一个是不确定行数的输入问题，这个问题在题目的hint中有所提示，所以应该不是个大问题，没能实现的可以看我下面的代码，自己仔细揣摩一下。第二个就是求幂的问题，这里直接使用数组实现的大整数乘法。

大整数乘法

先来使用列表法研究一下，可以说也是一种找规律吧：NOTICE 图片上最终结果算错了，应该是2051010

从上面的图片中找到规律后，给出相应的伪代码：

x=[]
y=[]
ans=[0]
for i from 1 to length_of_x:
    for j from 0 to length_of_y:
        ans[i+j]+=x[i]*y[j];
for i from length_of_ans to 1:
    if(ans[i]>=10)
       ans[i-1]+=ans[i]/10
       ans[i]%=10

c++的简单实现版：实现了上图中的例子

#include <iostream>

using namespace std;

int main()
{
    int x[]={0,8,7,6,5};
    int y[]={2,3,4};
    int ans[10]={0};
    
    for(int i=1;i<=4;i++)
        for(int j=0;j<=2;j++)
            ans[i+j]+=x[i]*y[j];
    for(int i=6;i>=1;i--)
        if(ans[i]>=10)
        {
            ans[i-1]+=ans[i]/10;
            ans[i]%=10;
        }
    for(int i=0;i<=7;i++)
        cout<<ans[i];
        return 0;
}

虽然已经实现了大整数的乘法，而且可以很方便的应用于实数乘法，但还有一些琐碎的问题需要去解决。

1输入问题

这里指的输入问题并不是hint里提示的，而是将实数作为字符串输入，然后进行运算前将其处理的问题。

先将上面的大整数乘法中x,y,ans的类型改成char型

char x[]={'0','8','7','6','5'};
    char y[]={'2','3','4'};
    char ans[10]={'0'};
    for(int i=1;i<=4;i++)
        for(int j=0;j<=2;j++)
            ans[i+j]+=(x[i]-'0')*(y[j]-'0');
    for(int i=6;i>=1;i--)
        if(ans[i]>=10)
        {
            ans[i-1]+=ans[i]/10;
            ans[i]%=10;
        }
    if(ans[0]-48!=0)
        cout<<ans[0]-48;
    for(int i=1;i<=6;i++)
        cout<<int(ans[i]);

经过这样的改写，就可以处理以char类型输入的数据了

2实数(主要是含小数点的数据处理)处理

我们通过先把小数点去除转化为整数除法，然后输出时加上小数点来解决

通过代码

#include <stdio.h>
#include <string.h>
char R[7];
int n,dot;
char a[150];
char ans[150];
void mul(char *b)
{
	int i,j,f,t;
	for(i=0;ans[i];i++)
	{
		f=0;
		for(j=0;b[j];j++)
		{
			t=((ans[i]-'0')*(b[j]-'0')+(a[i+j]-'0')+f)%10;
			f=((ans[i]-'0')*(b[j]-'0')+(a[i+j]-'0')+f)/10;
			a[i+j]=t+'0';
		}
		if(f) a[i+j]=f+'0';
	}
	//m位乘以n位结果为m+n或m+n-1位
	//如果为m+n-1位，此处多放了一个0，但最后print时会过滤掉
	a[i+j]='/0';
	strcpy(ans,a);
	memset(a,'0',sizeof(a));
}
void print()
{
	int front,back;
	front=strlen(ans)-1;
	back=0;
	while(ans[back]=='0') back++;
	while(front>=dot&&ans[front]=='0') front--;
	while(front>dot) printf("%c",ans[front--]);
	if(back<dot)
		printf(".");
	while(dot>=back) printf("%c",ans[dot--]);
	printf("/n");
}
int main()
{
	char *t;
	while(scanf("%s%d",R,&n)!=EOF)
	{
		memset(a,'0',sizeof(a));
		t=strchr(R,'.');
		dot=n*(strlen(R)-(t-R)-1)-1;
		while(*t) *t++=*(t+1);
		strrev(R);//用C提交Compile Error
		strcpy(ans,R);
		n--;//第一次调用mul时就算了两个数的乘积
		while(n--) mul(R);
		print();
	}
	return 0;
}