计算大数的阶乘

问题：

计算整数N的阶乘，要求计算的结果最长可以到50位！（50位十进制而不是二进制）。

 

先上代码，然后分析：

#include<stdio.h>
#define MAXLEN 50

int data[MAXLEN];

int *calcfact(int N)
{
	int i,j,k,r;
	int digit = 1;

	for(i = 1; i < MAXLEN + 1; i++)
		data[i] = 0;

	data[0] = 1;
	data[1] = 1;

	for(i = 1; i < N + 1; i++)
	{
		for(j = 1; j < digit + 1; j++)
			data[j] *= i;

		for(j = 1; j < digit + 1; j++)
		{
			if(data[j] > 10)
			{
				for(r = 1; r < digit + 1; r++)
				{
					if(data[digit] > 10)
						digit++;
					data[r + 1] += data[r] / 10;
					data[r] = data[r] %10;
				}
			}
		}

		printf("%d! = ",i);

		for(k = digit; k > 0; k--)
			printf("%d",data[k]);
		printf("\n");
		
	}
	return data;
}

int main()
{    
	int N;

	printf("Enter a number : ");
	scanf("%d",&N);
	*calcfact(N);
	
	getchar();
	getchar();
}

解析：主要是通过数组的思想来做。因为阶乘太大，一般的数据类型是放不下的！所以可以用数组来做，数组每个元素都在0~9之间，数组的一个元素表示一位！比如：

data[9] = 2;

data[8] = 3;

data[7] = 4;

data[6] = 6;

data[5] = 8;

data[4] = 3;

data[3] = 1;

data[2] = 3;

data[1] = 0;

data[0] = 2;

那么这个数就可以表示为：2346831302~~~

我们只需要把这个数组按顺序printf出来即可，不需要再进行合成。只是在计算的时候，如果这个数组里面的元素大于9，就应该将其进位！

 

总结：

1. 一般在32位机上，double类型占64位，一般能表示20位十进制。

2. 在最开始编程的时候，将return的位置写错了，囧，所以计算完1的阶乘就退出了。

3. 编写的程序都是参考《经典C面试 真题精讲》里面的，虽然里面有很多文字上的错误，不过还是一本很好的书。

4. 哪位大虾能推荐些更好的书，可以留言哈！~感谢~