最长公共上升子序列(输出序列)
内存限制: 256 MiB 时间限制: 1000 ms 标准输入输出 题目类型: 传统 评测方式: 文本比较
题目描述
研究发现,大猩猩的基因序列和人的基因序列只有1.3%的区别,更进一步,不仅仅离人最近的大猩猩和人的基因序列高度近似,就连以打洞为生的老鼠和人的基因序列也有高达95%的相同序列。于是有魔法师提出一个大胆设想,即改变人类的某些特定基因以期产生超级人类。
现在,他们要做的第一步是将两种不同生物的基因序列转换成两个整数序列,并试图确定他们的最大公共上升子序列的长度,例如有 A 序列为 4 3 2 1 7 8 9,B 序列为 7 8 9 4 3 2 1,其最长公共子序列是 4 3 2 1,而最长公共递增子序列应该是 7 8 9。
输入格式
输入每个序列由 M 个整数组成(1 ≤ M ≤ 500),M 个整数范围在(-2^31 ≤ Bi < 2^31)之间。
输出格式
第一行输出最长公共上升子序列长度 L,第二行输出该子序列,如果该序列有多个答案,输出任意序列 1 的编号字典序最小的一个。
样例
样例输入
复制5
1 4 2 5 -12
4
-12 1 2 4
样例输出
复制2
1 4
数据范围与提示
输出按第一个序列的数字的编号字典需最小的那一组最长公共子序列
#include <bits/stdc++.h>
using namespace std;
const int N = 2e3 + 5;
int e, n, m, ans;
int a[N], b[N], pre[N];
int dp[N][N];
stack<int> s;
int main() {
scanf("%d", &m);
for(int i = 1; i <= m; i++) scanf("%d", &b[i]);
scanf("%d", &n);
for(int i = 1; i <= n; i++) scanf("%d", &a[i]);
for(int i = 1; i <= n; i++) {
for(int j = 1; j <= m; j++) {
if(a[i] == b[j]) {
dp[i][j] = 1;
for(int k = 1; k <= j; k++) {
if(b[k] < b[j] && dp[i - 1][k] + 1 > dp[i][j]) {
dp[i][j] = dp[i - 1][k] + 1, pre[j] = k;
}
}
} else dp[i][j] = dp[i - 1][j];
}
}
for(int i = 1; i <= m; i++) {
if(dp[n][i] > ans) ans = dp[n][i], e = i;
}
printf("%d\n", ans);
while(e != 0) {
s.push(b[e]);
e = pre[e];
}
while(!s.empty()) {
printf("%d ", s.top());
s.pop();
}
return 0;
}