nyoj___大數階乘

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

大數階乘

描述

我們都知道如何計算一個數的階乘，可是，如果這個數很大呢，我們該如何去計算它並輸出它？

輸入

輸入一個整數m(0<m<=5000)

輸出

輸出m的階乘，並在輸出結束之後輸入一個換行符

樣例輸入

50

樣例輸出

30414093201713378043612608166064768844377641568960512000000000000

程序代碼:

     

 1 #include<stdio.h>
 2 #include<math.h>
 3 int main(){
 4     int n;
 5     while(scanf("%d",&n)!=EOF){
 6             long a[10000];
 7             int i,j,l,c,m=0,w;
 8             a[0]=1;
 9             for(i=1;i<=n;i++){
10                 c=0;
11                 for(j=0;j<=m;j++){
12                     a[j]=a[j]*i+c;
13                     c=a[j]/10000;
14                     a[j]=a[j]%10000;
15                 }
16                 if(c>0) {m++;a[m]=c;}
17             }
18 
19         w=m*4+log10(a[m])+1;
20         printf("%ld",a[m]);
21         for(i=m-1;i>=0;i--) printf("%4.4ld",a[i]);
22         printf("\n");
23     }
24     return 0;
25 }

大數階乘位數的計算方法是(int)log10(x)+1 
這個是計算x的位數的
log10(n!)=log10(1)+log10(2)+log10(3)+...+log10(n)