最长公共子序列动态规划算法实现

文章展示了如何使用Python编写算法来找到两个字符串的最长公共子序列。首先定义了一个`publ`函数尝试找出公共子序列,然后通过`lcs`函数及`printLcs`辅助函数,基于动态规划计算并打印最长公共子序列。示例中比较了字符串ababd和abgcd的最长公共子序列。
摘要由CSDN通过智能技术生成

败笔

'''        for j in range(len2):
            c = list2[j]
            print(q)
            if(b == c):
                if q == None:
                    q = c
                    print(q)
                else:
                    q = q+c
                    print('qc'+q)
            else:
                if q != None:
                    listfin.append(q)
                    q = None'''

def publ(list1,list2):
    listfin = []
    len1 = len(list1)
    len2 = len(list2)
    q = None
    for i in range(len1):
        b = list1[i]
        for j in range(len2):
            c = list2[j]
            if b == c:
                if q == None:
                    q = c
                    print(q)
                else:
                    q = q+c
                    print('qc='+q)
    print(listfin)
    len3 = len(listfin)
    for n in range(len3):
        a = None
        if a == None:
            a = listfin[n]
        else:
            if len(a) < len(listfin[n]):
                a = listfin[n]
    return 0

listb = ['a','b','c','f','r','h']
lista = ['b','c','f','r','a','b','c','h']
print(publ(lista,listb))

 是不是自己写的已经忘了

 

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]:
                print(c[i][j])
                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 c, 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 = 'ababd'
b = 'abgcd'
c, flag = lcs(a, b)


printLcs(flag, a, len(a), len(b))

结果:

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值