C++斯特林公式求n!的位数

最新推荐文章于 2024-07-24 20:25:39 发布

SkeletonKing233

最新推荐文章于 2024-07-24 20:25:39 发布

阅读量784

点赞数 2

文章标签： C++ 算法斯特林公式

本文链接：https://blog.csdn.net/SkeletonKing233/article/details/101313184

版权

斯特林公式简介

斯特林公式是一条用来取n的阶乘的近似值的数学公式。一般来说，当n很大的时候，n阶乘的计算量十分大，所以斯特林公式十分好用，而且，即使在n很小的时候，斯特林公式的取值已经十分准确。

形式为 n! ≈ sqrt(2πn) * (n / e) ^ 2
等式两边同时log10的话得到：
log10(n!) ≈ log10(2πn) / 2 + n * log10(n / e)

这样n!的位数就是log10(n!)之后向上取整，也就是log10(2πn) / 2 + n * log10(n / e)向上取整了。这样原本可能计算n!的位数时需要On级别的算法，n大了还要用到高精度。可是有了这个斯特林公式，不仅时间复杂度减少为了O1级别，代码也减少了很多行。

但这里请注意，如果说n = 0或者n = 1时，用斯特林公式计算出来的值是不准确的，所以要加一个特判，既 if (n == 0 || n == 1) ans = 1。

斯特林公式代码

代码直接把公式抄上去就OK了。

# include <cstdio>
# include <iostream>
# include <cmath>
# include <cstring>
# include <algorithm>
using namespace std;

# define FOR(i, a, b) for(int i = a; i <= b; i++)
# define _FOR(i, a, b) for(int i = a; i >= b; i--)

const double E = 2.7182818284;
const double PI = 3.1415926535;

int n;

int main()
{
	scanf("%d", &n);
	int ans = ceil(log10(2 * PI * n) / 2 + n * log10(n / E));
	if(n == 0 || n == 1) ans = 1;
	printf("%d\n", ans);
	return 0;
}