求10的7次方阶乘的位数

即使找到最先进的计算大数阶乘的算法， 
若完全计算10的7次阶乘，然后统计其位数，也需要几十分钟才行。 
所以不可行，况且题目只需要求其位数。

那该怎么计算呢？非常的简单，其实是纯粹的数学问题。 
现由斯特林[striling]公式可得： 
lnN!=NlnN－N+0.5ln(2N*pi)

而10的7次方阶乘的位数等于： 
log10(N!)取整后加1 
log10(N!)=lnN!/ln(10)

所以整个结果用c描述：

#include<stdio.h> 
#include<math.h> 
#define PI 3.14159265 
#define N 10000000 
int main(){ 
    int len; 

    len=ceil((N*log(N)-N+log(2*N*PI)/2)/log(10));
    //ceil求上界，即不小于某值的最小整数 
    //log在c语言里面是自然对数 
    printf("%d\n",len); 
    return 0;
} 

 

len=65657060