利用斯特林(Stirling)公式的进行求解大数阶乘的位数。下面是推导得到的公式:res=(long)( (log10(sqrt(4.0*acos(0.0)*n)) + n*(log10(n)-log10(exp(1.0)))) + 1 ); 当n=1的时候,上面的公式不适用,所以要单独处理n=1的情况!有关斯特林(Stirling)公式及其相关推导,这里就不进行详细描述,有兴趣的话自己可以做更多了解。这种方法速度很快就可以得到结果。
具体看一道例题:只有当数据很小的时候,才可以用下面的第一种方法做
Description
Fresher自从学会阶乘概念以后,越来越着迷,每天都想着阶乘的各种巧妙应用及其求解。有一天,Fresher正在和要好的同学讨论天底下最大的数据可能有多少位?于是就有一位数学较好的同学提出了一个问题:当N较大时,虽然求解其阶乘的具体值比较困难,但我们应该有方法知道N!的量级吧,也就是说N!在十进制表示法下有多少位。假设现在需要求解N!的十进制位数,您可以帮Fresher想个奇妙的方法求解吗?
Input
一个正整数m(m<=1000000)
Output
一个正整数n,代表m!的十进制表达法的位数
Sample Input
5
Sample Output
3
这个方法只能解决小一点的数,因此不推荐使用这种方法
#include<stdio.h>
#include<math.h>
int main()
{
int n,m,i;
double sum=0.0;
scanf("%d",&m);
for(i=1;i<=m;i++)
sum+=log(i)/log(10);
//C语言中没有以任意为底数的对数函数,所以log10是没有的.所以可以用
logx(y)=log(y)/log(x)来代替
printf("%d\n",(int)sum+1);
return 0;
}
但是当数据很大时,上面的方法就不能用了,下面有更厉害的方法