数学:题目69.数的长度

 
 
/* NYOJ69 阶乘数位长度
 * 方法一:
 * 可设想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比较大的时候,这种方法方法需要花费很多的时间。
 * 
 * 方法二:
 * 利用斯特林(Stirling)公式的进行求解。下面是推导得到的公式:
 * res=(long)( (log10(sqrt(4.0*acos(0.0)*n)) + n*(log10(n)-log10(exp(1.0)))) + 1 );
 * 当n=1的时候,上面的公式不适用,所以要单独处理n=1的情况!
 * 有关斯特林(Stirling)公式及其相关推导,这里就不进行详细描述,有兴趣的话可看这里。
 * 这种方法速度很快就可以得到结果。详细证明如下:
 * http://episte.math.ntu.edu.tw/articles/mm/mm_17_2_05/index.html
*/
#include<iostream>
#include <cmath>
using namespace std;
int normal(double n)
{
 double x=0;
 while(n)
 {
  x +=log10(n);
  n--;
 }
 return (int)x+1;
}
long stirling(double n)
{
 long x=0;
 if( n ==1 )
  x = 1;
 else
 {
  x = (long)( (log10(sqrt(4.0*acos(0.0)*n)) + n*(log10(n)-log10(exp(1.0)))) + 1 );
 }
 return x;
}
int main()
{
 int n;
 cin>>n;
 while(n--)
 {
  int x;
  cin>>x;
  cout<<stirling(x)<<endl;
 }
 return 0;
}      

此为转载自 http://acm.nyist.net/JudgeOnline/bestcode.php?pid=69        

转载于:https://www.cnblogs.com/maple-2013/archive/2013/04/23/3038054.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值