Python求两个字符串最长公共子序列代码实例


Posted in Python onMarch 05, 2020

一、问题描述

给定两个字符串,求解这两个字符串的最长公共子序列(Longest Common Sequence)。比如字符串1:BDCABA;字符串2:ABCBDAB。则这两个字符串的最长公共子序列长度为4,最长公共子序列是:BCBA

二、算法求解

这是一个动态规划的题目。对于可用动态规划求解的问题,一般有两个特征:①最优子结构;②重叠子问题

①最优子结构

设X=(x1,x2,...,xn)和Y=(y1,y2,...,ym)是两个序列,将X和Y的最长公共子序列记为LCS(X,Y)

找出LCS(X,Y)就是一个最优化问题。因为,我们需要找到X和Y中最长的那个公共子序列。而要找X和Y的LCS,首先考虑X的最后一个元素和Y的最后一个元素。

⑴如果xn=ym,即X的最后一个元素与Y的最后一个元素相同,这说明该元素一定位于公共子序列中。因此,现在只需要找:LCS(Xn-1,Ym-1)

LCS(Xn-1,Ym-1)就是原问题的一个子问题。为什么叫子问题?因为它的规模比原问题小。

为什么是最优的子问题?因为我们要找的是Xn-1和Ym-1的最长公共子序列啊。最长的!换句话说就是最优的那个。

⑵如果xn!=ym,这下要麻烦一点,因为它产生了两个子问题:LCS(Xn-1,Ym)和LCS(Xn,Ym-1)

因为序列X和序列Y的最后一个元素不相等,那说明最后一个元素不可能是最长公共子序列中的元素。

LCS(Xn-1,Ym)表示:最长公共序列可以在(x1,x2,...xn-1)和(y1,y2,...,ym)中找。

LCS(Xn,Ym-1)表示:最长公共序列可以在(x1,x2,...xn)和(y1,y2,...,ym-1)中找。

求解上面两个子问题,得到的公共子序列谁最长,那谁就是LCS(X,Y)。用数学表示就是:

LCS=max{LCS(Xn-1,Ym),LCS(Xn,Ym-1)}

由于条件⑴和⑵考虑到了所有可能的情况。因此,我们成功的把原问题转化成了三个规模更小的问题。

②重叠子问题

重叠子问题是什么?就是说原问题转化成子问题后,子问题中有相同的问题。

原问题是:LCS(X,Y)。子问题有❶LCS(Xn-1,Ym-1)❷ LCS(Xn-1,Ym)❸ LCS(Xn,Ym-1)

乍一看,这三个问题是不重叠的。可本质上它们是重叠的,因为它们只重叠了一大部分。举例:

第二个子问题:LCS(Xn-1,Ym)就包含了问题❶LCS(Xn-1,Ym-1),为什么?

因为,当Xn-1和Ym的最后一个元素不相同时,我们又需要将LCS(Xn-1,Ym-1)进行分解:分解成:LCS(Xn-1,Ym-1)和LCS(Xn-2,Ym)

也就是说:在子问题的继续分解中,有些问题是重叠的。

由于像LCS这样的问题,它具有重叠子问题的性质,因此:用递归来求解就太不划算了。国为采用递归,它重复地求解了子问题,而且需要注意的是,所有子问题加起来的个数是指数级的。

那么问题来了,如果用递归求解,有指数级个子问题,故时间复杂度是指数级的。这指数级个子问题,难道用了动态规划,就变成多项式时间了??

关键是采用动态规划时,并不需要去一一计算那些重叠了的子问题。或者说:用了动态规划之后,有些子问题是通过“查表”直接得到的,而不是重新又计算一遍得到的。举个例子:比如求Fib数列。

Python求两个字符串最长公共子序列代码实例

求fib(5),分解成了两个子问题:fib(4)和fib(3),求解fib(4)和fib(3)时,又分解了一系列的小问题...

从图中可以看出:根的左右子树:fib(4)和fib(3)下,是有很多重叠的!比如,对于fib(2),它就一共出现了三次。如果用递归来求解,fib(2)就会被计算三次,而用DP(Dynamic Programming)动态规划,则fib(2)只会计算一次,其他两次则是通过“查表”直接求得。而且,更关键的是:查找求得该问题的解之后,就不需要再继续去分解该问题了。而对于递归,是不断地将问题解,直到分解为基准问题(fib(0)或者fib(1))

说了这么多,还是写下最长公共子序列的递归式才完整。

Python求两个字符串最长公共子序列代码实例

C[i,j]表示:(x1,x2,...,xi)和(y1,y2,...,yj)的最长公共子序列的长度。公式的具体解释可参考《算法导论》动态规划章节

三、LCS Python代码实现

#! /usr/bin/env python3
# -*- coding:utf-8 -*-

# Author  : mayi
# Blog   : http://www.cnblogs.com/mayi0312/
# Date   : 2019/5/16
# Name   : test03
# Software : PyCharm
# Note   : 用于实现求解两个字符串的最长公共子序列


