python最长回文串算法


Posted in Python onJune 04, 2018

给定一个字符串,要求在这个字符串中找到符合回文性质的最长子串。所谓回文性是指诸如 “aba”,"ababa","abba"这类的字符串,当然单个字符以及两个相邻相同字符也满足回文性质。

看到这个问题,最先想到的解决方法自然是暴力枚举,通过枚举字符串所有字串的起点,逐一判断满足回文性的子串,记录长度并更新最长长度。显然这种算法的时间复杂度是很高的,最坏情况可以达到O(N*N)。所以呢,这里提出一个优化的方案,通过枚举字符串子串的中心而不是起点,向两边同时扩散,依然是逐一判断子串的回文性。这种优化算法比之前的算法在最坏的情况下(即只有一种字符的字符串)效率会有很大程度的上升。

由上述的优化方案,我们知道了枚举中心要比枚举起点效率要好,然而这并不是最优的算法。由于枚举中心的算法同时影响的是中心两边的字符,所以我们可以通过枚举中心的左边字符作为中心的子串的回文性判断枚举中心右边的字符作为中心得子串的回文性,这就是manacher算法。

manacher算法思想非常巧妙,首先遍历字符串,假设 i 为枚举中心,则 j (j<i) 为中心的最长回文子串长度发f[j] 便已经求出,此时 j 的影响范围便是[j-f[j]/2,j+f [j]] 。为了使左边的字符 j 对枚举中心右边的影响最大,需要使 j+f[j]/2 最大。找到满足j+f[j]/2最大的 j 之后,若 i 在[j,j+f[j]/2]中,则分两种情况:

1 . i 关于 j 对称的字符i'的影响范围完全包含在j的影响范围内,则由于回文性,i 的影响范围大于等于i'的影响范围,即f[i]>=f[i']

2. i 关于 j 对称的字符i'的影响范围不完全包含在j的影响范围内,此时i的右侧影响范围大于等于[j-f[j]/2,i'],即i+f[i]/2>=i'-j+f[j]/2

由于对称性,可得i+i" = 2*j。因此第一种情况下,f[i]>=f[2*j-i];第二种情况下,f[i]>=f[j]+2*j-2*i。

综上1,2,可得f[i]>=min(f[2*j-i],f[j]+2*j-2*i)。由于i右边存在未遍历的字符,因此在此基础上,继续向两边扩展,直到找到最长的回文子串。

若i依然在j+f[j]/2后面,则表示i没有被前面的字符的影响,只能逐一的向两边扩展。

这个算法由于只需遍历一遍字符串,扩展的次数也是有限的,所以时间复杂度可以达到O(N)。

下面是Pthon3的程序,为了检测算法的效率,依然提供最初的暴力枚举算法作为最坏算法的参照。

python代码:

#求最长回文串类 
class LPS:      
 #初始化,需要提供一个字符串 
 def __init__(self,string): 
  self.string = string 
  self.lens = len(self.string) 
  
 #暴力枚举:作为算法效率参照 
 def brute_force(self): 
  maxcount = 0 
  for j in range(self.lens):      
   for k in range(j,self.lens): 
    count = 0 
    l,m = j,k 
    while m>=l: 
     if self.string[l]==self.string[m]: 
      l,m = l+1,m-1 
     else: 
      break 
    if m<l: 
     count = k-j+1 
    if count>maxcount : 
     maxcount = count 
  return maxcount 
  
 #优化版:枚举子串中心 
 def brute_force_opti(self): 
  maxcount = 0 
  if self.lens == 1:        #只有一个字符直接返回1 
   return 1 
  for j in range(self.lens-1):     #枚举中心 
   count,u = 1,j 
   #对于奇数子串,直接扩展 
   for k in range(1,j+1):      #两边扩展 
    l,m = u+k,j-k 
    if (m>=0)&(l<self.lens): 
     if(self.string[l]==self.string[m]): 
      count += 2 
     else: 
      break 
   if count>maxcount :       #更新回文子串最长长度 
    maxcount = count 
   if self.string[j]==self.string[j+1]:  #处理偶数子串,将两个相邻相同元素作为整体 
    u,count= j+1,2 
   for k in range(1,j+1):      #两边扩展 
    l,m = u+k,j-k 
    if (m>=0)&(l<self.lens): 
     if(self.string[l]==self.string[m]): 
      count += 2 
     else: 
      break 
   if count>maxcount :       #更新回文子串最长长度 
    maxcount = count 
  return maxcount 
   
 #manacher算法 
 def manacher(self): 
  s = '#'+'#'.join(self.string)+'#'    #字符串处理,用特殊字符隔离字符串,方便处理偶数子串 
  lens = len(s) 
  f = []           #辅助列表:f[i]表示i作中心的最长回文子串的长度 
  maxj = 0          #记录对i右边影响最大的字符位置j 
  maxl = 0          #记录j影响范围的右边界 
  maxd = 0          #记录最长的回文子串长度 
  for i in range(lens):       #遍历字符串 
   if maxl>i:         
    count = min(maxl-i,int(f[2*maxj-i]/2)+1)#这里为了方便后续计算使用count,其表示当前字符到其影响范围的右边界的距离 
   else :          
    count = 1 
   while i-count>=0 and i+count<lens and s[i-count]==s[i+count]:#两边扩展 
    count +=1 
   if(i-1+count)>maxl:       #更新影响范围最大的字符j及其右边界 
     maxl,maxj = i-1+count,i               
   f.append(count*2-1) 
   maxd = max(maxd,f[i])      #更新回文子串最长长度 
  return int((maxd+1)/2)-1      #去除特殊字符

