六、动态规划 -- 最长公共子序列

题目描述

给定两个长度分别为 N 和 M 的字符串 A 和 B，求既是 A 的子序列又是 B 的子序列的字符串长度最长是多少。

输入格式

第一行包含两个整数 N 和 M。

第二行包含一个长度为 N 的字符串，表示字符串 A。

第三行包含一个长度为 M 的字符串，表示字符串 B。

字符串均由小写字母构成。

输出格式

输出一个整数，表示最大长度。

数据范围

1≤N,M≤1000

样例

输入样例

4 5
acbd
abedc

输出样例

3

样例解释

在以上的输入中，最长公共子序列为 a b d

---

思路–动态规划

状态表示：

• 集合：f[i][j]表示在a的前i个字母，和b的前j个字母的最长公共子序列长度
• 属性：Max

状态计算：

​ 划分依据：以a[i],b[j]是否包含在子序列当中为依据，因此可以分成四类：

1. a[i]不在，b[j]不在
   f[i][j] = max(f[i][j], f[i - 1][j - 1])

2. a[i]不在，b[j]在
   看似是max=f[i−1][j]max=f[i−1][j] , 实际上无法用f[i−1][j]f[i−1][j]表示，因为f[i−1][j]f[i−1][j]表示的是在a的前i-1个字母中出现，并且在b的前j个字母中出现,此时b[j]不一定出现，这与条件不完全相等，条件给定是a[i]一定不在子序列中，b[j]一定在子序列当中，但仍可以用f[i−1][j]f[i−1][j]来表示，原因就在于条件给定的情况被包含在f[i−1][j]f[i−1][j]中，即条件的情况是f[i−1][j]f[i−1][j]的子集，而求的是max，所以对结果不影响。

3. a[i]在，b[j]不在 同上

4. a[i] 在，b[j]在
   f[i][j] = max(f[i][j], f[i - 1][j - 1] + 1)

实际上，在计算时，①包含在②和③的情况中，所以①不用考虑

---

---

C++代码

#include <iostream>
#include <algorithm>

using namespace std;

const int N = 1010;

int n, m;
char a[N], b[N];
int f[N][N];

int main()
{
    scanf("%d%d", &n, &m);
    scanf("%s%s", a + 1, b + 1);

    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ )
        {
            f[i][j] = max(f[i - 1][j], f[i][j - 1]);
            if (a[i] == b[j]) f[i][j] = max(f[i][j], f[i - 1][j - 1] + 1);
        }

    printf("%d\n", f[n][m]);

    return 0;
}

参考链接

最长公共子序列