动态规划——导弹拦截（最长不上升子序列、最长上升子序列）

题目链接

题目描述

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

输入导弹依次飞来的高度（雷达给出的高度数据是 ≤ 50000 的正整数），计算这套系统最多能拦截多少导弹，如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

输入格式

1 行，若干个整数（个数 ≤ 100000）

输出格式

2 行，每行一个整数，第一个数字表示这套系统最多能拦截多少导弹，第二个数字表示如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

输入输出样例

输入

389 207 155 300 299 170 158 65

输出

6
2

方法一：O(n²)

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
int a[maxn],dp1[maxn],dp2[maxn];

int main()
{
	int n=0;
	while(cin>>a[n++]);
	n--; 
	int ans1=0,ans2=0;
	for(int i=n-1;i>=0;i--){	// 以 i 开头的最长不上升子序列 
		dp1[i]=1;
		for(int j=i+1;j<n;j++){
			if(a[j]<=a[i]){
				dp1[i]=max(dp1[i],dp1[j]+1);
			}
		}
		ans1=max(ans1,dp1[i]);
	}
	for(int i=0;i<n;i++){	// 以 i 结尾的最长上升子序列 
		dp2[i]=1;
		for(int j=0;j<i;j++){
			if(a[j]<a[i]){
				dp2[i]=max(dp2[i],dp2[j]+1);
			}
		}
		ans2=max(ans2,dp2[i]);
	}
	cout<<ans1<<endl;
	cout<<ans2<<endl;
	return 0;
}

方法二：O(nlogn)

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
int a[maxn],dp1[maxn],dp2[maxn];

int main()
{
	int n=0;
	while(cin>>a[n++]);
	n--;
	dp1[0]=a[0];	// 用于求最长不上升子序列 
	dp2[0]=a[0];	// 用于求最长上升子序列 
	int cnt1=0,cnt2=0;
	for(int i=1;i<n;i++){
		if(dp1[cnt1]>=a[i]){
			dp1[++cnt1]=a[i];
		}
		else{
			*upper_bound(dp1,dp1+cnt1,a[i],greater<int>() )=a[i];
		}
		if(dp2[cnt2]<a[i]){
			dp2[++cnt2]=a[i];
		}
		else{
			*lower_bound(dp2,dp2+cnt2,a[i])=a[i];
		}
	}
	cout<<cnt1+1<<endl;
	cout<<cnt2+1<<endl;
	return 0;
}