manacher算法原理和代码实现

马拉车算法

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];
		}
	}
}

参考资料：
怎么理解马拉车算法知乎