马拉车算法
manacher算法又称马拉车算法,是用来求解字符串中的最长连续回文子串的 方法,可以将时间复杂度降低到O(n)。马拉车算法的核心思想其实是中心扩展思想的变形,首先先来解释一下常见的求最长回文子串的方法,中心扩展思想 。
中心扩展思想
首先对于要求解的字符串ababaabc 来说,因为存在奇数长度的回文串和偶数长度的回文串,所以我们可以以一个字符为中心,进行左右同时扩展,这样得到的长度是奇数长度的回文串,同时我们也可以以两个字符串为中心进行扩展,这样得到的就是偶数长度的回文串。
这样代码就很容易实现:
public String getLongestPalindrome(String s){
if(s.length() == 0)
return "";
int beg = 0,end = s.length();
for(int i = 0;i <s.length();i++){
int l1 = help(s,i,i); //以一个字符为中心进行左右扩展
int l2 = help(s,i,i+1); //以两个字符为中心进行左右扩展
int maxL = Math.max(l1,l2);
if(maxL > end-beg){ // 如果当前最大长度大,进行更新
beg = i-(maxL-1)/2;
end = i+maxL/2;
}
}
return s.subString(beg,end+1);
}
public int help(String s,int beg,int end){ //实现中心扩展思想的算法代码
if(end >= s.length())
return 1;
int i = beg,j = end;
int cnt = 1;
while(i >=0 && j < s.length() && s.charAt(i) == s.charAt(j)){
cnt++;
i--;
j++;
}
return cnt;
}
manancher算法
说完中心扩展思想,我们来进入马拉车算法的解释,马拉车算法是更快速的求解最长连续回文子串的方法,具体的实现方式是对字符串进行变形,在每个字符串之间插入 # 字符,同时为了保证计算最长子串时的索引值的正确,还需要在开始和结束时插入两个不同的符合,例如可以插入 ^ $.
- 在两端分别插入 “^” 和 “$”,两个不可能在字符串中出现的字符,这样中心扩展的时候,判断两端字符是否相等的时候,如果到了边界就一定会不相等,从而出了循环。经过处理,字符串的长度永远都是奇数了。
对于字符串cbcbccde 来说,更新后的字符串应当为:
其中,P数组存放的是从当前位置进行中心扩展后的回文长度半径,
- 回文长度半径: 是指以当前位置为中心的回文长度的一半
T是进行处理之后的字符串。我们可以看出,这样保证了奇数长度的回文和偶数长度的回文都可以依次遍历。
马拉车算法的核心就是通过镜像的方法更新P[i]
- C,记录最大的P[i]的索引值
- R,所能到达的最右边界,即P[C]+C
对于当前的位置 i来说,如果 i< R,说明当前位置是以C为中心的回文串的一部分,根据对称性,存在i_mirror = 2*C-i,此时对于 i 来说,P[i]的值应当为(P[i_mirror],R-i)的最小值,为什么会这样呢:
首先我们知道,在C-P[C] 到R 的返回为,字符串关于C对称,此时i < R,i_mirror >C-P[C],存在两种情况:
- 以i_mirror为中心的回文串是以C为中心的回文串的一部分,所以此时P[i_mirror] < R-i
- 以i_mirror为中心的回文串包含C为中心回文串的一部分,所以此时镜像的法则是,以i为中心的回文串必然包含那一部分,即R-i 的这一部分
所以不管怎样,当前P[i]的初始值必然为Math.min(P[i_mirror],R-i);
我们以上面的图像为例,首先
对于位置索引值为6的地方来说,P[6] = 5,C = 6,R = 11,回文串的最右边界是R = 11,
而当索引值更新到 i = 8时,此时 i < R,存在镜像i_mirror = 4,P[i_mirror] = 3,P[i_mirror] <= R-I,P[i] 的也更新为3
上面说的可能很难理解,都是自己的一些浅显的解释,毕竟我也是刚刚感觉自己搞懂了,我认为结合代码和图像还是比较容易理解的,下面给出java的实现代码
public String preProcess(String s){
if(s.length() == 0)
return "^$";
String str = "";
str += "^#";
for(int i = 0;i < s.length();i++){
str += s.charAt(i);
str += "#";
}
str += "$";
return str;
}
public String manacher(Stirng s){
String T = preProcess(s);
int len = T.length();
int[] p = new int[len];
int C = 0,R = 0;
for(int i = 0;i < len;i++){
int i_mirror = 2*C-i;
if(i < R){
p[i] = Math.min(p[i_mirror],R-i);
}
else
p[i] = 0;
//根据中心扩展思想进行更新p[i]
while(T.charAt(i+P[i]+1) == T.charAt(i-1-P[i]))
p[i]++;
//更新 C和R
if(i+p[i] >R){
C = i;
R = i+p[i];
}
}
}
参考资料:
怎么理解马拉车算法知乎