Wallis公式&Stirling公式&Gamma函数

.DoubleBean.

已于 2022-05-27 22:00:06 修改

阅读量3.1k

点赞数 2

分类专栏：数学文章标签：算法几何学

于 2021-03-07 09:27:51 首次发布

本文链接：https://blog.csdn.net/weixin_45771864/article/details/114435753

版权

数学专栏收录该内容

5 篇文章 0 订阅

订阅专栏

本文探讨了如何使用Wallis公式来推导斯特林公式，斯特林公式在大数值阶乘计算中起到优化作用。通过分析lnx函数的图像面积，文章阐述了积分法、内接梯形分割法和外切梯形分割法，并指出积分法最为精确。斯特林公式简化了n阶乘的计算，尤其在n较大时，能有效估算n!的位数。此外，还介绍了利用斯特林公式估算数字位数的方法以及与nlogn的等价无穷大关系。最后，给出了C++代码示例来计算n!的位数。

摘要由CSDN通过智能技术生成

Wallis公式

(!!是双阶乘，不超过这个整数且具有相同奇偶的正整数相乘，比如5!!=1×3×5, 6!!=2×4×6)
在这里插入图片描述

在这里插入图片描述

利用Wallis公式推导Stirling公式

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

斯特林公式如下
在这里插入图片描述
接下来利用Wallis公式来推导斯特林公式

借助函数 f(x) = ln x 的图像面积，通常有三种求法，分别是积分法，内接梯形分割法，外切梯形分割法。实际上最准确的是第一种，后面两种都有一定的误差

在这里插入图片描述

如果要计算n的阶乘后得到的数字为几位数，则我们可以知道其位数等于lg n! + 1

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

但是当n很大的时候，我们可以通过斯特林公式进行优化：（即Stirling公式）

（e = 2.718）

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

再来证明log(n!) 与 nlog n是等价无穷大(n -> ∞)：
在这里插入图片描述

题目地址

http://acm.hdu.edu.cn/showproblem.php?pid=1018

题目描述

在这里插入图片描述

思路

当n很小时，直接取 lg n!
利用Stirling公式计算n!结果的位数时，则是两边取对数

故n!的位数为:

代码

语言支持：C++

C++ Code:

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const double pi = acos(-1.0), e = exp(1.0);
int main(){
	int n, T;
	scanf("%d", &T);
	while (T--){
		scanf("%d", &n);
		double temp = 0;
		if (n <= 100){ // lg(n!) = lg1 + lg2 + ... + lgn 
			for (int i = 1; i <= n; i++)
				temp += log10(i);
		}
		else // Stirling公式
			temp = floor(log10(2 * pi * n) / 2 + n * log10(n / e));
		printf("%lld\n", (long long)temp + 1);
	}
	return 0;
}