Python最长公共子序列（动态规划）

最长公共子序列（北大POJ1458)

题目描述：
给出两个字符串，求出这样的一个最长的公共子序列的长度：子序列中的每个字符都能在两个原串中找到， 而且每个字符的先后顺序和原串中的先后顺序一致。

Sample Input

abcfbc abfcab

programming    contest

abcd    mnp

Sample Output

4

2

0

解题思路：设输入的两个串为s1,s2, 设MaxLen(i,j)表示:：s1的左边i个字符形成的子串，与s2左边的j个字符形成的子串的最长公共子序列的长度(i,j从0 开始算），则MaxLen(i,j) 就是本题的“状态” （如果还不懂“状态”是什么意思，可以参考我上一篇文章）

假定 len1 = strlen(s1),len2 = strlen(s2），那么题目就是要求   MaxLen(len1,len2)

显然：

MaxLen(n,0)  = 0  ( n= 0…len1）

MaxLen(0,n)  = 0  ( n=0…len2）

于是，我们可以得到如下的递推公式：  

if ( s1[i-1] == s2[j-1] ) //s1的最左边字符是s1[0]   
	MaxLen(i,j) = MaxLen(i-1,j-1) + 1; 
else  
	MaxLen(i,j) = Max(MaxLen(i,j-1),MaxLen(i-1,j) );

*时间复杂度是O(len(s1)len(s2))

S1[i-1]!= s2[j-1]时，MaxLen(S1,S2)不会比MaxLen(S1,S2j-1) 和MaxLen(S1i-1,S2)两者之中任何一个小，也不会比两者都大。

参考代码：

#encoding:utf-8
"""
example:
Sample Input
    abcfbc abfcab
    programming    contest
    abcd    mnp
Sample Output
    4
    2
    0
"""
while True:
    s1,s2 = input().strip().split()
    #find the longest using DP
    l1 = len(s1)
    l2 = len(s2)
    #创建二维数组Maxlen，Maxlen[i][j]代表s1前i个子字符串与s2前j个子字符串的最长公共子序列
    Maxlen = [[0 for j in range(l2+1)]for i in range(l1+1)]
    for i in range(1,l1+1):
        for j in range(1,l2+1):
            if s1[i-1] == s2[j-1]:
                Maxlen[i][j] = Maxlen[i-1][j-1]+1
            else:
                Maxlen[i][j] = max(Maxlen[i-1][j],Maxlen[i][j-1])
    print(Maxlen[l1][l2])

注意：
我们在这里创建二维数组是用的

aList = [[0 for j in range(m)]for j in range(n)]

而不是

aList = [[0]*m]*n

这是因为后者是浅拷贝，在使用时会出现以下情况

>>> lists = [[]] * 3
>>> lists
[[], [], []]
>>> lists[0].append(3)
>>> lists
[[3], [3], [3]]

*[[]]是一个含有一个空列表元素的列表,所以[[]]3表示3个指向这个空列表元素的引用,修改任何一个元素都会改变整个列表

所以需要用前种方式进行创建多维数组,以免浅拷贝。

reference：
1、https://blog.csdn.net/baidu_28312631/article/details/47426445?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-2.channel_param&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-2.channel_param
2、https://blog.csdn.net/yockie/article/details/46127829