求N的阶乘的位数（大整数）

2011-04-11 20:48:22

朴素公式为：

    log10(n!)=log10(1*2*3…*n)=log10(1)+log10(2)+…+log10(n)

求位数：log10(n!)=log10(1*2*3…*n)=log10(1)+log10(2)+…+log10(n)，对log10(n!)的值取整加1就是n!的位数。

《计算机程序设计艺术》中给出了另一个公式

    n! = sqrt(2*π*n) * ((n/e)^n) * (1 + 1/(12*n) + 1/(288*n*n) + O(1/n^3))

    π = acos(-1)

    e = exp(1)

两边对10取对数

忽略log10(1 + 1/(12*n) + 1/(288*n*n) + O(1/n^3)) ≈ log10(1) = 0

得到公式

    log10(n!) = log10(sqrt(2 * pi * n)) + n * log10(n / e)

FOJ Problem 1050 Number lengths

题目地址：http://acm.fzu.edu.cn/problem.php?pid=1050

---

#include<iostream> #include<cmath> using namespace std; int main() { int n; double Pi,E; Pi=acos(-1.0); E=exp(1.0); while(~scanf("%d",&n)) { int len; len=(int)(log10l(sqrtl(2*Pi*n))+n*log10l(n/E)+1); // n>1是该公式成立！！！！！！！！！！！ if(n==1) len=1; cout<<len<<endl; } return 0; }