BZOJ 2048 数学(调和级数) 解题报告
来源:互联网 发布:js比较字符串大小 编辑:程序博客网 时间:2024/05/18 01:47
2048: [2009国家集训队]书堆
Description
Input
第一行正整数 N M
Output
一行(有换行符),L,表示水平延伸最远的整数距离 (不大于答案的最大整数)
Sample Input
Input: 1 100
Output: 49
Input: 2 100
Output: 74
Sample Output
N <= 10^18
数据保证答案 < 10^6
【解题报告】
初中物理难度,很明显,
这一层之上所有书的左端点到它的距离和这一层上所有书的右端点到它的距离是相等的。
这样列方程解出前4层的解.为1/2,1/4,1/6,1/8;显然规律为1/(2*i)。
答案就是和。
这样我们可以用到一个叫做调和级数 的东西,
这是Euler(欧拉)在1734年,利用Newton的成果,首先获得了调和级数有限多项和的值。结果是:
1+1/2+1/3+1/4+...+1/n= ln(n+1)+r(r为常量)
他的证明是这样的:
根据Newton的幂级数有:
ln(1+1/x) = 1/x - 1/2x^2 + 1/3x^3 - ..
于是:
1/x = ln((x+1)/x) + 1/2x^2 - 1/3x^3 + ...
代入x=1,2,…,n,就给出:
给出:1/1 = ln(2) + 1/2 - 1/3 + 1/4 -1/5 + ...1/2 = ln(3/2) + 1/2*4 - 1/3*8 + 1/4*16 - .........1/n = ln((n+1)/n) + 1/2n^2 - 1/3n^3 +...
相加,就得到:
得到:1+1/2+1/3+1/4+...1/n = ln(n+1) + 1/2*(1+1/4+1/9+...+1/n^2) - 1/3*(1+1/8+1/27+...+1/n^3) + ...
后面那一串和都是收敛的,我们可以定义
1+1/2+1/3+1/4+...1/n = ln(n+1) + r
Euler近似地计算了r的值,约为约为0.57721566。这个数字就是后来称作的欧拉常数。
考虑到在n比较小是精度误差较大,我们可以直接暴力算出、
代码如下:
#include<cstdio>#include<cmath>#define r 0.5772156649#define eps 1e-8long long n,m;double ans=0.0;int main(){ scanf("%lld%lld",&n,&m); if(n<=1000) for(int i=1;i<=n;++i) ans+=0.5/i; else {ans=log(n+1.0)+r;ans/=2.0;} ans*=m; printf("%d\n",(int)(ans-eps)); return 0; }
阅读全文
0 0
- BZOJ 2048 数学(调和级数) 解题报告
- BZOJ 1011 数学 解题报告
- BZOJ 1008 组合数学 解题报告
- Uva 1319 - Maximum 解题报告(数学)
- Uva 1388 Graveyard 解题报告(数学)
- 数学题目解题报告
- 数学题目解题报告
- 【数学数列】解题报告
- BZOJ 1588 TREAP 解题报告
- bzoj 3223 splay 解题报告
- BZOJ 3223 Splay 解题报告
- BZOJ 3224 TREAP 解题报告
- BZOJ 1820 DP 解题报告
- BZOJ 2144 LCA 解题报告
- BZOJ 4152 博弈论 解题报告
- BZOJ 1566 DP 解题报告
- BZOJ 3119 贪心 解题报告
- bzoj 1036 树链剖分 解题报告
- JAVA虚拟机系列(四)
- 【今日荐文】七大有效的编程习惯助你成为更好的程序员---Bartlomiej Karalus
- Vue2.0+Webpack项目环境构建到发布
- (个人学习笔记)string 与stringBuilder 的区别
- 2017年codeM资格赛 数码 思维
- BZOJ 2048 数学(调和级数) 解题报告
- B
- 如何动态显示li前面的点
- 使用POI导出excel
- Java学习day3
- MySQL 压缩解决方案
- html5中input的新类型(验证只有在点击提交按钮的时候才触发)
- 学习前端第五个月2017/7/13
- linux基础知识--目录结构