python实现求最长回文子串长度


Posted in Python onJanuary 22, 2018

给定一个字符串,求它最长的回文子串长度,例如输入字符串'35534321',它的最长回文子串是'3553',所以返回4。

最容易想到的办法是枚举出所有的子串,然后一一判断是否为回文串,返回最长的回文子串长度。不用我说,枚举实现的耗时是我们无法忍受的。那么有没有高效查找回文子串的方法呢?答案当然是肯定的,那就是中心扩展法,选择一个元素作为中心,然后向外发散的寻找以该元素为圆心的最大回文子串。但是又出现了新的问题,回文子串的长度即可能是基数,也可能好是偶数,对于长度为偶数的回文子串来说是不存在中心元素的。那是否有一种办法能将奇偶长度的子串归为一类,统一使用中心扩展法呢?它就是manacher算法,在原字符串中插入特殊字符,例如插入#后原字符串变成'#3#5#5#3#4#3#2#1#'。现在我们对新字符串使用中心扩展发即可,中心扩展法得到的半径就是子串的长度。

现在实现思路已经明确了,先转化字符串'35534321'  ---->  '#3#5#5#3#4#3#2#1#',然后求出以每个元素为中心的最长回文子串的长度。以下给出python实现:

#!/usr/bin/python
# -*- coding: utf-8 -*-

def max_substr(string):
  s_list = [s for s in string]
  string = '#' + '#'.join(s_list) + '#'
  max_length = 0
  length = len(string)
  for index in range(0, length):
    r_length = get_length(string, index)
    if max_length < r_length:
      max_length = r_length
  return max_length

def get_length(string, index):
  # 循环求出index为中心的最长回文字串
  length = 0
  r_ = len(string)
  for i in range(1,index+1):
    if index+i < r_ and string[index-i] == string[index+i]:
      length += 1
    else:
      break
  return length

if __name__ == "__main__":
  result = max_substr("35534321")
  print result

功能已经实现了,经过测试也没有bug,但是我们静下心来想一想,目前的解法是否还有优化空间呢?根据目前的解法,我们求出了‘35534321‘中每个元素中心的最大回文子串。当遍历到'4'时,我们已经知道目前最长的回文子串的长度max_length是4,这是我们求出了以4为中心的最长回文子串长度是3,它比max_length要小,所以我们不更新max_length。换句话说,我们计算以4为中心的最长回文字串长度是做了无用功。这就是我们要优化的地方,既然某个元素的最长的回文子串长度并没有超过max_length,我们就没有必要计算它的最长回文子串,在遍历一个新的元素时,我们要优先判断以它为中心的回文子串的长度是否能超越max_length,如果不能超过,就继续遍历下一个元素。以下是优化后的实现:

#!/usr/bin/python
# -*- coding: utf-8 -*-

def max_substr(string):
  s_list = [s for s in string]
  string = '#' + '#'.join(s_list) + '#'
  max_length = 0
  length = len(string)
  for index in range(0, length):
    r_length = get_length2(string, index, max_length)
    if max_length < r_length:
      max_length = r_length
  return max_length

def get_length2(string, index, max_length):
  # 基于已知的最长字串求最长字串
  # 1.中心+最大半径超出字符串范围, return
  r_ = len(string)
  if index + max_length > r_:
    return max_length

  # 2.无法超越最大半径, return
  l_string = string[index - max_length + 1 : index + 1]
  r_string = string[index : index + max_length]
  if l_string != r_string[::-1]:
    return max_length

  # 3.计算新的最大半径
  result = max_length
  for i in range(max_length, r_):
    if index-i >= 0 and index+i < r_ and string[index-i] == string[index+i]:
      result += 1
    else:
      break
  return result - 1

if __name__ == "__main__":
  result = max_substr("35534321")
  print result

那么速度到底提升了多少呢,以字符串1000个‘1'为例,优化前的算法执行时间为0.239018201828,优化后为0.0180191993713,速度提升了10倍左右

/usr/bin/python /Users/hakuippei/PycharmProjects/untitled/the_method_of_programming.py
0.239018201828
0.0180191993713

再给大家分享一个实例:

#!usr/bin/env python
#encoding:utf-8

'''
__Author__:沂水寒城
功能:寻找最长回文子序列
'''

def slice_window(one_str,w=1):
  '''
  滑窗函数
  '''
  res_list=[]
  for i in range(0,len(one_str)-w+1):
    res_list.append(one_str[i:i+w])
  return res_list


def is_huiwen(one_str_list): 
  '''
  输入一个字符串列表,判断是否为回文序列 
  ''' 
  if len(one_str_list)==1: 
    return True  
  else: 
    half=len(one_str_list)/2 
    if len(one_str_list)%2==0: 
      first_list=one_str_list[:half] 
      second_list=one_str_list[half:] 
    else: 
      first_list=one_str_list[:half] 
      second_list=one_str_list[half+1:] 
    if first_list==second_list[::-1]: 
      return True  
    else: 
      return False 


