博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
【bzoj3000】Big Number 【斯特林公式】
阅读量:5104 次
发布时间:2019-06-13

本文共 1059 字,大约阅读时间需要 3 分钟。

题意:给你两个整数N和K,要求你输出N!的K进制的位数。

题解:首先补了一下对数的运算法则。
1.loga(mn)=loga(m)+loga(n)loga(mn)=loga(m)+loga(n)
2.loga(mn)=loga(m)loga(n)loga(mn)=loga(m)−loga(n)
3.loga(mn)=nloga(m)loga(mn)=nloga(m)
4.loga(mn)=1nloga(m)loga(mn)=1nloga(m)
再有一个n!n!的近似值公式,就是斯特林公式。
n!2πn(ne)nn!≈2πn(ne)n
于是我们就可以求答案了。
ans=logk(n!)+1ans=logk(n!)+1
=>ans=logk(2πn(ne)n)+1ans=logk(2πn(ne)n)+1
=>ans=logk(2πn)+logk((ne)n)+1ans=logk(2πn)+logk((ne)n)+1
=>ans=12logk(2πn)+nlogk(ne)+1ans=12logk(2πn)+nlogk(ne)+1
注意当n比较小的时候结果误差较大,可以直接暴力计算。
代码

#include
#include
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;}

转载于:https://www.cnblogs.com/2016gdgzoi471/p/9476874.html

你可能感兴趣的文章
SQL Server获取月度列表
查看>>
python常用函数
查看>>
python 描点画圆
查看>>
FastDFS使用
查看>>
服务器解析请求的基本原理
查看>>
pycharm 如何设置方法调用字体颜色
查看>>
VUE源码解析心得
查看>>
[HDU3683 Gomoku]
查看>>
【工具相关】iOS-Reveal的使用
查看>>
整体二分——[Poi2011]Meteors
查看>>
数据库3
查看>>
delphi之事件
查看>>
windows server 2008 r2 安装
查看>>
Enigma –> Sadness
查看>>
存储分类
查看>>
下一代操作系统与软件
查看>>
【iOS越狱开发】如何将应用打包成.ipa文件
查看>>
[NOIP2013提高组] CODEVS 3287 火车运输(MST+LCA)
查看>>
Hat’s Words (分成两个字符串考虑)
查看>>
Yii2 Lesson - 03 Forms in Yii
查看>>