斯特林公式在 ACM中的使用

*斯特林公式：

*求阶乘的第（1~14）位：

frac(n):取n的小数部分；

这个式子精确度比较高，精确度达到第14位；不过前几个数字需要手动打表算出，因为当n较小时，答案不符合；

详细的推导在链接：https://share.weiyun.com/df75bafcd0832c034b22fe2355a86a31 （密码：CWPP）（感谢作者）

*求阶乘的位数：

最后计算所得结果+1；

第一个题，求阶乘的第二位。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define pi acos(-1.0)
#define e exp(1.0)
using namespace std;
typedef long long ll;

int main()
{
//    printf("%lf\n",pi);
    int m;
    int T[11]={0,0,0,0,4,2,2,0,0,6,6};
    scanf("%d",&m);
    while(m--)
    {
        double n;
        double s;
        scanf("%lf",&n);
        if(n<11)
        {
            printf("%d\n",T[(int)n]);
            continue;
        }
        s=0.5*log10(2*n*pi)+n*log10(n)-n*(log10(e));
        s=s-(int)s;
        s=pow(10,s);
        s=s*exp(1.0/12.0/n-1.0/360.0/n/n/n);
        printf("%d\n",((int)(s*10))%10);
    }
    return 0;
}

第二题：判断阶乘的位数。

#include<cstdio>
#include<cstring>
#include<cmath>
#define pi acos(-1.0)
#define e exp(1.0)//exp(x)是求e的x次方
int main()
{
    int t,n;
    double ans;
    while(~scanf("%d",&n))
    {
        ans=log10(sqrt(2*pi*n))+n*log10(n/e);
        printf("%d\n",(int)ans+1);
    }
    return 0;
}

//over!