bzoj3000 Big Number 数论，斯特林公式

Description

给你两个整数N和K，要求你输出N！的K进制的位数。

Input

有多组输入数据，每组输入数据各一行，每行两个数——N，K

Output

每行一个数为输出结果

Sample Input

2 5
2 10
10 10
100 200

Sample Output

1
1
7
69
对于100%的数据，有2≤N≤2^31， 2≤K≤200，数据组数T≤200。

题解

用Stirling公式求近似值

位数＝logk(n!)+1

≈ logk(sqrt(2πn)*(n/e)^n)+1

= logk( sqrt(2πn))+log[(n/e)^n]+1

=1/2*logk( 2πn)+nlog(n/e)+1

=0.5*logk ( 2πn)+nlog(n/e)+1

=0.5*logk ( 2πn)+nlog(n)-nlog(e)+1

PS:pi=acos(-1.0),e=exp(1)

PS2:eps的存在是为了防止n=2,k=2这样刚好的情况出现，这个时候向上取整要多取1位

斯特林公式是求解n!的近似解，对于较大的n数值是十分准确的。

所以可以通过数学方法解决。

 #include<cstring>
 #include<cmath>
 #include<iostream>
 #include<cstdio>
 #include<algorithm>

 #define ll long long
 using namespace std;
 const double eps=0.00000000001;
 const double pai=3.14159265359;
 const double e=exp();

 int n,k;

 int main()
 {
     freopen("fzy.in","r",stdin);
     freopen("fzy.out","w",stdout);
     while(~scanf("%d%d",&n,&k))
     {
         if (n<=)
         {
             double ans=;
             for (int i=;i<=n;i++)
                 ans+=log(i);
             ans/=log(k);
             int res=ceil(ans+eps);
             printf("%d\n",res);
         }
         else
         {
             double res=(log(sqrt(*pai*n))+n*log(n/e))/log(k);
             ll ans=ceil(res-eps);
             printf("%lld\n",ans);
         }
     }
 }

对了，c++小数处理的时候会有精度损失的问题，所以需要适当加上一个小数