POJ 1423 计算阶乘的位数

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#define PI 3.141592654
#define E 2.71828182846
using namespace std;
int Len(int n)
{
    int s=1;
    if(n>3) s=log10(2*PI*n)/2+n*log10(n/E)+1; //斯特林公式，估算阶乘位数
    return s;
}
int main()
{
    int i,j,n,t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        n=Len(n);
        printf("%d\n",n);
    }
    return 0;
}
//n!的位数为 log10(2*PI*n)/2+n*log10(n/E)+1（注意：当n=1时，算得的结果为0）
//高德纳的《计算机程序设计艺术》中，
//n! = sqrt(2*π*n) * ((n/e)^n) * (1 + 1/(12*n) + 1/(288*n*n) + O(1/n^3))