最长公共子序列问题详解（LCS）

NothingAtall.

已于 2024-01-06 21:01:17 修改

阅读量971

点赞数 1

分类专栏： problemS 文章标签：算法动态规划 c语言 c++ 抽象代数

于 2022-10-04 15:54:41 首次发布

本文链接：https://blog.csdn.net/qq_63786973/article/details/127162882

版权

problemS 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

最长公共子序列问题（LCS）

给定两个字符串 S (n) 和 T (m) ，求出两个字符串的最长公共子序列的长度

$\color{Red}！$ 限制条件

1 <= n， m <= 1000

子序列：

字符串S的子序列，是一组按从左到右顺序出现的字符，但不一定是连续的

例如，我们有一字符串
$\color{gray}ACTTGCG$

ACT ， ATTC，T 这些都是此字符串的子序列
TTA 不是此字符串的子序列

公共子序列

例如，我们有字符串 S1 和 S2
$\color{gray}S1 =AAACCGT GAGT T AT T CGT T CT AGAA$

$\color{gray}S2 = CACCCCT AAGGT ACCT T T GGT T C$

那么对于 S1和 S2来说，它们的 LCS 就是:
$\color{gray}S1 =\color{gray}AA\color{Red}{ACC}\color{gray}G\color{Red}T\color{gray}G\color{Red}AG\color{gray}T\color{Red} T A\color{gray}T T \color{Red}C\color{gray}G\color{Red}T T \color{gray}C\color{Red}T \color{gray}A\color{Red}G\color{gray}AA$

$\color{gray}S2 = C\color{Red}A\color{gray}CC\color{Red}CCT A\color{gray}AG\color{Red}GT AC\color{gray}C\color{Red}T T T G\color{gray}GT T C$

$\color{gray}LCS = \color{Red}ACCTAGTACTTTG$

$\color{Turquoise}算法分析$

我们的目的是要找出字符串 X = < x1，x2，···，xm> 和 Y = < y1 ， y2，···，yn > 最长公共子序列

假设 Z = < z1， z2，···，zk > 是它们的最长公共子序列

若 $x_m = y_n$ ；
- 此时 $z k = x m = y n$ ，我们有 $Z_{k-1}$ 是 $X_{m-1} 与 Y_{n-1 }$ 的 LCS
若 $x_m ! = y_n$
- 则 $z_k ! = x_m$ ，说明 $Z_k$ 是 $X_{m-1} 与 Y_{n}$ 的 LCS
- 则 $z_k ! = y_n$ ，说明 $Z_k$ 是 $X_{m} 与 Y_{n-1 }$ 的 LCS

因此，求解LCS的问题则变成递归求解两个子问题。由于在过程中，重复的子问题多，效率低下。故此我们可以采用记忆化搜索 - 空间换时间，用数组保存中间状态，方便以后的计算。这是动态规划(DP)的核心思想

DP 求解 LCS

用二维数组 $f [i] [j]$ 记录字符串 $x 1 ， x 2 ， \cdot\cdot\cdot ， x i$ 和 $y 1 ， y 2 ， \cdot\cdot\cdot ， y j$ 的LCS 长度，则可以得到状态转移方程：
在这里插入图片描述

代码实现

const int N = 1010;

int n, m; // 串 S ，T 的长度
char S[N], T[N];
int f[N][N];

void solve()
{
    for(int i = 1; i <= n; i ++ )
    {
        for(int j = 1; j <= m; j ++ )
        	if(s[i] == T[j])
                f[i][j] = f[i - 1][j - 1] + 1;
        	else
                f[i][j] = max(f[i][j - 1], f[i - 1][j]);
        
    }
    
    printf("%d\n", f[n][m]);
    
}

以上内容尚未完全，随着今后学习的推进，我会继续对其进行补充与完善。

参考资料

[1] cs2035, Longest Common Subsequence.

[2] 最长公共子序列与最长公共子串

[3] GeeksforGeeks, Dynamic Programming | Set 29 (Longest Common Substring).

NothingAtall.

关注

1
点赞
踩
5

收藏

觉得还不错? 一键收藏
0
评论
最长公共子序列问题详解（LCS）

我们的目的是要找出字符串 X = < x1，x2，···，xm> 和 Y = < y1 ， y2，···，yn > 最长公共子序列。由于在过程中，重复的子问题多，效率低下。假设 Z = < z1， z2，···，zk > 是它们的最长公共子序列。，用数组保存中间状态，方便以后的计算。以上内容尚未完全，随着今后学习的推进，我会继续对其进行补充与完善。的子序列，是一组按从左到右顺序出现的字符，但不一定是连续的。(m) ，求出两个字符串的最长公共子序列的长度。例如，我们有一字符串。
复制链接

扫一扫