大数阶乘的位数和精确值计算
2017-12-30 19:18
239 查看
转自:https://www.cnblogs.com/stonehat/p/3603267.html
我们知道整数n的位数的计算方法为:log10(n)+1
故n!的位数为log10(n!)+1
如果要求出n!的具体值,对很大的n(例如n=1000000)来说,计算会很慢,如果仅仅是求阶乘的位数,可以用斯特林(Stirling)公式求解
斯特林(Stirling)公式:
于是求n!的位数就是求log10((2*PI*n)^1/2*(n/e)^n)+1
即 1/2*log10(2*PI*n)+n*log10(n/e)+1
所以采用下面代码计算阶乘位数,会非常快
如果要计算阶乘的精确值,则可以采用下面代码。
复制代码
n: n 的阶乘
返回值: 阶乘结果的位数
注意:
本程序直接输出n!的结果,需要返回结果请保留long a[]
需要 math.h
我们知道整数n的位数的计算方法为:log10(n)+1
故n!的位数为log10(n!)+1
如果要求出n!的具体值,对很大的n(例如n=1000000)来说,计算会很慢,如果仅仅是求阶乘的位数,可以用斯特林(Stirling)公式求解
斯特林(Stirling)公式:
于是求n!的位数就是求log10((2*PI*n)^1/2*(n/e)^n)+1
即 1/2*log10(2*PI*n)+n*log10(n/e)+1
所以采用下面代码计算阶乘位数,会非常快
#define PI 3.141592654 #define E 2.71828182846 int l(int n) { int s=1; if(n>3) s=log10(2*PI*n)/2+n*log10(n/E)+1; return s; }
如果要计算阶乘的精确值,则可以采用下面代码。
复制代码
n: n 的阶乘
返回值: 阶乘结果的位数
注意:
本程序直接输出n!的结果,需要返回结果请保留long a[]
需要 math.h
int factorial(int n) { long a[10000]; int i,j,l,c,m=0,w; a[0]=1; for(i=1;i<=n;i++) { c=0; for(j=0;j<=m;j++) { a[j]=a[j]*i+c; c=a[j]/10000; a[j]=a[j]%10000; } if(c>0) {m++;a[m]=c;} } w=m*4+log10(a[m])+1; printf("\n%ld",a[m]); for(i=m-1;i>=0;i--) printf("%4.4ld",a[i]); return w; }
相关文章推荐
- 大数阶乘的位数和精确值计算
- 大数阶乘的位数和精确值计算
- 大数阶乘的位数和精确值计算【转】
- 计算一个大数n的阶乘的位数宽度(十进制)(log i累加法 )
- C#阶乘类,可以精确计算大数的阶乘,10000的阶乘只要1.3秒左右
- 菜鸟的大数阶乘的位数公式计算
- POJ 1423 Big Number 大数阶乘的位数计算
- C#阶乘类,可以精确计算大数的阶乘,10000的阶乘只要1.3秒左右
- POJ 1423 Big Number 大数阶乘的位数计算
- 计算一个大数n的阶乘的位数宽度(十进制)(log i累加法 )
- 计算一个大数n的阶乘的位数宽度(十进制)转载
- POJ-1423 计算出n的阶乘的位数大数问题[Stirling公式]
- 阶乘类,可以精确计算大数的阶乘,10000的阶乘只要1.3秒左右
- 大数,高精度计算---大数阶乘
- 计算大数阶乘--Linux Bash版
- 大数,高精度计算---大数阶乘
- 大数阶乘的计算(六)
- N的阶乘的长度 V2(斯特林近似) 求 某个大数的阶乘的位数 .
- [ACM Steps]2.2.8 Big Number 求N!位数 大数N阶乘的位数
- 阶乘位数计算