def find_longest_sub_palindrome_str(one_str):
  '''
  主函数,寻找最长回文子序列
  '''
  all_sub=[]
  for i in range(1,len(one_str)):
    all_sub+=slice_window(one_str,i)
  all_sub.append(one_str)
  new_list=[]
  for one in all_sub:
    if is_huiwen(list(one)):
      new_list.append(one)
  new_list.sort(lambda x,y:cmp(len(x),len(y)),reverse=True)
  print new_list[0]


if __name__ == '__main__':
  one_str_list=['uabcdcbaop','abcba','dmfdkgbbfdlg','mnfkabcbadk']
  for one_str in one_str_list:
    find_longest_sub_palindrome_str(one_str)

结果如下:

abcdcba 
abcba 
bb 
abcba 
[Finished in 0.3s]
Python 相关文章推荐
Python下的twisted框架入门指引
Apr 15 Python
python实现带声音的摩斯码翻译实现方法
May 20 Python
Python3.2模拟实现webqq登录
Feb 15 Python
python+opencv实现霍夫变换检测直线
Oct 23 Python
使用Python控制摄像头拍照并发邮件
Apr 23 Python
详解Python二维数组与三维数组切片的方法
Jul 18 Python
django自带serializers序列化返回指定字段的方法
Aug 21 Python
python3的UnicodeDecodeError解决方法
Dec 20 Python
python打包生成so文件的实现
Oct 30 Python
python如何读取.mtx文件
Apr 22 Python
tensorflow中的梯度求解及梯度裁剪操作
May 26 Python
matplotlib画混淆矩阵与正确率曲线的实例代码
Jun 01 Python
Python获取本机所有网卡ip,掩码和广播地址实例代码
Jan 22 #Python
Linux CentOS7下安装python3 的方法
Jan 21 #Python
简述Python2与Python3的不同点
Jan 21 #Python
手把手教你用python抢票回家过年(代码简单)
Jan 21 #Python
分析Python中解析构建数据知识
Jan 20 #Python
学习Python selenium自动化网页抓取器
Jan 20 #Python
python使用pil库实现图片合成实例代码
Jan 20 #Python
You might like
PHP获取文件绝对路径的代码(上一级目录)
2011/05/29 PHP
PHP文件读取功能的应用实例
2015/05/08 PHP
PHP实现导出excel数据的类库用法示例
2016/10/15 PHP
php中通过eval实现字符串格式的计算公式
2017/03/18 PHP
yii框架使用分页的方法分析
2019/07/25 PHP
S2SH整合JQuery+Ajax实现登录验证功能实现代码
2013/01/30 Javascript
如何使用jquery动态加载js,css文件实现代码
2013/04/03 Javascript
JS网页图片按比例自适应缩放实现方法
2014/01/15 Javascript
浏览器窗口大小变化时使用resize事件对框架不起作用的解决方法
2014/05/11 Javascript
上传图片预览JS脚本 Input file图片预览的实现示例
2014/10/23 Javascript
再谈JavaScript线程
2015/07/10 Javascript
js实现千分符和保留几位小数的简单实例
2016/08/01 Javascript
Vue中的v-cloak使用解读
2017/03/27 Javascript
jQuery选择器之表单元素选择器详解
2017/09/19 jQuery
如何为你的JS项目添加智能提示与类型检查详解
2019/03/12 Javascript
JS实现根据详细地址获取经纬度功能示例
2019/04/16 Javascript
JS实现json数组排序操作实例分析
2019/10/28 Javascript
js实现点击生成随机div
2020/01/16 Javascript
javascript实现时间日期的格式化的方法汇总
2020/08/06 Javascript
[01:34]2014DOTA2展望TI 剑指西雅图VG战队专访
2014/06/30 DOTA
Python正则表达式常用函数总结
2017/06/24 Python
python如何为创建大量实例节省内存
2018/03/20 Python
详谈套接字中SO_REUSEPORT和SO_REUSEADDR的区别
2018/04/28 Python
Pandas DataFrame数据的更改、插入新增的列和行的方法
2019/06/25 Python
扩展Django admin的list_filter()可使用范围方法
2019/08/21 Python
使用pandas实现连续数据的离散化处理方式(分箱操作)
2019/11/22 Python
python如何构建mock接口服务
2021/01/28 Python
MONNIER Frères英国官网:源自巴黎女士奢侈品配饰电商平台
2018/12/06 全球购物
Maxpeedingrods美国:高性能汽车零件
2020/02/14 全球购物
财务会计应届生求职信
2013/11/24 职场文书
销售团队激励口号
2014/06/06 职场文书
市场调研项目授权委托书范本
2014/10/04 职场文书
工伤事故处理协议书怎么写
2014/10/15 职场文书
2015年12.4全国法制宣传日活动总结
2015/03/24 职场文书
浅谈MySQL user权限表
2021/06/18 MySQL
Java Socket实现多人聊天系统
2021/07/15 Java/Android