POJ 1423

问题描述：输入一个整数n（1，10^7），要你输出该整数的阶乘后的位数是多少？

     问题分析：由于n的输入数可能较大，倘若按照常规算法（就是算出n的阶乘）设计的话，一方面会出现TLE，另一方面会出现数据溢出的情况。所以我们不得不充分利用数学知识里面的对数与指数间的亲密关系。

 

但 是有没有办法不求这个阶乘，所以得用到stirling逼近    n!=(n/e)^n*sqrt(2nπ) n→∞     ①

不得不说这是个很牛叉的公式，在n=1时就有很高的精度，最多相差不过0.2 

       但问题是时间解决了，数据范围没有解决。。。。

       补充知识 ：若n=a+10^b,(a<1,b∈N)，则n的位数为b。所以要用到lg神器。。。。lgn=b+lg(a)，所以说要求b只需求lgn+1(由于lga<0所以要+1)即可

      因而利用对数的运算法则把①式改写就可以了

上代码

#include <iostream>
#include <cmath>
using namespace std;
const double e=2.718281828459;
const double pi=3.141582626535;

double f(int n)
{
	return n*(log10(n)-log10(e))+0.5*(log10(2)+log10(n)+log10(pi));
}

int main()
{
	int t,m,s;
	cin>>t;
	while (t--)
	{
		cin>>m;
		s=(int)f(m)+1;
		cout<<s<<endl;
	}
	return 0;
}