最大回文字符串算法Manacher
在刷leetcode时有个求最长回文字符串的问题。
#官方题解提供了4中解决办法,分别是
1.暴力法, 2.动态规划, 3.中心扩展算法, 4.就是我们今天要介绍的Manacher方法。
在介绍算法之前,首先介绍一下什么是回文串,所谓回文串,简单来说就是正着读和反着读都是一样的字符串,比如”aba”,”上海自来水来自海上”等等,一个字符串的最长回文子串即为这个字符串的子串中,是回文串的最长的那个。
而Mnacher是最高效的算法,算法复杂度为O(n)。
Manacher算法原理与实现
Manacher算法提供了一种巧妙地办法,将长度为奇数的和偶数的回文串一起考虑,具体做法是,在原字符串两个字符中间插入一个分隔符#,同时在开头插入(@),尾部也要添加一个分隔符($)。
Len数组的计算
从左往右依次计算Len[i],当计算Len[i]时,Lenj已经计算完毕。设P为之前计算中最长回文子串的右端点的最大值,并且设取得这个最大值的位置为po,分两种情况。
第一种情况:i<=P
找到i相对于po的对称位置,设为j,那么如果Len[j]<P-i
那么说明以j为中心的回文串一定在以po为中心的回文串的内部,且j和i关于位置po对称,由回文串的定义可知,一个回文串反过来还是一个回文串,所以以i为中心的回文串的长度至少和以j为中心的回文串一样,即Len[i]>=Len[j]。因为Len[j]<P-i,所以说i+Len[j]<P。由对称性可知Len[i]=Len[j] 。
如果Len[j]>=P-i,由对称性,说明以i为中心的回文串可能会延伸到P之外,而大于P的部分我们还没有进行匹配,所以要从P+1位置开始一个一个进行匹配,直到发生失配,从而更新P和对应的po以及Len[i]。
第二种情况: i>P
如果i比P还要大,说明对于中点为i的回文串还一点都没有匹配,这个时候,就只能老老实实地一个一个匹配了,匹配完成后要更新P的位置和对应的po以及Len[i]。
Python3实现Manacher算法
def longestPalindrome(self, s: str) -> str:
s = "@#"+"#".join(s)+"#$"
mx,ans,p0=0,0,0 ##### 最右边位置,最大长度,中心位置
Len = [0]*len(s)
mxstr = ''
for i in range(len(s)):
if mx>i:
Len[i]=min(mx-i,Len[2*p0 - i])
else:
Len[i] = 1
while i+Len[i]<len(s) and s[i-Len[i]]==s[i+Len[i]]:
Len[i] += 1
if Len[i]+i > mx :
mx=Len[i]+i
p0 = i
if Len[i] >= ans:
mxstr = s[i-Len[i]+1:i+Len[i]]
ans=Len[i]
return mxstr.replace("#","").replace("$","")
- 原文作者:春江暮客
- 原文链接:https://www.bobobk.com/291.html
- 版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 4.0 国际许可协议进行许可,非商业转载请注明出处(作者,原文链接),商业转载请联系作者获得授权。