最长公共子序列的输出

 使用动态规划求出最长公共子序列的长度之后，当如何打印出这个最长公共子序列呢？

我们采取递归回溯的策略，从在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的。