大数阶乘

大数阶乘
返回：n!的位数

结果：a[10000]记录n!的数据

#include<stdio.h>
#include<math.h>
int factorial(int n);
int main()
{
int re;
re=factorial(1000); //测试用的
printf("\n%d\n",re);
return 0;
}
int factorial(int n)
{
long a[10000];//int类型也可以
int i,j,l,c,w,m=0;
a[0]=1;
for(i=1;i<=n;i++)//i小于等于n，注意边界问题
{
c=0;
for(j=0;j<=m;j++)
{
a[j]=a[j]*i+c;
c=a[j]/10000;//取4位以上的数据
a[j]%=10000;
}
if(c>0)
{
m++;
a[m]=c;
}
}
w=m*4+log10(a[m])+1;//math里面的函数
printf("\n%ld",a[m]);
for(i=m-1;i>=0;i--)
printf("%4.4ld",a[i]);//如果a[i]不够4位 在前面自动补零
return w;
}

by  吴尚奇   Devil_box   2014/06