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;
}