使用动态规划求出最长公共子序列的长度之后,当如何打印出这个最长公共子序列呢?
我们采取递归回溯的策略,从在dp数组上以dp[n][n]为起点,讨论上一步是怎样到达这一步的,然后再进行下一层递归;当到达递归边界条件时,递归中止,开始回溯,这样很自然的就能得到一个答案。
具体举例就是,假如我现在就讨论dp[n][n]这个位置,从上一步到达dp[n][n]可能的方式有三种:
1. 当a[n]==b[n]时,dp[n][n]的上一步就是dp[n-1][n-1]。
2. 当a[n]!=b[n]时,dp[n][n]的上一步有可能是dp[n-1][n],也有可能是dp[n][n-1]。具体是哪一个,根据dp[n][n]的值和它们二者之中的哪个相等来确定。
递归的边界条件便是当r或c之中有一个到达了dp数组的边界0.
很容易想到,只有在序号1的情况,我们才需要对元素进行输出,因为只有在序号1的情况下,a[r]或b[c]才能算作是公共子序列的一部分。
代码如下:
#include <iostream>
#include <algorithm>
using namespace std;
int a[101], b[101], dp[101][101], n, m;
//输出最长公共子序列
void print(int r, int c) {
if (r == 0 || c == 0) {
return;
}
if (a[r] == b[c]) {
//仅有此情况才需要输出
print(r - 1, c - 1);
cout << a[r] << ' ';
}
else if (dp[r][c] == dp[r - 1][c]) {
print(r - 1, c);
}
else if (dp[r][c] == dp[r][c - 1]) {
print(r, c - 1);
}
}
int main()
{
cin >> n;
for (int i = 1; i <= n; i++)
cin >> a[i];
for (int i = 1; i <= n; i++)
cin >> b[i];
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= n; j++) {
if (a[i] == b[j]) {
dp[i][j] = dp[i - 1][j - 1] + 1;
}
else {
dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
}
}
}
cout << dp[n][n] << '\n';
print(n, n);
return 0;
}
这个方法具有一定的局限性,如果想要打印出所有的最长公共子序列,那么这个方法是不work的。