AcWing 315. 旅行 最长公共子序列变形(dfs+线性dp)

这道题是最长公共子序列的变形,难点在输出解决方案上:

表示a序列以第i个元素结尾,b序列以第j个元素结尾的最长公共子序列长度.

我自己的想法是,可以利用acwing314题的思路去输出输出方案.首先,如果s1[i] = s2[j], 那么我们就找到他们前面第一个公共的字母,如果这两个字母对应的位置a,b满足f[a, b] + 1 =  f[i, j],那么这两个字母就可以将答案继承下来,以此类推,我们得到所有位置的答案后排序并离散化即可,但是这样时间复杂度过高了.

接下来,我们用递归输出解决方案,详见代码(y总代码):

#include <bits/stdc++.h>
// #define LOCAL
#define INF 0x3f3f3f3f3f3f3f3f
#define IOS ios::sync_with_stdio(false), cin.tie(0)
// #define int long long
#define debug(a) cerr << #a << "=" << a << endl;
using namespace std;
const int N = 105;
char s1[N], s2[N], f[N][N], path[N];
int n, m;

void dfs(int i, int j, int u, int len){
    if(u > len){
        puts(path + 1);
        return;
    }

    if(s1[i] == s2[j]){
        path[u] = s1[i];
        dfs(i + 1, j + 1, u + 1, len);
    }
    //枚举每个字母在两个序列中下一个出现的位置,我们假设为a,b,f[a][b]==f[i][j],说明这两个是当前情况下的最优解
    else{
        for (int k = 0; k < 26; ++k){
            int a = 0;
            int b = 0;
            for (int x = i; x <= n; ++x)
                if(s1[x] == k + 'a'){
                    a = x;
                    break;
                }
            for (int x = j; x <= m; ++x)
                if(s2[x] == k + 'a'){
                    b = x;
                    break;
                }
            if(a && b && f[a][b] == f[i][j])
                dfs(a, b, u, len);
        }
    }
}
signed main(){
#ifdef LOCAL
    freopen("input.in", "r", stdin);
    freopen("output.out", "w", stdout);
#endif
    IOS;
    cin >> s1 + 1 >> s2 + 1;
    n = strlen(s1 + 1), m = strlen(s2 + 1);
    for (int i = n; i; --i)
        for (int j = m; j; --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]);
        }

    dfs(1, 1, 1, f[1][1]);

}