推荐一篇博客:详解二分查找算法
这篇博客当中查找一个数,左边界,右边界都做出了解释
而且各种取等,加一减一的细节也解释的很全面,算是保姆级服务了
补充一下为什么左边界和右边界能反馈回来正确信息:
1.左边界二分得到的是最左侧的搜寻数,如果不存在,是小于搜寻数的第一个数
因为左边界二分锁死了右边界,如果搜寻数存在,那么一定始终在l,r中间卡着
会不会出错呢?如果会出错,只可能是l=mid+1导致的,那么反证法,如果这步致错了,那么致错前这步mid位置就应该是搜寻数
如果搜寻数不存在,那么r先前被死锁在搜寻数的左边界,现在搜寻数不存在,那么r被死锁在搜寻数左侧第一个数,l被往右逼,最后搜到结果会与预期一致
并且因为二分的搜寻过程本身是不犯错的,同时l和r在夹,所以答案也是不犯错的,想不通可以像刚才一样反证
2.右边界二分得到的无论存不存在搜寻数,都是大于搜寻数的第一个数
因为左边界被挖掉而不是锁住,其余原理基本同上
以上是根据引用的博文补充的,二分很灵活,也可以自己根据原理写
题目传送门:P1020 [NOIP1999 普及组] 导弹拦截
P1020 [NOIP1999 普及组] 导弹拦截
题目描述
某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度。某天,雷达捕捉到敌国的导弹来袭。由于该系统还在试用阶段,所以只有一套系统,因此有可能不能拦截所有的导弹。
输入导弹依次飞来的高度,计算这套系统最多能拦截多少导弹,如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。
输入格式
一行,若干个整数,中间由空格隔开。
输出格式
两行,每行一个整数,第一个数字表示这套系统最多能拦截多少导弹,第二个数字表示如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。
输入输出样例
输入 #1复制
389 207 155 300 299 170 158 65
输出 #1复制
6 2
说明/提示
对于前 50%50% 数据(NOIP 原题数据),满足导弹的个数不超过 104104 个。该部分数据总分共 100100 分。可使用O(n2) 做法通过。
对于后 50%50% 的数据,满足导弹的个数不超过 105105 个。该部分数据总分也为 100100 分。请使用 O(nlogn) 做法通过。
对于全部数据,满足导弹的高度为正整数,且不超过 5×1045×104。
此外本题开启 spj,每点两问,按问给分。
upd 2022.8.24upd 2022.8.24:新增加一组 Hack 数据。
首先确定性质是动态规划的题,第一问无疑最长不严格下降子序列
第二问比较特别,需要用到Dilworth定理,转化为求最长严格上升子序列
(这百度百科最后一句话说的不对)
那么接下来考虑时间复杂度
想要过这个需要n log n做法,如果dp暴力是n平方显然过不了,要用贪心+二分
直接dp,好处是能够得知LIS序列,缺点是时间复杂度O(n)
贪心+二分,好处是时间复杂度O(n log n),缺点是只知道LIS长度,不知道具体序列
注意看这句话“缺点是只知道LIS长度”,无论dp还是贪心,都要满足无后效性
我们只要考虑怎么从过去的烂摊子转移到现在的情况就可以了,如何转移?
每次转移只依赖于哪一个可能最优的LIS序列的结尾拼接过来,当然如果能拼接就拼接喽,这个过程满足动态规划,不满足贪心
根据上面的想法,完全可以把每一个“可能最优的LIS序列”结尾元素都存下来,存进一个数组里
在这个的基础上,我每次修改数组中元素时,当然是选择场面更好的“结局”更好,比如有一堆结尾300 200 100 50,还有一堆结尾200 150 100 50,那么我绝对选择前者,共有元素200、100、50,但是前者的300强于150,因为“兼容性”更高,也就是每次修改的时候挑选对导弹下降高度影响最小的位置,这是贪心的思想
仔细想想,其实贪心的具体实现就是插入一个end时,找到的第一个“逆序对”,刷新掉它
这样相当于办正事的同时排好序了,那么这个程序可以进一步优化,采用二分
此时问题1属于寻找右边界(但不缩紧左侧)的二分,问题2属于普通二分
当然,不想这么细也可以,弄清楚二分回馈回来的是哪个位置就可以了,问题1是不严格的,需要“逆序对”,问题2是严格的,需要“平替”
解释已经穷尽了,代码如下:
#include<iostream>
using namespace std;
int a[500001];int n=0;
int end1[500001];int cnt_end1=0;
int end2[500001];int cnt_end2=0;
int main()
{
while(cin>>a[++n]){;}//input
n--;
end1[1]=a[1];cnt_end1++;
for(int i=2;i<=n;i++) //lis
{
int l=1,r=cnt_end1+1;
while(l<r){
int mid=(l+r)/2;
if(end1[mid]==a[i]){
l=mid+1;
}
else if(end1[mid]>a[i]){
l=mid+1;
}
else{
r=mid;
}
}
end1[l]=a[i];
cnt_end1=cnt_end1>l? cnt_end1:l;
}
cout<<cnt_end1<<endl;
end2[1]=a[1];cnt_end2++;
for(int i=2;i<=n;i++)//reverse
{
int l=1,r=cnt_end2+1;
while(l<r){
int mid=(l+r)/2;
if(end2[mid]==a[i]){
l=mid;
break;
}
else if(end2[mid]<a[i]){
l=mid+1;
}
else{
r=mid;
}
}
end2[l]=a[i];
cnt_end2=cnt_end2>l? cnt_end2:l;
}
cout<<cnt_end2;
return 0;
}