字符串中最长公共子序列和最长公共子串

1、最长公共子序列

例：BDCABA和ABCBDAB，最长公共子序列为 BCBA。

DP解法：字符串1的长度为length1，字符串2的长度为length2，建立一个二维数组C[length1][length2]用来记录最长公共子序列的长度，状态转移方程为

整个状态转移以及二维数组矩阵C为

def find_lcseque(s1, s2):
    # 生成字符串长度加1的0矩阵，m用来保存对应位置匹配的结果
    length1 = len(s1)
    length2 = len(s2)
    result = ''
    C = [[0 for x in range(length2 + 1)] for y in range(length1 + 1)]
    for i in range(length1):
        for j in range(length2):
            if s1[i] == s2[j]:
                C[i+1][j+1] = C[i][j] + 1
            else:
                C[i+1][j+1] = max(C[i][j+1], C[i+1][j])
    print("max length: ", C[length1][length2])

find_lcseque('ABCBDAB', 'BDCABA')

2、最长公共子串

子串是一种特殊的子序列，同样可以由DP做。

二维数组的定义如上，状态转移方程为

最大长度为result，初始值为0，当时，result = max(C[i][j], result)

def find_lcsSubstr(s1, s2):
    # 生成字符串长度加1的0矩阵，m用来保存对应位置匹配的结果
    length1 = len(s1)
    length2 = len(s2)
    substr_length = 0
    C = [[0 for x in range(length2 + 1)] for y in range(length1 + 1)]
    for i in range(length1):
        for j in range(length2):
            if s1[i]==s2[j]:
                C[i][j] = C[i-1][j-1] + 1
                substr_length = max(C[i][j], substr_length)
            else:
                C[i][j] = 0

    print("Substr max length: ", substr_length)

find_lcsSubstr('ABCBDAB', 'BDCABA')