【题解】最长公共子序列（DP）

最新推荐文章于 2022-11-15 10:57:58 发布

_BOSS_

最新推荐文章于 2022-11-15 10:57:58 发布

阅读量2.8k

点赞数 1

分类专栏：题解动态规划（DP）

本文链接：https://blog.csdn.net/STRVE/article/details/104223727

版权

题解同时被 2 个专栏收录

50 篇文章 1 订阅

订阅专栏

动态规划（DP）

12 篇文章 1 订阅

订阅专栏

题面

【题目描述】
一个给定序列的子序列是在该序列中删去若干元素后得到的序列。确切地说，若给定序列 $X =$ { $x_1, x_2,…, x_m$ }，则另一序列 $Z =$ { $z_1, z_2,…, z_k$ }是 $X$ 的子序列是指存在一个严格递增的下标序列 { $i_1, i_2,…, i_k$ }，使得对于所有 $j = 1, 2, \dots, k$ 有 $X_{i_j}=Z_j$ 。
例如，序列 $Z =$ { $B, C, D, B$ }是序列 $X =$ { $A, B, C, B, D, A, B$ }的子序列，相应的递增下标序列为{ $2, 3, 5, 7$ }。给定两个序列 $X$ 和 $Y$ ，当另一序列 $Z$ 既是 $X$ 的子序列又是 $Y$ 的子序列时，称 $Z$ 是序列 $X$ 和 $Y$ 的公共子序列。例如，若$X= ${ $A, B, C, B, D, A, B$ }和 $Y =$ { $B, D, C, A, B, A$ }，则序列{ $B, C, A$ }是 $X$ 和 $Y$ 的一个公共子序列，序列{ $B, C, B, A$ }也是 $X$ 和 $Y$ 的一个公共子序列。而且，后者是 $X$ 和 $Y$ 的一个最长公共子序列，因为 $X$ 和 $Y$ 没有长度大于 $4$ 的公共子序列。给定两个序列 $X =$ { $x_1, x_2, …, x_m$ }和 $Y =$ { $y_1, y_2, … , y_n$ }，要求找出 $X$ 和 $Y$ 的一个最长公共子序列。
【输入】
输入共两行。每行为一个由大写字母构成的长度不超过 $200$ 的字符串，表示序列 $X$ 和 $Y$ 。
【输出】
输出第一行为一个非负整数。表示所求得的最长公共子序列的长度。若不存在公共子序列，则输出仅有一行输出一个整数 $0$ 。
【样例输入】

ABCBDAB
BDCABA

【样例输出】

算法分析

状态：
$d p [i] [j]$ ：字符串 $S$ 前 $i$ 个字符与字符串 $T$ 前 $j$ 个字符的最长公共子序列的长度。

状态转移方程：
$（ 1 ）$ S[i]$不在公共子序列中，那么字符串 $S$ 前 $i$ 个字符与字符串 $T$ 前 $j$ 个字符的最长公共子序列的长度 $=$ 字符串 $S$ 前 $i - 1$ 个字符与字符串 $T$ 前 $j$ 个字符的最长公共子序列的长度：
$d p [i] [j] = d p [i - 1] [j]$

$（ 2 ）$ $T [j]$ 不在公共子序列中，那么字符串 $S$ 前 $i$ 个字符与字符串 $T$ 前 $j$ 个字符的最长公共子序列的长度 $=$ 字符串 $S$ 前 $i$ 个字符与字符串 $T$ 前 $j - 1$ 个字符的最长公共子序列的长度：
$d p [i] [j] = d p [i] [j - 1]$

$（ 3 ）$ $S [i] = T [j]$ ，将 $S [i] 与 T [j]$ 加入到公共子序列中，那么字符串 $S$ 前 $i$ 个字符与字符串 $T$ 前 $j$ 个字符的最长公共子序列的长度 $=$ 字符串 $S$ 前 $i$ 个字符与字符串 $T$ 前 $j - 1$ 个字符的最长公共子序列的长度 $+ 1$ ：
$d p [i] [j] = d p [i - 1] [j - 1] + 1$
最后取一个最大值。

边界条件
$d p [i] [0] = 0$ ，字符串 $S$ 前 $i$ 个字符与字符串 $T$ 前 $0$ 个字符的最长公共子序列的长度为0。
$d p [0] [j] = 0$ ，字符串 $S$ 前 $0$ 个字符与字符串 $T$ 前 $j$ 个字符的最长公共子序列的长度为0。
定义时，默认为0，就不需要初始化。

参考程序

有两种写法：
$1 .$ 字符串下标0开始，为了避免i-1，j-1越界，需要分类讨论：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int dp[1100][1100],ans;
char a[1100],b[1100]; 
int main()
{
    gets(a);
    gets(b);
    int lena=strlen(a);
    int lenb=strlen(b);
    for(int i=0;i<lena;i++)
        for(int j=0;j<lenb;j++)
        {
            if(i==0&&j==0)
            {
                if(a[0]==b[0]) dp[0][0]=1;
                else dp[0][0]=0;
            }
            else if(i==0)
            {
                if(a[0]==b[j]) dp[0][j]=1;
                else dp[0][j]=max(dp[0][j],dp[0][j-1]); 
            } 
            else if(j==0)
            {
                if(a[i]==b[0]) dp[i][0]=1;
                else dp[i][0]=max(dp[i][0],dp[i-1][0]); 
            }
            else
            {
                dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
                if(a[i]==b[j]) dp[i][j]=max(dp[i][j],dp[i-1][j-1]+1);
            }
        }
    cout<<dp[lena-1][lenb-1]<<endl;
    return 0;
}

$2 .$ 字符串下标1开始：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int dp[1100][1100],ans;
char a1[1100],b1[1100],a[1100],b[1100]; 
int main()
{
    gets(a1);
    gets(b1);
    int lena=strlen(a1);
    int lenb=strlen(b1);
    for(int i=1;i<=lena;i++)
        a[i]=a1[i-1];
    for(int i=1;i<=lenb;i++)
        b[i]=b1[i-1];
    for(int i=1;i<=lena;i++)
        for(int j=1;j<=lenb;j++)
        {
            dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
            if(a[i]==b[j]) dp[i][j]=max(dp[i][j],dp[i-1][j-1]+1);
        }
    cout<<dp[lena][lenb]-1<<endl;
    return 0;
}