luogu1020导弹拦截LIS模版（lower_bound用法）

某国为了防御敌国的导弹袭击，发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷：虽然它的第一发炮弹能够到达任意的高度，但是以后每一发炮弹都不能高于前一发的高度。某天，雷达捕捉到敌国的导弹来袭。由于该系统还在试用阶段，所以只有一套系统，因此有可能不能拦截所有的导弹。

输入描述 Input Description
输入导弹依次飞来的高度（雷达给出的高度数据是不大于30000的正整数）

输出描述 Output Description
输出这套系统最多能拦截多少导弹，如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

题目思路

1.第一问
f[i]表示以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;
}