77. 最长公共子序列 动态规划 python

77_longest_common_subsequence


给出两个字符串,找到最长公共子序列(LCS),返回LCS的长度。

说明

最长公共子序列的定义:

  • 最长公共子序列问题是在一组序列(通常2个)中找到最长公共子序列(注意:不同于子串,LCS不需要是连续的子串)。该问题是典型的计算机科学问题,是文件差异比较程序的基础,在生物信息学中也有所应用。
  • https://en.wikipedia.org/wiki/Longest_common_subsequence_problem
样例

给出"ABCD" 和 "EDCA",这个LCS是 "A" (或 D或C),返回1

给出 "ABCD" 和 "EACB",这个LCS是"AC"返回 2


class Solution:
    """
    @param A: A string
    @param B: A string
    @return: The length of longest common subsequence of A and B
    """
    def longestCommonSubsequence(self, A, B):
        
        if len(A)==0 or len(B)==0:
            return 0 
        arr =[[0 for col in range(len(B)+1)] for row in range(len(A)+1)]
        for i in range(1,len(A)+1):
            for j in range(1,len(B)+1):
                if A[i-1] == B[j-1]:
                    arr[i][j] = arr[i-1][j-1] +1
                else:
                    arr[i][j] = max(arr[i-1][j],arr[i][j-1])
        return arr[i][j],arr
        
    def getlongestCommonSubsequence(self, A, B,arr):
        rlt = []
        i = len(arr) - 1
        j = len(arr[0]) - 1
        while(i!= 0 and j!=0):
            if(arr[i][j]!=arr[i-1][j]):
                rlt.append(A[i-1])
                i-=1
                j-=1
                while(arr[i][j]==arr[i][j-1]):
                    j-=1
            else:
                i-=1
        
        return rlt
        
    



import numpy as np
a = Solution()
v,arr= a.longestCommonSubsequence('bedaacbade','dccaeedbeb')
print v
print np.array(arr)
print a.getlongestCommonSubsequence('bedaacbade','dccaeedbeb',arr)
5
[[0 0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 1 1 1]
 [0 0 0 0 0 1 1 1 1 2 2]
 [0 1 1 1 1 1 1 2 2 2 2]
 [0 1 1 1 2 2 2 2 2 2 2]
 [0 1 1 1 2 2 2 2 2 2 2]
 [0 1 2 2 2 2 2 2 2 2 2]
 [0 1 2 2 2 2 2 2 3 3 3]
 [0 1 2 2 3 3 3 3 3 3 3]
 [0 1 2 2 3 3 3 4 4 4 4]
 [0 1 2 2 3 4 4 4 4 5 5]]
['e', 'd', 'a', 'c', 'd']

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值