nyoj___大数阶乘

http://acm.nyist.net/JudgeOnline/problem.php?pid=28

大数阶乘

时间限制：3000 ms  |  内存限制：65535 KB

难度：3

描述

我们都知道如何计算一个数的阶乘，可是，如果这个数很大呢，我们该如何去计算它并输出它？

输入

输入一个整数m(0<m<=5000)

输出

输出m的阶乘，并在输出结束之后输入一个换行符

样例输入

50

样例输出

30414093201713378043612608166064768844377641568960512000000000000

程序代码:

     

#include<stdio.h>
#include<math.h>
int main(){
    int n;
    while(scanf("%d",&n)!=EOF){
            long a[10000];
            int i,j,l,c,m=0,w;
            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;}
            }

        w=m*4+log10(a[m])+1;
        printf("%ld",a[m]);
        for(i=m-1;i>=0;i--) printf("%4.4ld",a[i]);
        printf("
");
    }
    return 0;
}

大数阶乘位数的计算方法是(int)log10(x)+1 
这个是计算x的位数的
log10(n!)=log10(1)+log10(2)+log10(3)+...+log10(n)