当前位置: 代码迷 >> 综合 >> 51nod 最长公共子序列+输出路径
  详细解决方案

51nod 最长公共子序列+输出路径

热度:57   发布时间:2023-09-20 20:06:18.0

当x = 0 或 y = 0时 f[x][y] = 0

当a[x] = b[y]时  f[x][y] = f[x-1][y-1]+1

当a[x] != b[y]时 f[x][y] = max(f[x][y-1], f[x-1][y])

注意这里字符串要从1开始,因为转移方程里面0表示这个字符串为空的时候。

动态规划涉及到字符串最好从1开始

还有就是路径的问题,可以把二维数组画出来,推出转移的方向。

这个方法我还是第一次见,很牛逼。

#include<cstdio>
#include<cstring>
#include<vector>
#define REP(i, a, b) for(int i = (a); i < (b); i++)
using namespace std;const int MAXN = 1123;
char a[MAXN], b[MAXN];
int f[MAXN][MAXN], pre[MAXN][MAXN];void print(int x, int y)
{if(pre[x][y] == 1){	print(x - 1, y - 1);putchar(a[x]);}else if(pre[x][y] == 2) print(x - 1, y);else if(pre[x][y] == 3) print(x, y - 1);
}int main()
{scanf("%s%s", a + 1, b + 1);int n = strlen(a + 1), m = strlen(b + 1);REP(i, 1, n + 1)REP(j, 1, m + 1){if(a[i] == b[j]){f[i][j] = f[i - 1][j - 1] + 1;pre[i][j] = 1;}else{if(f[i - 1][j] > f[i][j - 1]){f[i][j] = f[i - 1][j];pre[i][j] = 2;}else{f[i][j] = f[i][j - 1];pre[i][j] = 3;}}}print(n, m);puts("");return 0;
}