2000:最长公共子上升序列——动态规划

描述

给定两个整数序列，写一个程序求它们的最长上升公共子序列。
当以下条件满足的时候，我们将长度为N的序列S1 , S2 , . . . , SN 称为长度为M的序列A1 , A2 , . . . , AM 的上升子序列：

存在 1 <= i1 < i2 < . . . < iN <= M ，使得对所有 1 <= j <=N，均有Sj = Aij，且对于所有的1 <= j < N，均有Sj < Sj+1。
 

输入

每个序列用两行表示，第一行是长度M(1 <= M <= 500)，第二行是该序列的M个整数Ai (-231 <= Ai < 231 )

输出

在第一行，输出两个序列的最长上升公共子序列的长度L。在第二行，输出该子序列。如果有不止一个符合条件的子序列，则输出任何一个即可。

样例输入

5
1 4 2 5 -12
4
-12 1 2 4

样例输出

2
1 4

思路

这是复杂度为O(n^2)的写法。比较难理解。

代码

#include<bits/stdc++.h>
using namespace std;

int a[505],b[505];

struct node{
	vector<int> v;
	int val=0;//一定要记得初始化
}dp[505];

int main(){
	int n,m;
	cin>>m;
	for(int i=1;i<=m;i++){
		cin>>a[i];
	}
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>b[i];
	}
	
	for(int i=1;i<=n;i++){//被匹配数组放外面
		node maxn;
		for(int j=1;j<=m;j++){
			if(b[i]>a[j]&&maxn.val<dp[j].val){  //如果b[i]<a[j]一定不会更新，所以不需要处理
				maxn=dp[j];
			}
			if(b[i]==a[j]){
				dp[j].val=maxn.val+1;//记录最长上升子序列
				dp[j].v=maxn.v;//整个复制
				dp[j].v.push_back(a[j]); 
			}
		}
	}
	node maxn=dp[1];
	for(int i=2;i<=m;i++){
		if(dp[i].val>maxn.val){
			maxn=dp[i];
		}
	}
	cout << maxn.val <<endl;
	for(int i=0;i<maxn.v.size();i++){
		cout << maxn.v[i] << " ";
	}
	return 0;
}