高精度运算 程序崩溃

gdczsue 2013-02-03 06:12:17
求高精度幂
Time Limit: 500MS Memory Limit: 10000K
Total Submissions: 113086 Accepted: 27443
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
#include<iostream>
#include<cstring>
#define maxn 250
using namespace std;
char c[maxn];//存取结果
char temp[maxn];//临时存取结果
int len;
int point(char* oc,char* nc)//去掉小数点
{
int j=1;
for(int i=5;i>=0&&oc[i]!='.';i--)
nc[j++]=oc[i];
int k=(i==-1?6:i+1);//可能没有小数点
i--;
while(i>=0)
nc[j++]=oc[i--];
return 6-k;//小数位数
}
void multiply(char* a,char* b)//字符数组相乘
{
memset(c,'0',sizeof(c));
int la=strlen(a)-1;
int lb=strlen(b)-1;
len=la+lb;//结果不会超出len
for(int i=1;i<=la;i++)
for(int j=1;j<=lb;j++)
c[i+j-1]=c[i+j-1]+(a[i]-'0')*(b[j]-'0');
for(i=1;i<len;i++)
{
c[i+1]=c[i+1]+(c[i]-'0')/10;
c[i]=(c[i]-'0')%10+'0';
}
if(c[len]=='0')
len--;
}
int main()
{
char s[7];
int n;
while(cin>>s>>n)
{
char a[7];
int p=point(s,a);//记录小数位数
strcpy(c,a);//把底数赋给c
for(int i=1;i<n;i++)
{
strcpy(temp,c);//用temp暂时存储结果
multiply(temp,a);
}
for(i=len;i>p*n;i--)//输出小数点前面的数
cout<<c[i];
cout<<'.';
for(i=p*n-1;i>0;i--)//输出小数点后面的数
cout<<c[i];
cout<<endl;
}
return 0;
}


...全文
133 2 点赞 打赏 收藏 举报
写回复
2 条回复
切换为时间正序
当前发帖距今超过3年,不再开放新的回复
发表回复
gdczsue 2013-02-03
引用 1 楼 healer_kx 的回复:
Debug哪里出问题啊?
debug 中间程序终止了
  • 打赏
  • 举报
回复
healer_kx 2013-02-03
Debug哪里出问题啊?
  • 打赏
  • 举报
回复
相关推荐
发帖
C++ 语言
加入

6.0w+

社区成员

C++ 语言相关问题讨论,技术干货分享,前沿动态等
申请成为版主
帖子事件
创建了帖子
2013-02-03 06:12
社区公告
暂无公告