这是一道水题, 但是对我的二分有很大启发, 也可以 n2DP 但是显然单调栈 nlogn 更佳
题目描述
某国为了防御敌国的导弹袭击, 发展出一种导弹拦截系统. 但是这种导弹拦截系统有一个缺陷: 虽然它的第一发炮弹能够到达任意的高度, 但是以后每一发炮弹都不能高于前一发的高度. 某天, 雷达捕捉到敌国的导弹来袭. 由于该系统还在试用阶段, 所以只有一套系统, 因此有可能不能拦截所有的导弹.
输入导弹依次飞来的高度(雷达给出的高度数据是 \ le 50000≤50000 的正整数), 计算这套系统最多能拦截多少导弹, 如果要拦截所有导弹最少要配备多少套这种导弹拦截系统.
输入输出格式
输入格式:
11 行, 若干个整数(个数 \ le 100000≤100000)
输出格式:
22 行, 每行一个整数, 第一个数字表示这套系统最多能拦截多少导弹, 第二个数字表示如果要拦截所有导弹最少要配备多少套这种导弹拦截系统.
输入输出样例
输入样例 #1: 复制
389 207 155 300 299 170 158 65
输出样例 #1: 复制 6 2
说明
为了让大家更好地测试 n 方算法, 本题开启 spj,n 方 100 分, nlogn200 分
每点两问, 按问给分
这道题就是求最长不上升和上升子序列;
这是为什么呢, 第一问很好理解, 第二问为为什么求最长上升就行了呢??
这个可以脑补一下, 相当于是一个个的上升序列, 然后他们的队头构成了最后的最长上升
然后我就想, 为什么不能在每一个队里选出来更多的数组成更长的答案?
这是不成立的, 因为这个队列是单调递减, 只能一个队里选一个(因为如果选两个就不是上升了), 所以一个队里选一个, 最后答案还是一样的, 并不影响
然后呢, n2 是白给的(思路很好懂, 暴力贪心)
- #include<cstdio>
- #include<algorithm>
- #include<cstring>
- using namespace std;
- int cnt;
- int a[100011];
- int dp[100011];
- int main()
- {
- // freopen("testdata.in","r",stdin);
- while(scanf("%d",&a[++cnt])!=EOF)continue;cnt--;
- dp[1]=1;
- for(int i=2;i<=cnt;i++)
- {
- int tmp=0;
- for(int j=1;j<=i-1;j++)
- {
- if(a[j]>=a[i])tmp=max(tmp,dp[j]);
- }
- dp[i]=max(dp[i],tmp+1);
- }
- int ans=-1;
- for(int i=1;i<=cnt;i++)ans=max(ans,dp[i]);
- printf("%d\n",ans);
- memset(dp,0,sizeof(dp));
- dp[1]=1;
- for(int i=2;i<=cnt;i++)
- {
- int tmp=0;
- for(int j=1;j<=i-1;j++)
- {
- if(a[j]<a[i])tmp=max(tmp,dp[j]);
- }
- dp[i]=max(dp[i],tmp+1);
- }
- ans=-1;
- for(int i=1;i<=cnt;i++)ans=max(ans,dp[i]);
- printf("%d\n",ans);
- return 0;
- }
- /*1 2 3 4 5*/
- /*
- 389 207 155 300 299 170 158 65
- */
然后便是单调栈的写法(一开始以为是单调队列, 但是单调队列是维护需要出队的, 用来缓存数组, 优化时间复杂度, 显然这道题并不用)
单调栈的思路就是维护一个单调的栈, 每次寻找最底层比它大 (或小) 的位置将当前元素插入, 最后栈内有多少个元素就是答案
那么插入的过程呢, 对, 用二分来维护, 这样才能保证是每次 log
那么我们来看一下这个二分
- for(int i=2;i<=cnt;i++)
- {
- int l=1,r=top,mid,tt=top+1;
- while(l<=r)
- {
- mid=(l+r)/2;
- if(stk[mid]>=a[i])l=mid+1;
- else r=mid-1,tt=mid;
- }
- if(tt==top+1)top++;
- stk[tt]=a[i];
- }
- for(int i=2;i<=cnt;i++)
- {
- int l=1,r=top,mid,tt=top+1;
- while(l<=r)
- {
- mid=(l+r)/2;
- if(stk[mid]<a[i])l=mid+1;
- else r=mid-1,tt=mid;
- }
- if(tt==top+1)top++;
- stk[tt]=a[i];
- }
这个 tt=mid, 也就是记录 ans 的地方到底应该加在哪里呢, 之前是背模板, 以为就是在 r 那里, 但是不是这样的,
应该是哪个符合你想要的在哪里记录 ans, 也就是有一个满足就记录一下
然后就没什么问题了, 注意对 top 指针的维护就好;
- #include<cstdio>
- #include<algorithm>
- #include<cstring>
- using namespace std;
- int cnt,head,tail,ans=-1,top=1;
- int a[100011];
- int stk[100011];
- int main()
- {
- while(scanf("%d",&a[++cnt])!=EOF)continue;cnt--;
- stk[1]=a[1];
- for(int i=2;i<=cnt;i++)
- {
- int l=1,r=top,mid,tt=top+1;
- while(l<=r)
- {
- mid=(l+r)/2;
- if(stk[mid]>=a[i])l=mid+1;
- else r=mid-1,tt=mid;
- }
- if(tt==top+1)top++;
- stk[tt]=a[i];
- }
- printf("%d\n",top);
- memset(stk,0,sizeof(stk));
- stk[1]=a[1],top=1;
- for(int i=2;i<=cnt;i++)
- {
- int l=1,r=top,mid,tt=top+1;
- while(l<=r)
- {
- mid=(l+r)/2;
- if(stk[mid]<a[i])l=mid+1;
- else r=mid-1,tt=mid;
- }
- if(tt==top+1)top++;
- stk[tt]=a[i];
- }
- printf("%d\n",top);
- return 0;
- }
By 浅夜_MISAKI
LuoGuP1020 导弹拦截[单调栈 + 二分深刻理解] [做题报告]
来源: http://www.bubuko.com/infodetail-2874295.html