题目描述
伐木工人米尔科需要砍倒M米长的木材。这是一个对米尔科来说很容易的工作,因为他有一个漂亮的新伐木机,可以像野火一样砍倒森林。不过,米尔科只被允许砍倒单行树木。
米尔科的伐木机工作过程如下:米尔科设置一个高度参数H(米),伐木机升起一个巨大的锯片到高度H,并锯掉所有的树比H高的部分(当然,树木不高于H米的部分保持不变)。米尔科就行到树木被锯下的部分。
例如,如果一行树的高度分别为20,15,10和17,米尔科把锯片升到15米的高度,切割后树木剩下的高度将是15,15,10和15,而米尔科将从第1棵树得到5米,从第4棵树得到2米,共得到7米木材。
米尔科非常关注生态保护,所以他不会砍掉过多的木材。这正是他为什么尽可能高地设定伐木机锯片的原因。帮助米尔科找到伐木机锯片的最大的整数高度H,使得他能得到木材至少为M米。换句话说,如果再升高1米,则他将得不到M米木材。
输入输出格式
输入格式:
第1行:2个整数N和M,N表示树木的数量(1<=N<=1000000),M表示需要的木材总长度(1<=M<=2000000000)
第2行:N个整数表示每棵树的高度,值均不超过1000000000。所有木材长度之和大于M,因此必有解。
输出格式:
第1行:1个整数,表示砍树的最高高度。
输入输出样例
输入样例#1:
5 20 4 42 40 26 46
输出样例#1:
36
这题是二分答案题。它满足二分答案题的几个特征:
(1)求最大/最小值;
(2)答案离散(答案有有限种可能);
(3)容易判断答案是否正确(事实上,这题如果要写SPJ,只需要满足ans满足条件但ans+1不满足条件即可)。
二分答案题的做法即是:
(1)确定答案区间;
(2)在保证答案在区间内的前提下,逐步缩小区间;
(3)当区间缩小到仅包含一个可能解时,该可能解即为答案。
这里主要要注意第(2)点,就是如何缩小区间。这是二分的精华所在,也是二分搜索程序如此难写的原因(可以自行查一下历史,或者读一下《编程珠玑》)。
下面贴一下代码。注意数据类型要用unsigned long(为了节省空间,用了unsigned;程序中用usl代表unsigned long)。
#include <iostream>
using namespace std;
typedef unsigned long usl; // 本程序主要使用unsigned long
bool check(usl * a, usl n, usl m, usl pos); //检查解是否可行的函数
int main(void){
usl n,m;
usl l=0,r=1,mid;
cin >> n >> m;
usl a[n];
for (int i=0;i<n;i++){
cin >> a[i];
if (a[i]>r)
r=a[i]; //r保存目前读到的最高的树的高度
}
r--; //在最高的树的高度砍肯定不行啦(这样什么都没砍到),因此-1
//答案区间为:ans in [l,r];即可能的砍树高度从0(全部砍掉)到最高的树的高度-1
while (l<r){
//注意上面这个循环的条件;当l==r,也就是ans in [l,r]变成ans in [l,l]时,ans就确定下来了(ans=l)
mid=(l+r+1)>>1; //即mid=(l+r+1)/2;
//上面不能用mid=(l+r)/2! 因为如果r==l+1,而且check(l)为true,那么每次的mid都是l,结果l和r的值都没有改变,造成死循环
//换句话说,mid要偏向r
//对于这个问题,底下有图解
if (check(a,n,m,mid)) //如果mid是可行解
l=mid; //那么最低就砍到mid,不要再低了(再低就破坏森林生态环境)
else //mid不可行
r=mid-1; //那么必须得砍得比mid还要低,所以最高的高度为mid-1
}
cout << l; //此时输出l和r都是一样的
return 0;
}
bool check(usl * a, usl n, usl m, usl pos){
usl cnt=0;
for (int i=0;i<n;i++)
if (a[i]>pos)
cnt+=(a[i]-pos);
return cnt>=m; //这里相当于return (cnt>=m)?true:false;
}