求最长公共子串长度与最长公共子串问题/求最长公共子序列长度与最长公共子序列问题

求最长公共子串的长度与最长公共子串：

import numpy as np

str1 = 'abaccb'
str2 = 'aaccab'

tmp_str1 = str1
tmp_str2 = str2

str1 = '0' + str1
str2 = '0' + str2

len1 = len(str1)
len2 = len(str2)

dp = np.zeros((len1, len2), dtype=np.int)

for i in range(1, len1):
    for j in range(1, len2):
        if str1[i] == str2[j]:
            dp[i][j] = dp[i - 1][j - 1] + 1
        else:
            dp[i][j] = 0

print(dp)

# 最长公共子串的长度
max_length = dp.max()

# 获取矩阵中最大值的行下标
index = np.where(dp == np.max(dp))[0][0]

print('最长公共子串的长度：', max_length)
print('最长公共子串：', tmp_str1[index - max_length:index])

求最长公共子序列的长度与最长公共子序列：

import numpy as np

str1 = 'abcdf'
str2 = 'abcdaf'

tmp_str1 = str1
tmp_str2 = str2

str1 = '0' + str1
str2 = '0' + str2

len1 = len(str1)
len2 = len(str2)

dp = np.zeros((len1, len2), dtype=np.int)

for i in range(1, len1):
    for j in range(1, len2):
        if str1[i] == str2[j]:
            dp[i][j] = dp[i - 1][j - 1] + 1
        else:
            dp[i][j] = max(dp[i - 1][j], dp[i][j - 1])

print(dp)

print('最长公共子序列长度', dp[-1][-1])

# (6, 7)
# print(dp.shape)

li = []

i = dp.shape[0] - 1
j = dp.shape[1] - 1

while i > 0 and j > 0:
    if dp[i][j] > dp[i - 1][j] and dp[i][j] > dp[i][j - 1]:
        # 来自坐上角的，则需要将该字符加入，表示该字符是公共子字符
        li.append(str1[i])
        i -= 1
        j -= 1
    elif dp[i][j] == dp[i][j - 1]:
        # 来自左边的
        j -= 1
    elif dp[i][j] == dp[i - 1][j]:
        # 来自上边的
        i -= 1

print('最长公共子序列：', li[::-1])

总结思路：

求公共子序列思路：
从后面往前推，
1、如果当前值=当前值左上角那个值+1，则该字符就是公共子序列中的一个字符
2、如果当前值=当前值的左边那个值，则进行左移
3、如果当前值=当前值的上边那个纸，则进行上移


求最长公共子串的思路：
1、获取矩阵中的最大值（也即是最长公共子串的长度值）
2、通过numpy方法中的np.where(np==np.max(dp))[0][0]，获取矩阵中最大值的行下标
3、最后进行截取字符串，截取字符串的截止下标是2中获得行下标，截取长度为最长公共子串的长度