poj1934 Trip

$题目传送门poj1934$

---

题目描述

统计LCS方案并输出

---

题解

LCS很简单，输出方案很难。
很简单的LCS：

    for (int i = 1; i <= l1; i++) {
        for (int j= 1; j <= l2; j++) {
            if (s1[i] == s2[j]) f[i][j] = f[i - 1][j - 1] + 1;
            else f[i][j] = max(f[i - 1][j], f[i][j - 1]);
        }
    }

很难的输出方案：
根据动规方程，我们可以发现，对于答案有贡献的只有前面一个方程。所以我们需要枚举每个字母，判断其对答案有没有贡献。
我们需要预处理出$g[i][j]/h[i][j]$，表示$A/B$串的前$i$个字母中，$j$最后出现的位置。
然后我们需要处理$A$串的前$l_1$位和$B$串的前$l_2$位的$LCS$，记为$(l_1, l_2, l_{LCS})$。我们枚举每个字母，假设当前字母$i$是$LCS$的最后一位，为了让剩下的$LCS$尽可能长，我们需要找到两串中的最后一个当前字母$i$，然后递归处理$(g[l_1][i], h[l_2][i], l_{LCS} - 1)$。

---

标程

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

const int N = 1050;
const int M = 85;
int f[M][M], g[M][M], h[M][M];
int tot; string ans[N];

void dfs(int l1, int l2, int l, string s) {
    if (!l) {
        ans[++tot] = s;
        return;
    }
    if (!l1 || !l2) return;
    for (int i = 1; i <= 26; i++) {
        int t1 = g[l1][i], t2 = h[l2][i];
        if (f[t1][t2] != l) continue;
        dfs(t1 - 1, t2 - 1, l - 1, char(i + 'a' - 1) + s);
    }
}

int main() {
    int l1, l2;
    string s1, s2;
    cin >> s1 >> s2;
    l1 = s1.size();
    l2 = s2.size();
    s1 = ' ' + s1;
    s2 = ' ' + s2;
    for (int i = 1; i <= l1; i++) {
        for (int j= 1; j <= l2; j++) {
            if (s1[i] == s2[j]) f[i][j] = f[i - 1][j - 1] + 1;
            else f[i][j] = max(f[i - 1][j], f[i][j - 1]);
        }
    }
    for (int i = 1; i <= l1; i++) {
        for (int j = 1; j <= 26; j++) {
            if (s1[i] == j + 'a' - 1) g[i][j] = i;
            else g[i][j] = g[i - 1][j];
        }
    }
    for (int i = 1; i <= l2; i++) {
        for (int j = 1; j <= 26; j++) {
            if (s2[i] == j + 'a' - 1) h[i][j] = i;
            else h[i][j] = h[i - 1][j];
        }
    }
    dfs(l1, l2, f[l1][l2], "");
    sort(ans + 1, ans + tot + 1);
    for (int i = 1; i <= tot; i++)
        cout << ans[i] << endl;
    return 0;
}

---