最长公共子上升序列

NOI / 2.6基本算法之动态规划

2000:最长公共子上升序列

总时间限制: 10000ms 内存限制: 65536kB
描述
给定两个整数序列，写一个程序求它们的最长上升公共子序列。
当以下条件满足的时候，我们将长度为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

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;

#define inf 0x3f3f3f3f

struct Node {
	int len = 0;//len
	vector<int>v;//序列 
};

int main() {
	int a[501], b[501];
	Node dp[501];
	int m, n;
	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 Max;//b[1..i]和a[1..n]的公共上升子序列长度 
		for (int j = 1; j <= m; j++) {
			//每轮的dp[j]表示b[1..i]和a[1..j]的公共上升子序列
			if (b[i] > a[j] && dp[j].len > Max.len)
				Max = dp[j]; //Max标记b[1..i]和a[1..j]的公共上升子序列长度最长的dp[j] 
			if (b[i] == a[j]) {
				//在之前的Max基础上更新 
				dp[j].len= Max.len + 1;
				dp[j].v = Max.v;
				dp[j].v.push_back(b[i]);
			}
		}
	}
	Node Max = dp[1];
	for (int i = 2; i <= m; i++) {
		if (dp[i].len > Max.len)
			Max = dp[i];
	}
	cout << Max.len << endl;
	for (int i = 0; i < Max.v.size(); i++)
		cout << Max.v[i] << " ";
	cout << endl;
	return 0;
}