求大数阶乘的位数ACM Big Number问题

哎~~~开始A题玩玩，ACM还是蛮有味么~~

问题描述：求一个大数阶乘的位数

解决策略：

方法一：log10(n!)=log10(1*2*3…*n)=log10(1)+log10(2)+…+log10(n)+1

即对log10(n!)的值取整加1就是n!的位数——朴素方法。

//Author:YunMengZe
//DateTime:2013.10.30
//Description:求给定数阶乘的位数
#include <iostream>
#include <cmath>
using namespace std;
int main()
{
	int n;
	double d=0;
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		d+=log10(double(i));
	}
	cout<<"n阶乘的位数位: "<<(int)(d+1)<<endl;
	return 0;
}

 

方法二：近似计算的方法，斯特林公式。n!近似等于sqrt(2*pi*n) * (n/e)^n

//Author:YunMengZe
//DateTime:20131030
//Method:斯特林公式 n！=sqrt(2*pi*n)*(n/e)^n
#include <iostream>
#include <cmath>
using namespace std;
#define PI 3.1415926
#define e  2.71828
int main()
{
 int n;
 double d=0;
 cin>>n;
 if(n==1||n==0)
  d=1;
 else
 {
  d=log10(double(2*PI*n))/2 + n*log10(double(n/e));
  d=ceil(d); 
 }
 cout<<"n阶乘的位数为："<<d<<endl;
 return 0;
}