PHP实现 Manacher 最大回文子串算法
来源:互联网 发布:mui框架 js 构建dom 编辑:程序博客网 时间:2024/06/01 10:29
题目:给一个字符串,找出它的最长的回文子序列的长度。
例如,如果给定的序列是“BBABCBCAB”,则输出应该是7,“BABCBAB”是在它的最长回文子序列。
输入:aaaa1212asdfdsa1144121
输出:47
这里我们还是将其封装成函数调用
何谓回文序列
回文序列就是正向和反向完全一样的序列,比如 asdfdsa
和 aaaa
接下来我们由浅及深,一步一步来说一下 Manacher
算法,这里我们只说 PHP 的实现
判断回文序列
通过 PHP 很容易实现,只需要判断正向反向是否相同就行了
function huiwen($str){ $str2 = implode(array_reverse(str_split($str)), ""); if ($str == $str2) { echo "$str, yes"; } else { echo "$str, no"; }}
接下来我们就来讲解最大回文子串如何去求
第一版代码
求回文子串,毫无疑问需要每个字符遍历一遍,分别求出来各个字符的回文长度,然后选出最长的那一个
代码如下:
function palindrome($str){ $n = strlen($str); $pos = 0; $max = 0; for ($i = 0; $i < $n; $i++) { for ($j = 0; ($i - $j >= 0) && ($i + $j < $n); $j++) { if ($str[$i - $j] != $str[$i + $j]) { break; } if ($j > $max) { $max = $j; $pos = $i; } } } var_dump(substr($str, $pos - $max, $max * 2 + 1));}
可以看到,这个代码是有bug的,因为回文子串可能是奇数长度,也可能是偶数长度,因为我们是以一个字符为中心来求的,所以用这种方法只能求出奇数长度的回文子串,接下来我们来改进一下
第二版改进
因为我们只能求出奇数长度的回文子串,因此我们需要把字符串改进一下
首先通过在每个字符的两边都插入一个特殊的符号,将所有可能的奇数或偶数长度的回文子串都转换成了奇数长度。比如 abba 变成 #a#b#b#a#, aba变成 #a#b#a#,这样我们再来看一下代码:
function palindrome($str){ $pos = 0; $max = 0; $newStr = "#" . implode(str_split($str), "#") . "#"; $n = strlen($newStr); for ($i = 0; $i < $n; $i++) { for ($j = 0; ($i - $j >= 0) && ($i + $j < $n); $j++) { if ($newStr[$i - $j] != $newStr[$i + $j]) { break; } if ($j > $max) { $max = $j; $pos = $i; } } } $r = substr($newStr, $pos - $max, $max * 2); $res = str_replace("#", "", $r); var_dump($res);}
这个样子我们就写好了基本的回文子串算法了,但是在这个算法中,我们用了两层 for
循环,并且需要判断当前字符的位置是否越界,效率较低,接下来我们来看 Manacher
算法
Manacher 算法实现
首先,为了进一步减少编码的复杂度,可以在字符串的开始和结尾加入另一个特殊字符,这样就不用特殊处理越界问题,这里我们在开头和结尾分别加入 @
和 \0
,如 abba
变成 @#a#b#b#a#\0
接下来,我们引入一个辅助序列 $p[]
来记录各个位置的回文长度(注意:我们这里记录的回文长度是单向的长度,比如 12321
我们记录的回文长度为 3,实际回文长度是 $p[$i] * 2 - 1
)
比如字符串 $s[]
与辅助序列 $p[]
的对应关系如下:
- S # 1 # 2 # 2 # 1 #
- P 1 2 1 2 5 2 1 2 1
最后,核心代码在于这一句: $p[$i] = $mx > $i ? min($p[$j], $mx - $i) : 1;
通过这步操作我们可以避免很多不必要的匹配,我们结合下面的代码来理解这句操作 $mx
是最大回文序列最右侧边界的坐标,$i
是当前要计算的位置,$j
是 $i
相对于最大回文序列中间坐标 $pos
的对称点
由于回文序列的性质,回文序列是对称的,也就是说
如果: $mx > $i
那么 $p[$i] >= $mx > $i ? min($p[$j], $mx - $i) : 1;
可以这么理解:
如果 $mx > $i
,这时 当前位置 在 当前最大回文序列 的右半部分里面,根据回文序列的对称性,可以得出,当前位置 $i
的回文长度一定大于等于与之对称 $j
的回文长度,所以说直接从 $j
的回文长度开始计算
但是如果 $mx > $i
,这时 当前位置 在 当前最大回文序列 之外,无法判断 $mx
以后字符的对称性,因此从 1 开始
function palindrome($str){ // 最大回文序列中间坐标 $pos = 0; // 最大回文长度 $max = 0; // 回文序列最右边界坐标 $mx = 0; $p = array("0" => 1, "1" => 1); $newStr = "@#" . implode(str_split($str), "#") . "#\0"; $n = strlen($newStr); for ($i = 2; $newStr[$i] != "\0"; $i++) { // $i 相对于最大回文序列中间坐标 $pos 的对称点 $j = $pos - $i > 0 ? $pos - $i : 1; $p[$i] = $mx > $i ? min($p[$j], $mx - $i) : 1; while ($newStr[$i - $p[$i]] == $newStr[$i + $p[$i]]) { $p[$i]++; } if ($p[$i] > $max) { $max = $p[$i]; $pos = $i; $mx = $i + $max; } } $r = substr($newStr, $pos - $max + 1, $max * 2 - 1); $res = str_replace(array("#", "@", "\0"), "", $r); var_dump($res);}
Manacher 算法的时间复杂度为O(n),优势在于避免了奇偶数讨论的问题,简化了边界判断,还记录了当前字符串的“回文状态”,利用之前的回文状态来求当前回文状态 ,体现了动态规划的思想
- PHP实现 Manacher 最大回文子串算法
- Manacher算法(最大回文子串),JAVA实现
- Manacher算法。最大回文子串。
- poj3974 manacher算法,最大回文子串
- 最长回文(manacher算法,最大回文子串)
- 线性时间求最大回文子串的Manacher算法
- Manacher——最大回文子串的线性算法
- 求最大回文子串(Manacher算法)
- hdu3294 Girls' research (manacher算法,最大回文子串)
- 刷题基础----Manacher算法求解最大回文子串
- 最长回文子串(manacher算法实现)
- 最长回文子串(manacher算法实现)
- Manacher算法实现求最长回文子串的长度
- Manacher算法求最大回文串
- Manacher算法求最长回文子串
- Manacher算法求最长回文子串
- 最长回文子串(Manacher算法)
- Manacher算法 最长回文子串
- 判断三角形
- Postgresql 开启huge pages
- 清明三天放假归来 jquery的使用
- 集合框架中取出元素的方式for循环增强型与一般形式for循环的对比
- 欢迎使用CSDN-markdown编辑器
- PHP实现 Manacher 最大回文子串算法
- Java enum的用法详解
- 模板匹配
- DEDECMS 列表无限加载
- [转载]【Matlab】 print输出图像大小调整
- 缓存架构设计二三事
- 分布式计算——Daytime协议的实现(TCP版)
- Android 常用 adb 命令总结
- php-fpm 启动参数及重要配置详解