大数阶乘

大数阶乘
返回：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