用Stirling逼近近似计算阶乘

用Stirling逼近近似计算阶乘n!

1、求n！的前几位

先看公式：

θ十分接近1，而且在逐渐地逼近1，实际上，即使是求1的阶乘，θ也会达到0.9727376027，这是一个本身就是一个很“精确”的数字了！当n＝1000时，θ将0.99999996665875876427498746773752，与1的差别只有0.000000033341241235725012532263（约等于3.33412×10^-8）！

 

这样可以推出求阶乘前几位的最后公式为：

 

2、求位数

lg(n)+1=n的位数  （注意：lg是指以10为底的对数）

继而推出

 

--------------------------------------------------------------------------------------------------下面附上用C++写的代码：

 

#include<iostream>

#include<cmath>

using namespace std;

int fac(int n){//递归求阶乘（很有限）

         if(n>1){

           return n*fac(n-1);

         }else{

           return n;

         }

}

 

double frac(double n){

   return n-floor(n);

}

 

int stirlingp(int n){//计算阶乘结果的位数

   return floor(0.5*log10(2*n*3.1415926)+n*log10(n/exp(1.0)))+1;

}

 

double stirling(int n,int p){//计算阶乘n！结果的前面（从左到右）p位

  return pow(10.0,frac(0.5*log10(2*n*3.1415926)+n*log10(n/exp(1.0))))*exp(1.0/12/n)*pow(10.0,p-1);

}

 

 

 

void main()

{

         cout<<stirling(4,1)<<"E+"<<stirlingp(4)<<endl;

         getchar();

}

--------------------------------------------------------------------------------------------------

文本公式推导源于《用Stirling逼近近似计算阶乘的探讨与应用》仲晨 2005年写的一篇文章那时还是高中生，实在厉害，文本只做结果使用！