精度计算-大数阶乘

精度计算-大数阶乘

本算法的目的在于计算一个比较大的数的阶乘，由于得到的结果比较大，是现有的数据类型无法存储的，所以我决定将结果存储在一个long a[]数组中。 我们的思路是把每4位数看做数组的一个元素来存储，例如：个、十、百、千存在a[0]，万、十万、百万、千万存在a[1]以此类推。 我们用10的阶乘来模拟一下求结果大于4位数阶乘的过程，9的阶乘为362880，而10的阶乘为9的阶乘乘以10，在计算完9的阶乘时a[0] = 2880,a[1]=36,因为362880*10 = (36*10+2880*10/10000)*10000+28800%10000 = 3628800。同理其他高阶的数的计算过程和这个大同小异了 下面是我的C语言实现过程：

void factorial(int n)
{
long a[10000];
int i,j,c,m=0; 
a[0]=1; 
for(i=1;i<=n;i++)
{ 
c=0; 
for(j=0;j<=m;j++)
{
a[j]=a[j]*i+c;
c=a[j]/10000;
a[j]=a[j]%10000;
}
printf("\n");
if(c>0) 
{
m++;
a[m]=c;
}
}
printf("\n%ld",a[m]);
for(i=m-1;i>=0;i--)
printf("%4.4ld",a[i]);
}