Stirling公式求阶乘位数 51nod1058 poj1423

输入N求N的阶乘的10进制表示的长度。例如6! = 720，长度为3。

Input

输入N(1 <= N <= 10^6)

Output

输出N的阶乘的长度

Input示例

6

Output示例

3        

公式 

log10(sqrt(2*acos(-1.0)*n))+n*(log10(n/exp(1.0)))+1

   

注意精度 

或者 π和e acos(-1.0)和exp(1.0)

 注意当n=1的时候公式不适用

    poj代码

#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;
int main()
{
	int n;
	int t;
	cin>>t;
	double p=3.141592653589793239,e=2.7182818284590452354;
	for(int i = 0;i<t;i++){
		cin>>n;
		if(n==1){
			cout<<1<<endl; continue;
		}
		double len=0.5*log10(2*p*n)+n*log10(n/e)+1; 
		long long ans = len;
		printf("%lld\n",ans);}
	return 0;
}