power oj 1787 数论（一）：求n的阶乘的位数

今天又是被数学教做人的一天。

话不多说，上题！

---

Description
计算n！的位数。

Input
多组数据输入，每次输入一个自然数n，n<=10^9

Output
n！的位数

Sample Input
10
20

Sample Output
7
19

---

求n的阶乘的位数，暴力绝对解不了。n的范围高达1e9，暴力肯定爆。

另外，用不断向下除10的方法也绝对不行。因为，即使是用long long，最多可以计算的位数是19；

这道题，只可以用斯特林公式：

可是这个公式求出来的是n阶乘的大小，不能求位数。

我看到一个厉害的博主，直接用了log10（2pin）/2.0+n*log10(n/e)+1;

其实，就是用的log10可以求位数的特性。

比如说，log10（1000）=3
log10（22）=1
都加上1之后，就是log10（n）的位数。

所以，log10（斯特林公式），用对数化简一下，就得到了这个公式。

还有，这个式子算不来1的阶乘

为什么呢？因为当n很小的时候，斯特林公式计算阶乘还是有点误差的。

所以，1的话，需要另外说明。

下面粘上ac代码：

#include<iostream>
#include <stdio.h>
#include<math.h>

#define e 2.718281828459
#define pi 3.1415926535898
typedef long long ll;
using namespace std;

int main()
{

    ll n;
    while(scanf("%d",&n)!=EOF)
    {
        if(n==1||n==0)
        {
        	printf("1\n");
        	continue;//Do you like higer brother?
		}
		if(n==-1)
		{
			printf("0\n");//When you say you like Higer too;
        	continue;
		}
        ll ans;
        ans=log10(2.0*pi*n)/2.0+n*log10(n/e)+1;//I am going crazy with happiness.
        printf("%lld\n",ans);
    }
    return 0;
}

好的，就是这样。

今天这道题目的收获：

log10（n^n）,如果数学差一些，就会想到用快速幂，其实，可以直接利用对数的性质，化为n*log(n)
。这样极大的减少了计算量。

所以，数学公式的用法，还是要掌握熟练的。

行了，继续加油吧。h