7-2 最长公共子串 (20 分)动态规划详解PTA

给定2个字符串，试求出这2个字符串的最长公共子串的长度。

输入格式

输入共2行，每行一个字符串。字符均为小写英文字母。

输出格式

仅一行，包含一个正整数，表示2个字符串的最长公共子串长度。

输入样例

在这里给出一组输入。例如：

ababc
cbaab

输出样例

在这里给出相应的输出。例如：

2

数据范围与提示:

对于30%的数据，保证字符串长度不超过10；
对于60%的数据，保证字符串长度不超过100；
对于90%的数据，保证字符串长度不超过1000；
对于100%的数据，保证字符串长度不超过5000

解题思路

最近学了动态规划，所以打算使用动态规划的方式来求解：首先看问题是否具有最优子结构性质，就是看这个问题的最优解是否包含了子问题的最优解。好比评选成绩最好的学生，先是班级评选，然后班级最好的推到校级评选，那么校级成绩最好的同学一定是产生在班级推选出来的同学之中。对与现在这个查找公共子串问题，我们可以想象暴力解法就是一 一比较，如果两个字母相同，那么子串长度+1。我们可以反过来思考，如果两个子串比较的子串相同，那么它当前的最长子串就一定等于之前最长子串长度+1。比如AABCD和DECF比较，比较当C相等时子串长度 = Max{AAB和CD子串长度} = 0，即长度为0+1 = 1.
现在我们尝试使用数学符号来描述上面的这一规律，这里我们假设第一组输入的字符串为A子串，第二次输入的字符串为B子串那么

接下来我们来看如何把他进一步符号化，我们可以绘制如下表格，当比对的字符串相等时，记为1否则为0：

这个表格只是粗略的记录了两个字符串的一 一对比结果，现在我们需要用表格中的内容来带代表扫描截至当前位置的最长公共子串长度，举个例子来说，我现在用b字符串的中出现的第一个“ a " 去和a字符串进行一 一比较，我们发现a字符串的第一个字符就是” a “,因为这是a字符串的第一个字符，所以前面不可能有子串，该处记为1，再看a字符串的第二个” a “，这时我们可以发现，” a “ 之前的b也和b字符串匹配上了，所以该处的最长公共子串长度就是1 + 1 = 2.通过比较我们能不难发现，其实就是当两个字符相等时，填入该位置的数 = 表格的左斜上位置的值 + 1，如果左斜上超出表格，直接填入1就行 . 类似可以得到下表，箭头即为更新处：
我们也可以进一步更新公式，假设Dp数组二维存放了一开始的表格：

在具体实现的时候其实不需要建立二维表格，建立一个一维表，每次更新是从右到左就行，因为更新到当前行数只和数组左斜上的值有关，这样空间复杂度也会更低，具体代码如下：

#include <iostream>
#include <cmath>
#include <iomanip>
#include <string.h>
using namespace std;

int main()
{
    string a, b;
    cin >> a >> b;
    int ans[a.size()]{}; //记录用的表格
    int max{}; //记录表格中的最大值
    for (int i = 0; i < b.size(); i++)
    {
        for (int j = a.size() - 1; j >= 0; j--)
        {
            if (b[i] == a[j])
            {
                if (j != 0 && i != 0) //防止数组越界
                {
                    ans[j] = ans[j - 1] + 1;
                    if (ans[j] > max)
                        max = ans[j];
                }
                else
                    ans[j] = 1;       //特殊情况直接置一就行
            }
            else
                ans[j] = 0;     //注意要置零
        }
    }
    cout << max;
    return {};
}