题意:给你两个整数N和K,要求你输出N!的K进制的位数。
题解:首先补了一下对数的运算法则。
1.
loga(mn)=loga(m)+loga(n)
l
o
g
a
(
m
n
)
=
l
o
g
a
(
m
)
+
l
o
g
a
(
n
)
2.
loga(mn)=loga(m)−loga(n)
l
o
g
a
(
m
n
)
=
l
o
g
a
(
m
)
−
l
o
g
a
(
n
)
3.
loga(mn)=nloga(m)
l
o
g
a
(
m
n
)
=
n
l
o
g
a
(
m
)
4.
loga(m−−√n)=1nloga(m)
l
o
g
a
(
m
n
)
=
1
n
l
o
g
a
(
m
)
再有一个
n!
n
!
的近似值公式,就是斯特林公式。
n!≈2πn−−−√(ne)n
n
!
≈
2
π
n
(
n
e
)
n
于是我们就可以求答案了。
ans=logk(n!)+1
a
n
s
=
l
o
g
k
(
n
!
)
+
1
=>
ans=logk(2πn−−−√(ne)n)+1
a
n
s
=
l
o
g
k
(
2
π
n
(
n
e
)
n
)
+
1
=>
ans=logk(2πn−−−√)+logk((ne)n)+1
a
n
s
=
l
o
g
k
(
2
π
n
)
+
l
o
g
k
(
(
n
e
)
n
)
+
1
=>
ans=12logk(2πn)+nlogk(ne)+1
a
n
s
=
1
2
l
o
g
k
(
2
π
n
)
+
n
l
o
g
k
(
n
e
)
+
1
注意当n比较小的时候结果误差较大,可以直接暴力计算。
代码
#include<cstdio>
#include<cmath>
using namespace std;
typedef long long ll;
const double pi=acos(-1),e=exp(1);
ll n,k;
double ans;
int main(){
while(~scanf("%lld%lld",&n,&k)){
if(n<=10000){
ans=1;
for(int i=1;i<=n;i++){
ans=ans+log(i)/log(k);
}
printf("%lld\n",(ll)ans);
}else{
ans=0.5*log(2*pi*n)/log(k)+n*log(n/e)/log(k)+1;
printf("%lld\n",(ll)ans);
}
}
return 0;
}