最长公共子序列(动态规划)

目录

1.问题描述

2.问题分析

3.手工运算

4.核心源码

---

1.问题描述

给出两个字符串，求出这样的一个最长的公共子序列的长度:子序列中的每个字符都能在两个原串中找到而且每个字符的先后烦序和原串中的先后顺序一致。

Sample Input

abcfbc  abfcab

programming  contest

abcd  mnp

Sample Output

4

2

0

---

2.问题分析

输入两个串s1,s2,

  设MaxLen(i,j)表示: 

        s1的左边i个字符形成的子串，与s2左边的j个字符形成的子串的最长公共子序列的长度            (i,j从0开始算)

        MaxLen(i,j)就是本题的“状态”

假定len1 = strlen(s1),len2 = strlen(s2)那么题目就是要求MaxLen(len1,len2)

显然:

MaxLen(n,0) = o ( n= 0...len1)

MaxLen(0,n) = o ( n=0...len2)

递推公式:

if ( s1[i-1] == s2[j-1] )//ls1的最左边字符是s1[0]

        MaxLen(i.j)= MaxLen(i-1,j-1)+ 1;

else

        MaxLen(i,j)=Max(MaxLen(i,j-1),MaxLen(i-1,j));

时间复杂度O(mn)，m,n是两个字串长度

 

S1[i-1]!= s2[j-1]时，MaxLen(S1,S2)不会比MaxLen(S1,S2,-1)和MaxLen(S1;-1,S2)两者之中任何一个小，也不会比两者都大。

---

3.手工运算

给大家解释一下这个规则，下面的例题中将会反复运用到这个规则

（1）LU：当对角线中右下角的数字比左上角的数字大1时，在rec表格中填入LU

                                                 

 

（2）U：当上面的数字与下面的数字一样大的时候，在rec表格中填入U 

                                                           

（3）L：当左边的数字与右边的数字一样大的时候，在rec表格中填入L

                                                   

 

基本思想：C表怎么填

（1）当扫描的两个字符串的字母不相同则在他的左边和上面的两个表格中找较大值，将其填入这个表格

 

 比如说这问号的表格中应该填入3和4的较大值4

（2）当扫描的两个字符串的字母相同则在他的左上方空格的基础上加1

比如说这个问号的表格中应该填入5 

接下来我们用一道例题给大家讲解一下

 

我们先对长度为0的C表格进行初始化

  

接下来我们把两个指针都置于两个字符串的开头进行扫描。两个字符串的字母不相同

指针向后移动一位

 

两个字符串中字母相同，在左上角的基础上加1，在rec表格中填入LU

 

两个字符串中字母不相同，填入左边和上面的较大值，最后填入的是左边的1，在rec表格中填入L

 

按照这个指导思想可以得出最后的结果

 

---

4.核心源码

class Solution {
public:
    int longestCommonSubsequence(string text1, string text2) {
        if(text1.size()<text2.size())
        {
            string tmp = text1;
            text1 = text2;
            text2 = tmp;
        }
        int m = text1.size();
        int n = text2.size();
        vector<int> dp(m + 1, 0);
        vector<int> pre(m + 1, 0);
        for (int i = 1; i <= m; i++) {
            char c1 = text1.at(i - 1);
            for (int j = 1; j <= n; j++) {
                char c2 = text2.at(j - 1);
                if (c1 == c2) {
                    dp[j] = pre[j - 1] + 1;
                } else {
                    dp[j] = max(dp[j-1], pre[j]);
                }
            }
            pre = dp;
        }
        return dp[n];
    }
};