def longestCommonSequence(str_one, str_two, case_sensitive=True):
  """
  str_one 和 str_two 的最长公共子序列
  :param str_one: 字符串1
  :param str_two: 字符串2(正确结果)
  :param case_sensitive: 比较时是否区分大小写,默认区分大小写
  :return: 最长公共子序列的长度
  """
  len_str1 = len(str_one)
  len_str2 = len(str_two)
  # 定义一个列表来保存最长公共子序列的长度,并初始化
  record = [[0 for i in range(len_str2 + 1)] for j in range(len_str1 + 1)]
  for i in range(len_str1):
    for j in range(len_str2):
      if str_one[i] == str_two[j]:
        record[i + 1][j + 1] = record[i][j] + 1
      elif record[i + 1][j] > record[i][j + 1]:
        record[i + 1][j + 1] = record[i + 1][j]
      else:
        record[i + 1][j + 1] = record[i][j + 1]

  return record[-1][-1]

if __name__ == '__main__':
  # 字符串1
  s1 = "BDCABA"
  # 字符串2
  s2 = "ABCBDAB"
  # 计算最长公共子序列的长度
  res = longestCommonSequence(s1, s2)
  # 打印结果
  print(res) # 4

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Python中针对函数处理的特殊方法
Mar 06 Python
使用pdb模块调试Python程序实例
Jun 02 Python
python使用电子邮件模块smtplib的方法
Aug 28 Python
tensorflow实现对图片的读取的示例代码
Feb 12 Python
Python字符串逆序的实现方法【一题多解】
Feb 18 Python
Python中使用遍历在列表中添加字典遇到的坑
Feb 27 Python
python实现飞机大战游戏
Oct 26 Python
ubuntu上安装python的实例方法
Sep 30 Python
Python绘制动态水球图过程详解
Jun 03 Python
Python如何爬取b站热门视频并导入Excel
Aug 10 Python
Python confluent kafka客户端配置kerberos认证流程详解
Oct 12 Python
python中delattr删除对象方法的代码分析
Dec 15 Python
Python操作MongoDb数据库流程详解
Mar 05 #Python
Python文字截图识别OCR工具实例解析
Mar 05 #Python
win10下opencv-python特定版本手动安装与pip自动安装教程
Mar 05 #Python
python+OpenCV实现图像拼接
Mar 05 #Python
windows下Pycharm安装opencv的多种方法
Mar 05 #Python
解决pycharm中opencv-python导入cv2后无法自动补全的问题(不用作任何文件上的修改)
Mar 05 #Python
浅谈Python的方法解析顺序(MRO)
Mar 05 #Python
You might like
PHP学习笔记之数组篇
2011/06/28 PHP
PHP SplObjectStorage使用实例
2015/05/12 PHP
php getcwd与dirname(__FILE__)区别详解
2016/09/24 PHP
visual studio code 调试php方法(图文详解)
2017/09/15 PHP
PHP实现无限极分类的两种方式示例【递归和引用方式】
2019/03/25 PHP
浅谈laravel中的关联查询with的问题
2019/10/10 PHP
laravel实现图片上传预览,及编辑时可更换图片,并实时变化的例子
2019/11/14 PHP
PHP笛卡尔积实现原理及代码实例
2020/12/09 PHP
js 数值转换为3位逗号分隔的示例代码
2014/02/19 Javascript
jQuery中:nth-child选择器用法实例
2014/12/31 Javascript
javascript中offset、client、scroll的属性总结
2015/08/13 Javascript
AngularJS 让人爱不释手的八种功能
2016/03/23 Javascript
javaScript中的原型解析【推荐】
2016/05/05 Javascript
ES6中Array.includes()函数的用法
2017/09/20 Javascript
基于 Vue 实现一个酷炫的 menu插件
2017/11/14 Javascript
javascript中的replace函数(带注释demo)
2018/01/07 Javascript
vue+webpack实现异步组件加载的方法
2018/02/03 Javascript
原生JS实现获取及修改CSS样式的方法
2018/09/04 Javascript
layui默认选中table的CheckBox复选框方法
2019/09/19 Javascript
vue单文件组件无法获取$refs的问题
2020/06/24 Javascript
[03:17]DOTA2-DPC中国联赛1月29日Recap集锦
2021/03/11 DOTA
[41:52]DOTA2-DPC中国联赛 正赛 CDEC vs Dynasty BO3 第二场 2月22日
2021/03/11 DOTA
让python在hadoop上跑起来
2016/01/27 Python
python机器学习案例教程——K最近邻算法的实现
2017/12/28 Python
Python 中PyQt5 点击主窗口弹出另一个窗口的实现方法
2019/07/04 Python
python把ipynb文件转换成pdf文件过程详解
2019/07/09 Python
解决Pytorch自定义层出现多Variable共享内存错误问题
2020/06/28 Python
美国最顶级的精品店之一:Hampden Clothing
2016/12/22 全球购物
Missguided美国官网:英国时尚品牌
2018/01/18 全球购物
诗普兰迪官方网站:Splendid
2018/09/18 全球购物
上课看小说检讨书
2014/02/22 职场文书
年度考核自我鉴定
2014/03/19 职场文书
大学生应聘求职信
2014/05/26 职场文书
组织生活会发言材料
2014/12/15 职场文书
安全先进个人材料
2014/12/29 职场文书
Python与C++中梯度方向直方图的实现
2022/03/17 Python