Расчет больших факториалов (до 100000 !!!!!)

← →
Alx2 (2002-02-15 10:53) [7]

Вот приближение десятичного логарифма факториала:
#include <math.h> double log10factorial(n) double n; { double t1; double t15; double t17; double t2; double t23; double t25; double t3; { t1 = n*n; t2 = t1*n; t3 = t1*t1; t15 = log(902961561600.0*t3*t2+75246796800.0*t3*t1+3135283200.0*t3*n -2421135360.0*t3-207204480.0*t2+707957280.0*t1+62961828.0*n-534703531.0); t17 = log(n); t23 = log(0.3141592653589793E1/0.4076697908635553E24); t25 = log(10.0); return((2.0*t15-13.0*t17+2.0*n*t17-2.0*n+t23)/t25/2.0); } }
Естественно, целая часть этого покажет количество знаков числа+1, а с мантиссой можно работать.
Исходная формула (аппроксимация факториала):
n! ~ 1/902961561600*sqrt(2)*sqrt(Pi)*(902961561600*n^7+75246796800*n^6+3135283200*n^5-2421135360*n^4-207204480*n^3+707957280*n^2+62961828*n-534703531)*n^(-13/2+n)*exp(-n)

Расчет больших факториалов (до 100000 !!!!!) Найти похожие ветки