2013-11-21 70 views
0

我有2个序列,需要找到最长的公共子序列。不知道为什么我的功能恢复不起作用。最长的公共子序列(还原序列)

#sequenses 
A=[1,2,3] 
B=[2,3,1,5] 
#table AxB 
rows=[0]*(len(B)+1) 
table=[rows for l in range(len(A)+1)] 
for i in range(len(A)): 
    for k in range(len(B)): 
     if A[i]==B[k]: 
      table[i+1][k+1]=table[i][k]+1 
     else: 
      table[i+1][k+1]=max(table[i][k+1], table[i+1][k]) 
print(table) 
lst=[]#subsequence 
#function to restore subsequence by walking through table 
def restore(s1,s2):#s1=len(A)-1, s2=len(B)-1 
    if s1==-1 or s2==-1: 
     pass 
    elif A[s1]==B[s2]: 
     print (1) 
     lst.append(A[s1]) 
     restore(s1-1, s2-1) 
    else: 
     if table[s1][s2+1]==table[s1+1][s2+1]: 
      restore(s1-1,s2) 
     else: 
      restore(s1, s2-1) 
    return lst 

如果我做

print (restore(2,3)) 

函数返回

[] 

我认为,这个问题是索引,但无法找到它在哪里。

回答

0

这是你的第一个问题:

rows=[0]*(len(B)+1) 
table=[rows for l in range(len(A)+1)] 

Python使用引用的名单。所以当你制作一堆rows时,它实际上是对同一个rows列表的多个引用。看看这个测试代码:

>>> rows = [0]*5 
>>> table = [rows for l in range(5)] 
>>> table 
[[0, 0, 0, 0, 0], [0, 0, 0, 0, 0], [0, 0, 0, 0, 0], [0, 0, 0, 0, 0], [0, 0, 0, 0, 0]] 
>>> table[0][2] = 4 
>>> table 
[[0, 0, 4, 0, 0], [0, 0, 4, 0, 0], [0, 0, 4, 0, 0], [0, 0, 4, 0, 0], [0, 0, 4, 0, 0]] 

看看如何重复单个更改?这是因为我有5个引用相同的列表,而不是5个列表。试试这个:

table = [None]*(len(A)+1) 
for i in range(len(table)): 
    table[i] = [0]*(len(B)+1) 
+0

谢谢!我的表错了。现在所有的脚本都正常工作 –