精度计算-大数阶乘

精度计算-大数阶乘

本算法的目的在于计算一个比较大的数的阶乘，由于得到的结果比较大，是现有的数据类型无法存储的，所以我决定将结果存储在一个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]);
}