poj 1423 stirling公式（阶乘的位数）

题意：给定一个数（小于10^7），求这个数阶乘的位数。

思路：首先，求一个数n的位数可以用log10（n）向上取整来得到（10的整数次幂要特殊判断）。所以所求由log10（1）+log10（2）+log10（3）+……+log10（n）可得。n很大，如此速度太慢，用stirling公式求n！。

斯特林公式（Stirling's approximation）是一条用来取n的阶乘的近似值的数学公式。n！近似等于（2*pi*n)^1/2  * (n/e)^n *e

gcc调用math库还要加上 -lm选项

用下面注释掉的那行输出用c提交AC，用gcc便wa……

#include <stdio.h>
#include <string.h>
#include <math.h>
#define E exp(1)
#define PI acos(-1)
int T,n;
int main(){
	freopen("a.txt","r",stdin);
	scanf("%d",&T);
	while(T--){
		double a,res;
		scanf("%d",&n);
		if(n <= 3){
			printf("1\n");
			continue;
		}
		a = sqrt(2*PI*n);
		res = log10(a)+n*(log10(1.0*n)-log10(E));
		printf("%d\n",(int)res+1);
		//printf("%.0lf\n",ceil(res));
	}
	return 0;
}