数学问题--大数的阶乘

#include <stdio.h>
#include <math.h>

int 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;
}
if(c>0) {m++;a[m]=c;}
}

printf("%ld",a[m]);
for(i=m-1;i>=0;i--) printf("%4.4ld",a[i]);

int w=m*4+(int)log10(a[m])+1;
return w;
}

void main()
{
factorial(20);
}