超出基本数据类型范围的n!问题

#include <iostream>
#include <string.h> 
using namespace std;
#define len 3000
int a[len];
//将得到的数字都保存在数组中，i由0~len表示位数由低到高，例如个、十、百、千、万 

int main()
{ 
	int n,j;
	cin>>n;
	memset(a,0,sizeof(a));//对数组进行清0 
	a[0]=1;//最低位为1，作为n！中的1进行运算 
	for(int i=2;i<=n;i++)
	{//i作为阶乘中的n进行运算 
	//模拟乘法 
		int c=0;//作为对高位的进位 
		for(j=0;j<len;j++)
		{
			int s=a[j]*i+c;//已有的阶乘i！的每一位乘以i+1，并加上进位的数 
			a[j]=s%10;//如果s>10,向高位进位，只把小于10的部分保留到本位 
			c=s/10; //将大于10的部分进位到高位，假设i=32，而a[j]=4,c=0，则应向高位进位12，而12会在下次循环中在高位再次进行处理 
		} 
	}
	//因为开始时对数组进行了置0，而高位在后，因此需要从后往前输出，还有跳过仍为0的位 
	j=len;
	while(!a[--j]);//到达第一个非0处停止，即为结果的最高位
	while(j>=0) cout<<a[j--]; 
} 

PS:在该程序中我个人定义的范围是在1000以内的数的阶乘，1000！计算之后可大概用长度为3000的int数组进行保存，若需要更大的数进行计算的话，更改宏定义即可