斯特林公式 ——Stirling公式(取N阶乘近似值)

   斯特灵公式是一条用来取n阶乘近似值数学公式。一般来说,当n很大的时候,n阶乘的计算量十分大,所以斯特灵公式十分好用。从图中可以看出,即使在n很小的时候,斯特灵公式的取值已经十分准确。

                                                    

公式为   n! \approx \sqrt{2\pi n}\, \left(\frac{n}{e}\right)^{n}.

 从图中看出,对于足够大的整数n,这两个数互为近似值。更加精确地:   

      \lim_{n \rightarrow \infty} {\frac{e^n\, n!}{n^n \sqrt{n}}} = \sqrt{2 \pi}.       或者        \lim_{n \rightarrow \infty} {\frac{n!}{\sqrt{2\pi n}\, \left(\frac{n}{e}\right)^{n}}} = 1

 这个公式,以及误差的估计,可以推导如下。我们不直接估计n!,而是考虑它的自然对数

   \ln(n!) = \ln 1 + \ln 2 + \cdots + \ln n.

按一般方法计算N的阶乘,其时间复杂度为O(N):    N!= 1 * 2 * 3 * 4 * 5 * ............ * N;


如果要计算N后得到的数字为几位数,则我们可以知道其位数等于lgN!+1;

则: 
\ln(n!) = \ln 1 + \ln 2 + \cdots + \ln n.

但是当N很大的时候,我们可以通过斯特林公式进行优化:(即Stirling公式)

n! \approx \sqrt{2\pi n}\, \left(\frac{n}{e}\right)^{n}.(e = 2.718)

斯特林公式可以用来估算某数的大小,结合lg可以估算某数的位数,或者可以估算某数的阶乘是另一个数的倍数。

例题:  http://acm.hdu.edu.cn/showproblem.php?pid=1018
题目给出的N的范围是: 1<= N <= 107  
用普通方法肯定算不出N的阶乘后的出的数字位数,但运用斯特林公式则很好解决

#include < stdio.h >
#include < iostream >
#include < math.h >
using namespace std;

int main ()
{
    int w;
    scanf("%d",&w);
    while(w--)
    {
        int z, i;
        double sum = 0;
        scanf("%d",&z);
        for(i = 1; i <= z; i++)
            sum += log10( double(i) );
        printf("%d\n",int(sum) + 1);
    }
    return 0;
}



  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值