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;
}