【nowcoder 算法训练3】不凡的夫夫（斯特灵公式，阶乘的位数）

链接：https://www.nowcoder.net/acm/contest/75/A
来源：牛客网

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 32768K，其他语言65536K
64bit IO Format: %lld
题目描述
夫夫有一天对一个数有多少位数感兴趣，但是他又不想跟凡夫俗子一样，
所以他想知道给一个整数n，求n！的在8进制下的位数是多少位。
输入描述:
第一行是一个整数t(0<t<=1000000)(表示t组数据)
接下来t行，每一行有一个整数n(0<=n<=10000000)
输出描述:
输出n！在8进制下的位数。
示例1
输入
3
4
2
5
输出
2
1
3

分析：斯特灵公式用于求n!

利用Stirling公式求解n!的位数：易知整数n的位数为[lgn]+1。利用Stirling公式计算n!结果的位数时，可以两边取对数，得:

n!的位数为：

斯特灵公式见：http://blog.csdn.net/liangzhaoyang1/article/details/51145807

求n！在8进制下的位数 故用换底公式：
logn(m)=log(n)/log(m) log(n) 是以e为底的对数函数

注意 n==0 时 特判

代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <map>
#include <set>
#include <bitset>
#include <cctype>
#include <cstdlib>
#include <queue>
#include <cmath>
#include <stack>
#include <ctime>
#include <string>
#include <vector>
#include <sstream>
#include <functional>
#include <algorithm>
using namespace std;

#define mem(a,n) memset(a,n,sizeof(a))
#define memc(a,b) memcpy(a,b,sizeof(b))
#define rep(i,a,n) for(int i=a;i<n;i++) ///[a,n)
#define dec(i,n,a) for(int i=n;i>=a;i--)///[n,a]
#define pb push_back
#define fi first
#define se second
#define IO ios::sync_with_stdio(false)
#define fre freopen("in.txt","r",stdin)
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
typedef long long ll;
typedef unsigned long long ull;
const double PI=acos(-1.0);
const double e=2.718281828459045;
const double eps=1e-8;
const int INF=0x3f3f3f3f;
const int MOD=258280327;
const int N=1e5+5;
const ll maxn=5e4;
const int dir[4][2]= {-1,0,1,0,0,-1,0,1};
///换底公式 logn(m)=log(n)/log(m)  log(n) 是以e为底的对数函数

/***
斯特灵公式：求n!的近似值的数学公式
 n! ~ sqrt(2*PI*n)*(n/e)^n

log10(n!) = log10(2*n*PI)/2+n*log10(n/e)
则答案为：
ans = log10(2*n*Pi)/2+n*log10(n/e) + 1
***/
int main()
{
    int T,n;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        if(n==0) {puts("1");continue;}
        double tmp=log(2*PI*n)/log(8)/2+n*log(n/e)/log(8);
        printf("%d\n",int(tmp)+1);
    }
    return 0;
}