给一个数X,len=log10(X)+1就是X这个数的长度
当N的值不超过10^6时
当N的值不超过10^6时
那么N的阶乘长度=log10(1*2*3*4*……*N)+1=log10(1)+log10(2)+log10(3)+log10(4)+……+log10(N)+1;
当然当N比较大时,时间复杂度也很高,这时有
斯特林公式N!=sqrt(2*PI*N)*(N/e)^N (其中PI=3.141592653,e=2.718281828459)
N!的长度len=log10(sqrt(2*PI*N)*(N/e)^N)+1 =0.5*log10(2.0*PI*N)+N*log10(N/e)+1;
输入任意一个数,阶乘位数代码如下:
//斯特林公式n!=sqrt(2*PI*n)*(n/e)^n
#include <stdio.h>
#include <algorithm>
#include <math.h>
#include <iostream>
using namespace std ;
#define PI 3.1415926
#define e 2.718281828459
int main()
{
__int64 len,n;
scanf ( "%I64d", &n ) ;
if(n == 1) printf("1\n");
else
{
len = 0.5*log10(2.0*PI*n)+n*log10(n*1.0/e)+1;
printf("%I64d\n",len) ;
}
return 0;
}