这一道题一开始看着挺简单的,就是求阶乘然后再求位数。但是理解题意之后发现这个想法有一个很大的问题就是当输入数字很大时就会出现问题 int long 之类的东西无法存储这些数字。所以在网上搜索了一下发现还有一个斯特林公式可以求解这个问题。
斯特林公式
根据上述公式可以得到位数的求解
lnN!=NlnN-N+0.5*ln(2*N*pi)
要想求有多少位,将他换成以10为底便可!
利用换底公式得
lnN!/ln10=log10N!
最后再加上 1就可以得到结果了。
#include<iostream>
using namespace std;
#define pi acos(-1.0)
int main()
{
int m, n;
cin >> m;
while (m--)
{
int t = 0;
cin >> n;
t=(int)((n*log(n) - n + (log(2 * pi*n)) / 2) / log(10)) + 1;
cout <<t<<endl;
}
return 0;
这道题基本也就结束了。最后在网上发现一个大牛xiaonian简单了分析了一下整个原理最后发现原来还可以这么简单的理解这个问题。下面是它的分析过程。
在这之前,我们必须要知道一个知识,任意一个正整数a的位数
我们知道了一个正整数a的位数等于(int)log10(a) + 1,
现在来求n的阶乘的位数:
假设A=n!=1*2*3*……*n,那么我们要求的就是
(int)log10(A)+1,而:
log10(A)=log10(1*2*3*……n) (根据log10(a*b) = log10(a) + log10(b)有)
=log10(1)+log10(2)+log10(3)+……+log10(n)
现在我们终于找到方法,问题解决了,我们将求n的阶乘的位
数分解成了求n个数对10取对数的和,并且对于其中任意一个数,
都在正常的数字范围之类。
总结一下:n的阶乘的位数等于
(int)(log10(1)+log10(2)+log10(3)+……+log10(n)) + 1
根据这个思路我们很容易写出程序
#include<iostream>
using namespace std;
#define pi acos(-1.0)
int main()
{
int m, n;
cin >> m;
while (m--)
{
double t = 0;
cin >> n;
for (int j = 1; j <= n; j++)
{
t += log10(j);
}
cout << int(t)+1<<endl;
}
return 0;
}