C++版 - Lintcode 77-Longest Common Subsequence最长公共子序列(LCS) - 题解

本文介绍使用C++实现最长公共子序列(Longest Common Subsequence, LCS)算法的方法。通过动态规划解决两个字符串间的LCS问题,分析了不同情况下的递推公式,并提供已通过测试的代码。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

版权声明:本文为博主Bravo Yeung(知乎UserName同名)的原创文章,欲转载请先私信获博主允许,转载时请附上网址
http://blog.csdn.net/lzuacm

C++版 - Lintcode 77-Longest Common Subsequence最长公共子序列(LCS) - 题解

在线提交(不支持C#):
https://www.lintcode.com/problem/longest-common-subsequence/

题目描述

一个字符串的一个子序列是指,通过删除一些(也可以不删除)字符且不干扰剩余字符相对位置所组成的新字符串。例如,”ACE” 是 “ABCDE” 的一个子序列,而 “AEC” 不是)。
给出两个字符串,找到最长公共子序列(LCS),返回LCS的长度。

说明

样例
给出“ABCD”“EDCA”,这个LCS是 “A” (或 D或C),返回1
给出 “ABCD”“EACB”,这个LCS是“AC”返回 2

注意
序列可以不连续。


  ● Difficulty: Medium
  • Total Accepted: 18202

  • Total Submitted: 45985

  • Accepted Rate: 39%

Tags:
Longest Common Subsequence
LintCode Copyright
Dynamic Programming(DP)

分析:
将算式的计算结果记录在内存中,需要时直接调用该结果,从而避免无用的重复计算,提高处理效率,这在程序和算法设计中是一种行之有效的手法。动态规划就是这类手法之一。

事实上动态规划是一种记忆化递归(memorized recursive),缓存部分重要数据。另外,动态规划法可以建立递归式,通过循环顺次求出最优解。

为方便说明,这里我们用XiXi代表{x1,x2,,xix1,x2,⋯,xi},用YjYj代表{y1,y2,,yjy1,y2,⋯,yj }。那么,求长度分别为m、n的两个序列XY的LCS,就相当于求XmXmYnYn的LCS。我们将其分割为局部问题进行分析。

首先,求XmXmYnYn的LCS时要考虑以下两种情况:

  • xm=ynxm=yn时,在Xm1Xm−1Yn1Yn−1的LCS后面加上xm(=yn)xm(=yn)就是xmxmynyn的LCS。

    举个例子,X=(a,b,c,c,d,a),Y={a, b, c, b, a}时xm=ynxm=yn,所以在Xm1Xm−1Yn1Yn−1的LCS({a, b,c})后面加上xmxm {=a) 即为XmXmYnYn的LCS。

  • xmynxm≠yn时,Xm1Xm−1YnYn的LCS和XmXmYn1Yn−1的LCS中更长的一方就是XmXmYnYn的LCS。

    举个例子,X={a,b,c,c,d}, Y={a,b,c,b,a}时,Xm1Xm−1YnYn的LCS为{a,b,c),XmXmYnYn的LCS为{a,b,c,b},因此XmXmYn1Yn−1的LCS就是XmXmYnYn的LCS。

这个算法对XiXiYjYj同样适用。于是可准备下述函数,用来求解LCS的局部问题。

c[m+1][n+1]: 该二维数组中,c[i][j] 代表XiXiYjYj的LCS的长度

c[i][j] 的值可由下述递推公式(Recursive Formula)求得。

c[i][j]=0c[i1][j1]+1max(c[i][j1],c[i1][j])i=0 || j=0i,j>0 and xi=yji,j>0 and xiyjc[i][j]={0i=0 || j=0c[i−1][j−1]+1i,j>0 and xi=yjmax(c[i][j−1],c[i−1][j])i,j>0 and xi≠yj

基于上述变量和公式,可以用动态规划法求序列XY的LCS。

已AC代码如下:

class Solution {
public:
    int longestCommonSubsequence(string &A, string &B) {
        int m = A.size();
        int n = B.size();

        int **c = (int **)malloc((m+1) * sizeof(int *));
        for (int i = 0; i < m + 1; i++)
            c[i] = (int *)malloc((n+1) * sizeof(int));

        int max1 = 0;
        A = ' ' + A;
        B = ' ' + B;
        for (size_t i = 1; i <= m; i++)
            c[i][0] = 0;
        for (size_t j = 1; j <= n; j++)
            c[0][j] = 0;

        for (size_t i = 1; i <= m; i++)
        {
            for (size_t j = 0; j <= n; j++)
            {
                if (A[i] == B[j])
                {
                    c[i][j] = c[i - 1][j - 1] + 1;
                }
                else
                    c[i][j] = max(c[i][j - 1], c[i - 1][j]);
                max1 = max(max1, c[i][j]);
            }
        }
        return max1;
    }
};

Rank:
您的提交打败了 92.60% 的提交.

扩展阅读:
最长公共子序列问题 - Fogsail Chen - SegmentFault 思否
https://segmentfault.com/a/1190000008521545

转载于:https://www.cnblogs.com/enjoy233/p/10408692.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值