BZOJ 2048 2009国家集训队 书堆 数学算法

题目大意：经典的物理上的桌边堆书问题，初中物理老师曾经还讲过，不过只记住了结论。。。没关系，简单证明一下就好

首先我们设由上至下第i本书比它下面那本书多伸出去的长度为a[i]，前缀和为s[i]，那么我们要求的就是s[n]

为了简化问题我们设一本书的长度为1

假设n=1

a[1]=1/2，毫无疑义

然后考虑两本书

两本书的时候，重心明显在距下面那本书左端点的3/4处，故a[2]=1-3/4=1/4

好的我知道了！第一本是1/2，第二本是1/4，那么第三本就是1/8！

这样想的同学都只过了样例。

我们考虑三本书的情况

这个是怎么算的呢？

首先先算三本书的重心距第三本书左端点的距离

第一本书的重心为1/4 + 1/2 + 1/2

第二本书的重心为1/4           + 1/2

第三本书的重心为                    1/2

于是我们可以得到三本书的重心位置为(1/4*2+1/2+1/2*3)/3=1/2*2/3+1/2

于是a[3]=1-(1/2*2/3+1/2)=1/2-1/2*2/3=1/2*(1-2/3)=1/2*1/3=1/6

然后四本书的情况就不用多说了吧

讨论四本书的重心距第四本书左端点的距离

第一本书的重心位置为：1/6 + 1/4 + 1/2 + 1/2

第二本书的重心位置为：1/6 + 1/4           + 1/2

第三本书的重心位置为：1/6                     + 1/2

第四本书的重心位置为：                              1/2

于是四本书重心位置为(1/6*3+1/4*2+1/2+1/2*4)/4=1/2*3/4+1/2

a[4]=1-(1/2*3/4+1/2)=1/2-1/2*3/4=1/2*(1-3/4)=1/2*1/4=1/8

以此类推，我们有a[i]=1/2i

那么这个东西怎么求和呢？n<=10^18，O(n)肯定是不现实的

我们考虑调和级数极限公式

lim(n->+∞)1/1+1/2+1/3+...+1/n = ln(n+1)+r

其中r为欧拉常数 近似值约为0.57721566490153286060651209

但这是极限公式 对于n比较小的情况误差会比较大 所以当n比较小的时候O(n)暴力出解 n比较大的时候套用公式

分界线理论上O(n)能过去就行 1000W左右就可以 但是这题精度实在太低 所以1W就能过去 0MS出解

然后就水过去了。。。一不小心手贱压了下代码 然后就RANK1了。。。

此外100%达成 补番去了

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define eps 1e-10
#define r 0.5772156649
using namespace std;
typedef long long ll;
ll n,m;
double ans;
int main()
{
	int i;
	cin>>n>>m;
	if(n<=10000)
		for(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) );
}