某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度。某天,雷达捕捉到敌国的导弹来袭。由于该系统还在试用阶段,所以只有一套系统,因此有可能不能拦截所有的导弹。
输入描述 Input Description
输入导弹依次飞来的高度(雷达给出的高度数据是不大于30000的正整数)
输出描述 Output Description
输出这套系统最多能拦截多少导弹,如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。
题目思路
1.第一问
f[i]表示以h[i]为结尾最长的不增子序列
f
[
i
]
=
m
a
x
(
f
[
j
]
+
1
)
j
<
=
i
,
h
[
j
]
>
=
h
[
i
]
f[i]=max(f[j]+1) j<=i, h[j]>=h[i]
f[i]=max(f[j]+1)j<=i,h[j]>=h[i]
复杂复O(n^2)
2. 第二问
下降子序列的个数等于最长上升子序列的长度
所以题目等价于求最长上升子序列的长度
O(n^2) 求LIS
for(i=1;i<=n;i++)
{
f[i]=1;
for(int j=1;j<i;j++)
if(a[j]<a[i])
f[i]=max(f[i],f[j]+1);
}
ans=0;
for(i=1;i<=n;i++) ans=max(ans,f[i]);
cout<<ans<<endl;
O(nlogn) 求LIS
low[1]=a[1];
l=1;
for(int i=2;i<=n;i++)
{
if(a[i]>low[l]) low[++l]=a[i];
else
{
int pos=lower_bound(low+1,low+1+l,a[i])-low;
low[pos]=a[i];
}
}
写在前面
假设我们查找x,那么:
lower_bound会找出序列中第一个大于等于x的数
upper_bound会找出序列中第一个大于x的数
序列必须是升序
bound的内置比较符是’<’
如果是降序呢?
lower_bound(a + 1, a + 1 + n, x, greater <int> () );
总结
对一个下降序列a
int p = lower_bound(a + 1, a + 1 + n, x, greater () ) - a;
a[p]即a[1]到a[n]中第一个小于等于x的数
https://www.luogu.org/blog/w1049/solution-p1020
题目代码
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int inf=0x3f3f3f3f;
const int N=1e5+100;
int a[N];
int d1[N],d2[N];
int n;
int main()
{
int i=0,x;
while(cin>>x)
{
a[++i]=x;
}
n=i;
int l1=1,l2=1;
d1[1]=a[1];
d2[1]=a[1];
for(int i=2;i<=n;i++)
{
if(a[i]<=d1[l1]) d1[++l1]=a[i];
else
{
int p1=upper_bound(d1+1,d1+1+l1,a[i],greater<int> () )-d1;
d1[p1]=a[i];
}
if(a[i]>d2[l2]) d2[++l2]=a[i];
else
{
int p2=lower_bound(d2+1,d2+1+l2,a[i])-d2;
d2[p2]=a[i];
}
}
cout<<l1<<endl<<l2<<endl;
return 0;
}