最长公共子序列（LCS）

Posted 2021-02-15 steven2020

【问题】 求两字符序列的最长公共字符子序列

 1 def lcs_length(x,y):
 2     m = len(x)
 3     n = len(y)
 4     c = [[0 for _ in range(n+1)] for _ in range(m+1)]
 5     for i in range(1,m+1):
 6         for j in range(1,n+1):
 7             if x[i-1] == y[j-1]:
 8                 c[i][j] = c[i-1][j-1] + 1
 9             else:
10                 c[i][j] = max(c[i-1][j],c[i][j-1])
11     return c[m][n]
12 
13 def lcs(x,y):
14     m = len(x)
15     n = len(y)
16     c = [[0 for _ in range(n+1)] for _ in range(m+1)]
17     b = [[0 for _ in range(n+1)] for _ in range(m+1)]
18     for i in range(1,m+1):
19         for j in range(1,n+1):
20             if x[i-1] == y[j-1]:
21                 c[i][j] = c[i-1][j-1] +1
22                 b[i][j] = 1
23             elif c[i-1][j] > c[i][j-1]:
24                 c[i][j] = c[i-1][j]
25                 b[i][j] = 2
26             else:
27                 c[i][j] = c[i][j-1]
28                 b[i][j] = 3
29     return c[m][n],b 
30 
31 def lcs_trackback(x,y):
32     c,b = lcs(x,y)
33     i = len(x)
34     j = len(y)
35     res = []
36     while i>0 and j>0:
37         if b[i][j] == 1:  #来自左上方=>匹配
38             res.append(x[i-1])
39             i-=1
40             j-=1
41         elif b[i][j] == 2:#来自于上方=>不匹配
42             i-=1
43         else: #==3来自于左方=>不匹配
44             j-=1
45     return "".join(reversed(res))

View Code