题目:给一个字符串,找出它的最长的回文子序列的长度。
例如,如果给定的序列是“BBABCBCAB”,则输出应该是7,“BABCBAB”是在它的最长回文子序列。
输入:
aaaa
1212asdfdsa1144121
输出:
4
7
这里我们还是将其封装成函数调用
何谓回文序列
回文序列就是正向和反向完全一样的序列,比如 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),优势在于避免了奇偶数讨论的问题,简化了边界判断,还记录了当前字符串的“回文状态”,利用之前的回文状态来求当前回文状态 ,体现了动态规划的思想