斯特林公式

原创 2015年11月21日 14:59:42

斯特林[striling]公式(求阶乘(n!)的位数)

斯特林公式是一条用来取n阶乘近似值的数学公式,由公式
可得n!的位数为0.5*log10(2*pi)+(m+0.5)*log10(m)-m*log10(e)+1;
另外还有一个计算方法,n!=1*2*3*...*n;

所以log10(n!)=log10(1*2*3*...*n);
即log10(n!)=log10(1)+log10(2)+log10(3)+...+log10(n);
n!的位数为log10(n!)+1。

/*例如1000阶乘位数:
log10(1)+log10(2)+···+log10(1000)取整后加1
*/
#include
<stdio.h>

#include<math.h>


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


#include<stdio.h>
#include<math.h>
#define PI 3.14159265
int main(){
int len,N;
while(scanf("%d",&N)!=EOF)
{
if(N==1)
len=1;
else
len=(int)ceil((N*log(N)-N+log(2*N*PI)/2)/log(10));////ceil求上界,即不小于某值的最小整数
//string公式lnN!=NlnN-N +0.5*ln( 2*N*pi)
//而N次方阶乘的位数等于:
//log10(N!)取整后加1
//log10(N!)=lnN!/ln(10)

//ceil为求上界,即不小n的最小整数
//log取自然对数
printf("%d\n",len);
}
return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

HDU-1018-Big Number( 这段数学推导我给满分! && 然而并不想用斯特林公式! )

Big Number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Tot...

HDU 1018 斯特林公式应用

链接:http://acm.hdu.edu.cn/showproblem.php?pid=1018 题意:输入一个整数,输出这个数的阶乘的位数 思路:    ...

HDU1018(斯特林公式)

Problem Description In many applications very large integers numbers are required. Some of these ...

hdu 1018 Big Number (斯特林公式)

(先看题目) Big Number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Ot...

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

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

斯特林引擎小车制作

  • 2016-02-04 22:47
  • 9.28MB
  • 下载

Cpp环境【Vijos1060】斯特林数:盒子与球

【问题描述】   n 个盒子排成一行(编号为1..n)。你有A个红球和B个蓝球。球除了颜色没有任何区别。你可以将球放进盒子。一个盒子可以同时放进两种球,也可以只放一种,也可以空着。球不必全...

斯特林发动机模型

  • 2013-05-19 12:02
  • 135KB
  • 下载

hdu3625 Examining the Rooms && hdu4372 Count the Buildings(第一类斯特林数)

先吐槽下,这两道题的数据很坑。。 参考:Acdream,之所以这里写出来只是为了让自己巩固,没别的意思。。 首先解释第一类斯特灵数:将p个物体排成k个非空循环排列的方法数。s(p,0)=0...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)