最长公共子序列

题目:求两个str的公共子序列长度,在基因对中用的多

输入:str1 = “ABCDEF” str2 = “AEBDEDF”

输出:5

思路:动态规划

状态转移方程:

1
2
3
4
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])

状态:dp[i,j]表示str[:i]和str[:j]的公共子序列长度

边缘化处理:dp的0行和0列为0

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
class :
def find_lcs(self, str1, str2):
size1 = len(str1)
size2 = len(str2)

dp = [[0 for i in range(size2+1)] for j in range(size1+1)]
# 动态规划过程
for i in range(1,size1+1):
for j in range(1,size2+1):
if str1[i-1] == str2[j-1]:
dp[i][j] = dp[i - 1][j - 1] + 1
else:
dp[i][j] = max(dp[i - 1][j], dp[i][j - 1])
return dp[-1][-1]

# 例子
str1 = "ABCDEF"
str2 = "AEBDEDF"
print(Solution().find_lcs(str1, str2))