阶乘计算

题目
输入一个正整数n，输出n!的值。（n<=1000）

描述
　　其中n!=123*…*n。
　　n!可能很大，而计算机能表示的整数范围有限，需要使用高精度计算的方法。使用一个数组A来表示一个大整数a，A[0]表示a的个位，A[1]表示a的十位，依次类推。
　　将a乘以一个整数k变为将数组A的每一个元素都乘以k，请注意处理相应的进位。
　　首先将a设为1，然后乘2，乘3，当乘到n时，即得到了n!的值。

注意：
此题与高精度加法类似，都要考虑进位的问题，但是本题中对于最终的结果a，需要通过计算才能知道长度，因此采取将数组A全初始化为0的方法，对A[i]中的每个数进行计算并进位，结束后，用一个for循环来判断A[i]中的非0数（结果a的最高位），循环输出结果a。

代码：

#include <iostream>
#include <string.h>
using namespace std;
#define M 100000//M值设大一些，否则对于n大的数而言，得到的结果储存不下 
int a[10];
int main()
{
	int n;
	cin >> n;
	int a[M],r = 0;   //进位r	
	memset(a,0,sizeof(a));//若把a[M]在main()函数外定义，使其成为全局变量，则无需此行来清空数组使值为0
	a[0] = 1;
	int i, j;
	for ( i = 2;i <= n;i++)
	{
		for ( j = 0;j < M;j++)//每一次从a[0]开始循环至a[M-1]，若进位r=0，则十位数及更高位数皆为0； 
		{
			a[j] = a[j] * i + r;//i=2:a[0]=2,r=0,a[1]~a[M-1]=0;
			r = a[j] / 10;//i=3:a[0]=6,r=0;a[1]~a[M-1]=0;
			a[j] = a[j] % 10;//i=4:a[0]=4,r=2,a[1]=2,a[2]~a[M-1]=0;......直至循环到n 
		}
	}
 
	for (i = M - 1;i > 0;i--)//从a[M-1]开始循环，因为未给a数组赋值前，数组内皆为0
		if (a[i])//因此，a[i]为真，即a[i]非0时，跳出这个for循环 
			break;
	for (j = i;j >= 0;j--)//跳出后将a[i]！=0中的i值赋给j，此时的j为所求结果的最高位，从j循环输出a数组值 
		cout << a[j] ;
	cout<<endl;
}