算法分析与设计,大整数的存储与运算,运算结果超过数据表示范围
例如:100!一百的阶乘,结果为 100!=93 326215 443944 152681 699238 856266 700490 715968 264381 621468 592963 895217 599993 229915 608941 463976 156518 286253 697920 827223 758251 185210 916864 00000 00000 00000 00000 结果已经远远超过longlong的表示范围
方法为:用一个长整型数组存储运算结果,每个元素存储6位数字,以减少累乘次数和存储空间的消耗。
源码:
#include<stdio.h>
void main(){
long a[256],b=0,d=0;
int m=0,n=0,i=0,j=0,r=0;
//循环变量i表示要累乘的数据,循环变量j表示当前累乘结果的数组下标,数据b存储计算的中间结果,d为存储超出6位后的进位
scanf("%d",&n);m=log(n)*n/6+2;//估计占据多少位数组空间
a[1]=1;
for(i=2;i<=m;i++)
a[i]=0;
for(i=2;i<=n;i++){
d=0;
for(j=1;j<=m;j++){
b=a[j]*i+d;
a[j]=b%1000000;
d=b/1000000;
}
if(d<0||d>0)
a[j]=d;
}
for(i=m;i>=1;i--){
if(a[i]==0)
continue;
else{
r=i;
break;
}
}//此循环计算结果的准确位数
//输出结果
printf("%d!=",n);
printf("%ld",a[r]);
printf(" ");
for(i=r-1;i>=1;i--)
{
if(a[i]>99999)
{
printf("%ld",a[i]);
printf(" ");
continue;
}
if(a[i]>9999){
printf("0",a[i]);
printf(" ");
continue;
}
if(a[i]>999){
printf("00",a[i]);
printf(" ");
continue;
}
if(a[i]>99){
printf("000",a[i]);
printf(" ");
continue;
}
if(a[i]>9){
printf("0000",a[i]);
printf(" ");
continue;
}
printf("00000",a[i]);
printf(" ");
}
}
运行截图: