POJ 1423 计算阶乘的位数

#include<cstdio>#include<cstring>#include<iostream>#include<cmath>#define PI 3.141592654#define E 2.71828182846using 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))