Python 最长公共子串算法代码,Python代码库,德仔网

代码语言

代码分类

【Python】 Python 最长公共子串算法

作者:Sephiroth / 发布于2011/1/10/ 1220

<div>!/usr/bin/env python 
# find an LCS (Longest Common Subsequence). 
# *public domain*</div>
<div>def find_lcs_len(s1, s2): 
 m = [ [ 0 for x in s2 ] for y in s1 ] 
 for p1 in range(len(s1)): 
 for p2 in range(len(s2)): 
 if s1[p1] == s2[p2]: 
 if p1 == 0 or p2 == 0: 
 m[p1][p2] = 1 
 else: 
 m[p1][p2] = m[p1-1][p2-1]+1 
 elif m[p1-1][p2] < m[p1][p2-1]: 
 m[p1][p2] = m[p1][p2-1] 
 else: # m[p1][p2-1] < m[p1-1][p2] 
 m[p1][p2] = m[p1-1][p2] 
 return m[-1][-1]</div>
<div>def find_lcs(s1, s2): 
 # length table: every element is set to zero. 
 m = [ [ 0 for x in s2 ] for y in s1 ] 
 # direction table: 1st bit for p1, 2nd bit for p2. 
 d = [ [ None for x in s2 ] for y in s1 ] 
 # we don't have to care about the boundery check. 
 # a negative index always gives an intact zero. 
 for p1 in range(len(s1)): 
 for p2 in range(len(s2)): 
 if s1[p1] == s2[p2]: 
 if p1 == 0 or p2 == 0: 
 m[p1][p2] = 1 
 else: 
 m[p1][p2] = m[p1-1][p2-1]+1 
 d[p1][p2] = 3 # 11: decr. p1 and p2 
 elif m[p1-1][p2] < m[p1][p2-1]: 
 m[p1][p2] = m[p1][p2-1] 
 d[p1][p2] = 2 # 10: decr. p2 only 
 else: # m[p1][p2-1] < m[p1-1][p2] 
 m[p1][p2] = m[p1-1][p2] 
 d[p1][p2] = 1 # 01: decr. p1 only 
 (p1, p2) = (len(s1)-1, len(s2)-1) 
 # now we traverse the table in reverse order. 
 s = [] 
 while 1: 
 print p1,p2 
 c = d[p1][p2] 
 if c == 3: s.append(s1[p1]) 
 if not ((p1 or p2) and m[p1][p2]): break 
 if c & 2: p2 -= 1 
 if c & 1: p1 -= 1 
 s.reverse() 
 return ''.join(s)</div>
if __name__ == '__main__': 
 print find_lcs('abcoisjf','axbaoeijf') 
 print find_lcs_len('abcoisjf','axbaoeijf')</div>

评论列表

本站所提供的代码，版权归原作者所有，若有侵犯作者版权，请与我们联系，我们将立即删除或修改。谢谢!
本站所有代码发布及提供者。

试试其它关键字

　最长公共子串　

同语言下

可能有用的

Sephiroth贡献的其它代码(13)