哎~~~开始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;
}