通过上面的程序,使用字符串为长度1000的纯‘a'字符串作为样例,经过测试:

暴力枚举:49.719844s

中心枚举:0.334019s

manacher:0.008000s

由此可见,长度为1000时,暴力枚举的耗时已经无法忍受了,而相比而言,中心枚举在效率上已经有很大幅度的提升,最优的manacher耗时则为更短。

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

Python 相关文章推荐
Python3使用requests登录人人影视网站的方法
May 11 Python
深入理解 Python 中的多线程 新手必看
Nov 20 Python
Python实现计算两个时间之间相差天数的方法
May 10 Python
Python实现的双色球生成功能示例
Dec 18 Python
Python字符串拼接六种方法介绍
Dec 18 Python
学习python的前途 python挣钱
Feb 27 Python
Python 绘制酷炫的三维图步骤详解
Jul 12 Python
Python 函数list&amp;read&amp;seek详解
Aug 28 Python
django从后台返回html代码的实例
Mar 11 Python
Python图像阈值化处理及算法比对实例解析
Jun 19 Python
详解Python描述符的工作原理
Jun 11 Python
端午节将至,用Python爬取粽子数据并可视化,看看网友喜欢哪种粽子吧!
Jun 11 Python
python中字符串的操作方法大全
Jun 03 #Python
Python Logging 日志记录入门学习
Jun 02 #Python
python实现寻找最长回文子序列的方法
Jun 02 #Python
python实现对求解最长回文子串的动态规划算法
Jun 02 #Python
Python 网络爬虫--关于简单的模拟登录实例讲解
Jun 01 #Python
用Python一键搭建Http服务器的方法
Jun 01 #Python
python 编写简单网页服务器的实例
Jun 01 #Python
You might like
探讨:使用XMLSerialize 序列化与反序列化
2013/06/08 PHP
PHP运行SVN命令显示某用户的文件更新记录的代码
2014/01/03 PHP
php之readdir函数用法实例
2014/11/13 PHP
php分页原理 分页代码 分页类制作教程
2016/09/23 PHP
PHP的HTTP客户端Guzzle简单使用方法分析
2019/10/30 PHP
JavaScript延迟加载
2021/03/09 Javascript
JS读取cookies信息(记录用户名)
2012/01/10 Javascript
一个js导致的jquery失效问题的解决方法
2013/11/27 Javascript
浅析js中substring和substr的方法
2015/11/09 Javascript
js简单实现调整网页字体大小的方法
2016/07/23 Javascript
AngularJS 实现弹性盒子布局的方法
2016/08/30 Javascript
js运动事件函数详解
2016/10/21 Javascript
JS实现课堂随机点名和顺序点名
2017/03/09 Javascript
Vue.js展示AJAX数据简单示例讲解
2017/03/29 Javascript
jQuery表单设置值的方法
2017/06/30 jQuery
NodeJS实现图片上传代码(Express)
2017/06/30 NodeJs
Vue动态组件实例解析
2017/08/20 Javascript
Vue数据监听方法watch的使用
2018/03/28 Javascript
详解Vue2.0配置mint-ui踩过的那些坑
2018/04/23 Javascript
解决vue-quill-editor上传内容由于图片是base64的导致字符太长的问题
2018/08/20 Javascript
微信小程序学习笔记之获取位置信息操作图文详解
2019/03/29 Javascript
Python实现的生成格雷码功能示例
2018/01/24 Python
python 将列表中的字符串连接成一个长路径的方法
2018/10/23 Python
对pandas中两种数据类型Series和DataFrame的区别详解
2018/11/12 Python
使用Django2快速开发Web项目的详细步骤
2019/01/06 Python
python简单实现AES加密和解密
2019/03/28 Python
Tensorflow中的图(tf.Graph)和会话(tf.Session)的实现
2020/04/22 Python
涉外经济法专业毕业生推荐信
2013/11/24 职场文书
租房协议书范本
2014/04/09 职场文书
2014机关党员干部“正风肃纪”思想汇报
2014/09/15 职场文书
学习教师敬业奉献模范事迹材料思想汇报
2014/09/19 职场文书
有子女的离婚协议书怎么写(范本)
2014/09/29 职场文书
学校机关党总支领导班子整改工作方案
2014/10/26 职场文书
初中教师个人总结
2015/02/10 职场文书
爱护公物主题班会
2015/08/17 职场文书
SQL Server使用T-SQL语句批处理
2022/05/20 SQL Server