最长公共子序列问题（LCS）——Python实现

 

# 最长公共子序列问题

# 作用：求两个序列的最长公共子序列

# 输入：两个字符串数组：A和B

# 输出：最长公共子序列的长度和序列

def LCS(A,B):

    print('输入字符串数组A',A)
    print('输入字符串数组B',B);print('\n')
    n = len(A)
    m = len(B)

    # 在字符串数组A、B之前插入字符0，目的是使后面下标统一
    A.insert(0,'0')
    B.insert(0,'0')

    # 二维表L存放公共子序列的长度
    L = [ ([0]*(m+1)) for i in range(n+1) ]
    # 二维表C存放公共子序列的长度步进
    C = [ ([0]*(m+1)) for i in range(n+1) ]

    for x in range (0,n+1):
        for y in range (0,m+1):
            if (x==0 or y==0):
                L[x][y] = 0
            elif A[x] == B[y]:
                L[x][y] = ( L[x-1][y-1] + 1 )
                C[x][y] = 0
            elif L[x-1][y] >= L[x][y-1]:
                L[x][y] = L[x-1][y]
                C[x][y] = 1
            else:
                L[x][y] = L[x][y-1]
                C[x][y] = -1

        print('二维表行标：',x)
        print(L[x])

    print('\n');print('公共子序列长度二维表：');print (L)
    print('\n');print('公共子序列长度步进表：');print (C);print('\n')
    return L[n][m],C,n,m

其中返回的 L[n][m] 就是最长公共子序列的长度，以下打印序列：

def printLCS(C,A,x,y):

    if ( x == 0 or y == 0):
        return 0
    if C[x][y] == 0:
        printLCS(C,A,x-1,y-1)
        print (A[x])
    elif C[x][y] == 1:
         printLCS(C,A,x-1,y)
    else:
        printLCS(C,A,x,y-1)

输入字符串数组A、B并进行函数调用：

A = ['z', 'x', 'y', 'x', 'y', 'z']
B = ['x', 'y', 'y', 'z', 'x']
length,C,x,y = LCS(A,B)
print('最长公共子序列长度为：',length)
print('最长公共子序列为：')
printLCS(C,A,x,y)

运行结果：

二维表行标： 0
[0, 0, 0, 0, 0, 0]
二维表行标： 1
[0, 0, 0, 0, 1, 1]
二维表行标： 2
[0, 1, 1, 1, 1, 2]
二维表行标： 3
[0, 1, 2, 2, 2, 2]
二维表行标： 4
[0, 1, 2, 2, 2, 3]
二维表行标： 5
[0, 1, 2, 3, 3, 3]
二维表行标： 6
[0, 1, 2, 3, 4, 4]

公共子序列长度二维表：
[[0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 1, 1], [0, 1, 1, 1, 1, 2], [0, 1, 2, 2, 2, 2], [0, 1, 2, 2, 2, 3], [0, 1, 2, 3, 3, 3], [0, 1, 2, 3, 4, 4]]

公共子序列长度步进表：
[[0, 0, 0, 0, 0, 0], [0, 1, 1, 1, 0, -1], [0, 0, -1, -1, 1, 0], [0, 1, 0, 0, -1, 1], [0, 0, 1, 1, 1, 0], [0, 1, 0, 0, -1, 1], [0, 1, 1, 1, 0, -1]]

最长公共子序列长度为： 4
最长公共子序列为：
x
y
y
z