python实现最长公共子序列的求解

最新推荐文章于 2024-09-19 08:38:15 发布

勉旃

最新推荐文章于 2024-09-19 08:38:15 发布

阅读量2.3k

点赞数 2

分类专栏：数据结构（python）

本文链接：https://blog.csdn.net/yushupan/article/details/82350748

版权

数据结构（python）专栏收录该内容

9 篇文章 4 订阅

订阅专栏

（待完善...）

最长公共子序列是动态规划基本题目，下面按照动态规划基本步骤解出来。

1.找出最优解的性质，并刻划其结构特征

序列a共有m个元素，序列b共有n个元素，如果a[m-1]==b[n-1]，那么a[:m]和b[:n]的最长公共子序列长度就是a[:m-1]和b[:n-1]的最长公共子序列长度+1；如果a[m-1]!=b[n-1]，那么a[:m]和b[:n]的最长公共子序列长度就是MAX（a[:m-1]和b[:n]的最长公共子序列长度，a[:m]和b[:n-1]的最长公共子序列长度）。

2.递归定义最优值

$c[i,j] = \left\{\begin{matrix} 0 & i=0 \, \, or\, \, j =0\\ c[i-1, j-1]+1& i,j>0\! \; \; and\, \: x_{i} = y_{j}\\ max(c[i,j-1],c[i-1,j])&i,j>0\, \: \: and\: \: x_{i} \neq y_{j} \end{matrix}\right.$

3.自底向上根据最优值的信息来构造最优解

def lcs(a, b):
    lena = len(a)
    lenb = len(b)
    c = [[0 for i in range(lenb + 1)] for j in range(lena + 1)]
    flag = [[0 for i in range(lenb + 1)] for j in range(lena + 1)]
    for i in range(lena):
        for j in range(lenb):
            if a[i] == b[j]:
                c[i + 1][j + 1] = c[i][j] + 1
                flag[i + 1][j + 1] = 'ok'
            elif c[i + 1][j] > c[i][j + 1]:
                c[i + 1][j + 1] = c[i + 1][j]
                flag[i + 1][j + 1] = 'left'
            else:
                c[i + 1][j + 1] = c[i][j + 1]
                flag[i + 1][j + 1] = 'up'
    return flag

def printLcs(flag, a, i, j):
    if i == 0 or j == 0:
        return
    if flag[i][j] == 'ok':
        printLcs(flag, a, i - 1, j - 1)
        print(a[i - 1], end='')
    elif flag[i][j] == 'left':
        printLcs(flag, a, i, j - 1)
    else:
        printLcs(flag, a, i - 1, j)

a = 'ABCBDAB'
b = 'BDCABA'
flag = lcs(a, b)
printLcs(flag, a, len(a), len(b))