51nod 1130 N的阶乘的长度（阶乘近似公式--斯特林公式）

大致题意：

求n！有多少位，n<= 1e9 ，1000组数据

大致思路：

显然是企图求 lg（n！）+1 向下取整

显然求lg(n!) 展开的复杂度是On的

但是，他是求多少位，有n!的近似公式斯特林公式（Stirling's approximation）

维基百科：https://en.wikipedia.org/wiki/Stirling%27s_approximation

显然对于求对数函数来说精确度在n比较大的情况下必然满足。

//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <iostream>
#include <cstring>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <string>
#include <vector>
#include <cstdio>
#include <ctime>
#include <bitset>
#include <algorithm>
#define SZ(x) ((int)(x).size())
#define ALL(v) (v).begin(), (v).end()
#define foreach(i, v) for (__typeof((v).begin()) i = (v).begin(); i != (v).end(); ++ i)
#define reveach(i, v) for (__typeof((v).rbegin()) i = (v).rbegin(); i != (v).rend(); ++ i)
#define REP(i,n) for ( int i=1; i<=int(n); i++ )
#define rep(i,n) for ( int i=0; i< int(n); i++ )
using namespace std;
typedef long long ll;
#define X first
#define Y second
#define PB push_back
#define MP make_pair
typedef pair<int,int> pii;

template <class T>
inline bool RD(T &ret) {
    char c; int sgn;
    if (c = getchar(), c == EOF) return 0;
    while (c != '-' && (c<'0' || c>'9')) c = getchar();
    sgn = (c == '-') ? -1 : 1;
    ret = (c == '-') ? 0 : (c - '0');
    while (c = getchar(), c >= '0'&&c <= '9') ret = ret * 10 + (c - '0');
    ret *= sgn;
    return 1;
}
template <class T>
inline void PT(T x) {
    if (x < 0) {
        putchar('-');
        x = -x;
    }
    if (x > 9) PT(x / 10);
    putchar(x % 10 + '0');
}

const double pi = acos(-1);

int main(){

        int T;
        RD(T);
        while( T--){
                ll n;
                RD(n);
                ll ans = log10(2*pi*n)/2 + n*(log10(n/exp(1.0)));
                PT(ans+1); puts("");
        }
}

Input

第1行：一个数T，表示后面用作输入测试的数的数量。（1 <= T <= 1000)
第2 - T + 1行：每行1个数N。（1 <= N <= 10^9)

Output

共T行，输出对应的阶乘的长度。

Input示例

3
4
5
6

Output示例

2
3
3