码迷,mamicode.com
首页 > 其他好文 > 详细

最长公共子序列lcs实现

时间:2021-04-22 16:14:26      阅读:0      评论:0      收藏:0      [点我收藏+]

标签:reverse   字符   span   ret   sed   orm   and   inf   lse   

def lcs(s1, s2):
    m = len(s1)         # 记录s1长度
    n = len(s2)         # 记录s2长度
    a = [[0 for j in range(n+1)]for i in range(m+1)]        # 得分数组
    b = [[0 for j in range(n+1)]for i in range(m+1)]        # 路径方向数组
    for i in range(1, m+1):
        for j in range(1, n+1):
            if s1[i-1] == s2[j-1]:          # 当前字符相同,左上对角加1
                a[i][j] = a[i-1][j-1] + 1
                b[i][j] = 1
            elif a[i][j-1] > a[i-1][j]:     # 之前最长序列由左边得到
                a[i][j] = a[i][j-1]
                b[i][j] = 2
            else:
                a[i][j] = a[i-1][j]         # 之前最长序列由上边得到
                b[i][j] = 3
    return a[m][n], b                       # 返回公共序列的长度及方向矩阵

def lcs_traceback(s1, s2):
    """ 路径回溯
    :param s1:
    :param s2:
    :return:
    """
    a, b = lcs(s1, s2)
    i = len(s1)
    j = len(s2)
    res = []
    while i > 0 and j > 0:      # 有一个字符为0时,循环停止
        if b[i][j] == 1:
            res.append(s1[i-1]) # 能够取到的下标是i-1
            i -= 1
            j -= 1
        elif b[i][j] == 2:
            j -= 1
        else:
            i -= 1
    print(b)
    print("{}".format(list(reversed(res))))


if __name__ == __main__:
    lcs_traceback("abcde", "bcdfe")

结果显示:

技术图片

 

最长公共子序列lcs实现

标签:reverse   字符   span   ret   sed   orm   and   inf   lse   

原文地址:https://www.cnblogs.com/demo-deng/p/14687174.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!