hdu1018 Big Number stirling公式

最新推荐文章于 2019-02-26 20:56:29 发布

落雨飞辰

最新推荐文章于 2019-02-26 20:56:29 发布

阅读量505

点赞数 1

CC 4.0 BY-SA版权

分类专栏： ACM 文章标签： Big Number

如有错误，欢迎指正批评 && 转载请注明谢谢

本文链接：https://blog.youkuaiyun.com/bigsungod/article/details/47296883

ACM 专栏收录该内容

70 篇文章

订阅专栏

本文介绍了斯特林公式用于近似计算阶乘的方法，并通过代码实现了一个函数来计算对数形式的阶乘近似值。适用于需要快速估算大数阶乘场合。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Stirling公式：n!与sqrt(2πn) * n^n * e^(-n)的值十分接近

所以log10(n!) = log(n!) / log(10) = ( n*log(n) - n + 0.5*log(2*π*n))/log(n);

代码：

#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstdlib>
#include <cstring>
#include <iomanip>
#include <cstdio>
#include <string>
#include <bitset>
#include <vector>
#include <queue>
#include <stack>
#include <cmath>
#include <list>
#include <map>
#include <set>
#define sss(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define mem1(a) memset(a,-1,sizeof(a))
#define mem(a) memset(a,0,sizeof(a))
#define ss(a,b) scanf("%d%d",&a,&b)
#define s(a) scanf("%d",&a)
#define p(a) printf("%d\n", a)
#define INF 0x3f3f3f3f
#define w(a) while(a)
#define PI acos(-1.0)
#define LL long long
#define eps 10E-9
#define N 100000
#define mod 100000000
using namespace std;
void mys(int& res)
{
    int flag=0;
    char ch;
    while(!(((ch=getchar())>='0'&&ch<='9')||ch=='-'))
        if(ch==EOF)  res=INF;
    if(ch=='-')  flag=1;
    else if(ch>='0'&&ch<='9')  res=ch-'0';
    while((ch=getchar())>='0'&&ch<='9')  res=res*10+ch-'0';
    res=flag?-res:res;
}
void myp(int a)
{
    if(a>9)
        myp(a/10);
    putchar(a%10+'0');
}
/*************************THE END OF TEMPLATE************************/
double stirling(int n){
    return ceil((n*log(double(n))-n+0.5*log(2.0*n*PI))/log(10.0));
}
int main()
{
	int cas,n;
	scanf("%d",&cas);
	while(cas--)
	{
		scanf("%d",&n);
		if(n<=1)printf("1\n");
		else printf("%.0lf\n",stirling(n));
	}
	return 0;
}