第一章 动态规划 最长公共子序列

1. 题目最长公共子序列

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

输入格式
第一行包含两个整数 N 和 M。
第二行包含一个长度为 N 的字符串，表示字符串 A。
第三行包含一个长度为 M 的字符串，表示字符串 B。
字符串均由小写字母构成。

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

数据范围
1≤N,M≤1000
输入样例：
4 5
acbd
abedc
输出样例：
3

2. 分析

我们将第一个序列记为a[i],第二个序列记为b[j]。根据最后一个a[i]和最后一个b[j]进行集合划分。
集合划分讲究个不重不漏，但是不重不是所有属性都需要保证，例如max，不需要各种集合是不重复的，即使重复了，也能选出max。
我们看集合计算。

1. 不含a[i],也不含b[j]
   _ _ _ _ _ _ _ _ _ _ a[i - 1]
   _ _ _ _ _ _ _ _ _ _ b[j - 1]
   按照定义就是f[i - 1,j - 1]
2. 含a[i],含b[j]
   _ _ _ _ _ _ _ _ _ _ a[i]
   _ _ _ _ _ _ _ _ _ _ b[j]
   满足这个情况需要ai == bj,是f[i - 1,j - 1] + 1
3. 含a[i],不含b[j]
   _ _ _ _ _ _ _ _ _ _ a[i]
   _ _ _ _ _ _ _ _ _ _ b[j - 1]
   按照我们直接的想法就是应该记为f[i ,j - 1]，可是这个真的只含a[i],不含b[j]吗？其实按照定义f[i,j - 1]来说，这里面其实也包含了 不含a[i],也不含b[j]的集合，但是就像我们上面讨论的，重复对于求max不是硬性要求，不漏才是关键，所以我们仍用f[i,j - 1]来表示这个集合
4. 不含a[i],含b[j]
   _ _ _ _ _ _ _ _ _ _ a[i - 1]
   _ _ _ _ _ _ _ _ _ _ b[j]
   f[i - 1,j]
   我们得到
   $$f[i,j]=max(f[i-1,j-1],f[i,j-1],f[i-1,j])$$
   $$if(a[i]==b[j])f[i,j]=max(f[i][j],f[i-1,j-1]+1)$$

3. 代码

#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1010;
char a[N],b[N];
int f[N][N];

int main()
{
    int n,m;
    cin >> n >> m;
    for(int i = 1; i <= n; i ++) cin >> a[i];
    for(int i = 1; i <= m; i ++) cin >> b[i];
    
    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);
        }
    }
    cout << f[n][m];
}

参考资料

Acwing基础班