最长上升子序列问题(LIS)

题目

有一个长为n的数列a0,a1,...,a(n-1)。请求出这个序列中最长的上升子序列的长度。上升子序列指的是对于任意的i<j都满足ai<aj的子序列。

限制条件：1<=n<=1000, 0<=ai<=1000000

样例输入：n=5         a={4,2,3,1,5}

输出：3（a1,a2,a4构成的子序列）

 分析

最长上升子序列（LIS，Longest Increasing Subsequence）这是个经典的问题，也代表了一个类型的题目。寻找LIS，很容易想到使用动态规划，那么如何使用呢？

首先定义一个一维数组dp[n]来用于记录对应位置的LIS，然后寻找该问题的边界，即可分解的最小的状态。因为就单个数字而言，其可以被看为是一个长度为1的序列，所以边界应该是dp[0]=1。

接着来寻找状态转移方程。对一个位置的dp[i]而言，从序列开始到该位置的LIS的长度有两种可能，分别是1(ai<aj，其中0<=j<i）和dp[j]+1(ai>aj，其中0<=j<i)。转换成公式如下

完整代码

#include<iostream>
#include<algorithm>
using namespace std;

int main(){
//	读入数据 
	int n;
	cin>>n;
	int a[n],dp[n];
	for(int i=0;i<n;i++)
		cin>>a[i];
	
//	存放当前位置的LIS 
	int Max=0; 
	for(int i=0;i<n;i++){
//		每个位置的初始长度都为1 
		dp[i]=1;
		for(int j=0;j<i;j++){ 
			if(a[j]<a[i]){
				dp[i]=max(dp[i],dp[j]+1);
			}
		}
		Max=max(Max,dp[i]);
	} 
	cout<<Max<<endl;
	
	return 0;
}

同类型问题

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

输入导弹依次飞来的高度（雷达给出的高度数据是不大于30000的正整数），计算这套系统最多能拦截多少导弹。

【输入格式】一行，为导弹依次飞来的高度

【输出格式】最多能拦截的导弹数

【样例输入】389 207 155 300 299 170 158 65

【样例输出】

6

分析

这个问题的本质是找“最长下降子序列”，所以用同样的思路就可以解答。

完整代码

#include<iostream>
#include<algorithm>
using namespace std;

int dp[100];
int a[100],t=1;

int main(){
//	读入数据 
    cin>>a[0];
    while(cin.get()!='\n')
	    cin>>a[t++];

//	存放当前最长下降子序列 
	int Max=0;
	for(int i=0;i<t;i++){
		dp[i]=1;
		for(int j=0;j<i;j++){
			if(a[i]<a[j])	dp[i]=max(dp[i],dp[j]+1);
		}
		Max=max(Max,dp[i]);
	}
	cout<<Max<<endl;	
		
	return 0;
}

 

 参考了《挑战程序设计竞赛》P64-65