poj-1423 NYOJ_69 数字长度斯特林公式对数应用-CSDN博客

数的长度

时间限制： 3000 ms | 内存限制： 65535 KB

难度： 1

描述

N！阶乘是一个非常大的数，大家都知道计算公式是N!=N*(N-1）······*2*1.现在你的任务是计算出N！的位数有多少（十进制）？

输入

首行输入n，表示有多少组测试数据(n<10)
随后n行每行输入一组测试数据 N( 0 < N < 1000000 )

输出

对于每个数N，输出N!的（十进制）位数。

样例输入

样例输出

1
1
130271



刚拿到这个题，因为是在数学分类里的，我想着是不是又公式呀！！一般阶乘的题都有特殊方法，首先想到的是用数组存 高位数 的方法，但是当大于5000就开始慢了，要超时。没办法，看了留言
板，说是 可以用对数 或者斯特林公式。


对数的方法：
可设想n!的结果是不大于10的M次幂的数,即n!<=10^M(10的M次方),则不小于M的最小整数就是 n!的位数，对 该式两边取对数，有 M =log10^n! 即：M = log10^1+log10^2+log10^3...+log10^n 循环求和,就能算得M值，该M是n!的精确位数。当n比较大的时候，这种方法方法需要花费很多的时间。

View Code

#include <stdio.h>
#include <math.h>
int main()
{
int n,N,i;
double sum;
scanf("%d",&N);
while (N--)
{
scanf("%d",&n);
sum=1;
for(i=1;i<=n;i++)
sum+=log10((double)i);
printf("%d\n",(int)sum);
}
return 0;
}

斯特林公式：

log(n!) = log10(sqrt(2*pi*n)) + n*log10(n/e);

其中pi是圆周率，e是自然对数。

1时要特殊处理：

View Code

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cmath>
#define e 2.718281828459045
#define pi 3.141592653589793239
using namespace std;
int main ()
{
    int cas,n;
    scanf("%d",&cas);
    while (cas --)
    {
          scanf("%d",&n);
          double t = log10(sqrt(2*pi*n)) + n * log10(n/e);
          printf ("%d\n",(int)t + 1);
    }
return 0;
}