彻底搞懂马拉车（Manacher）

Manacher算法，又叫“马拉车”，它可以在时间复杂度和空间复杂度都是O(n)的情况下，求出一个字符串的最长回文串长度。

回文串的基本解法

       以每一个点为中心对称点，每次保留最长回文子串的长度，最后得到的就是最长回文子串的长度，但是这样的时间复杂度为O（n^2），并且奇偶回文串是不一样的，比如“aba”的对称中心为“b”，而“abba”的对称中心在两个“b”之间，所以，为了解决这样的问题，可以在每两个字符之间添加一个特殊字符，比如“#”，字符串变成这样“#a#b#b#a#”，求得中间“#”的回文串长度为9，再除2，就得到回文串的长度为4。

    基础解法的时间复杂度太高了，所以就有了Manacher算法，他的复杂度为O(n)，用起来很方便。

    刚开始看马拉车，用了很长时间，只要搞明白原理就能懂了。

下面开始详细讲解Manacher算法

对于奇偶字符串的处理，manacher采用的是填充特殊字符的方法（上文提到过），并且在字符串两端都加入不同的字符，防止越界（因为不同字符一定不是回文串），比如字符串“abbccbba”，增加字符后变成“@#a#b#b#c#c#b#b#a#0”.

重定义字符串代码

重定义字符串的长度为2*n+2。

void getstr() {//重定义字符串
	int k = 0;
	str[k++] = '@';//开头加个特殊字符防止越界
	for (int i = 0; i < len; i++) {
		str[k++] = '#';
		str[k++] = s[i];
	}
	str[k++] = '#';
	len = k;
	str[k] = 0;//字符串尾设置为0，防止越界
}

在重新定义字符串之后，就要开始进行Manacher算法了，先来看一下核心算法。

int manacher() {
	int mx = 0, id;//mx为最右边，id为中心点
	int maxx = 0;
	for (int i = 1; i < len; i++) {
		if (mx > i) Len[i] = min(mx - i, Len[2 * id - i]);//判断当前点超没超过mx
		else Len[i] = 1;//超过了就让他等于1，之后再进行查找
		while (str[i + Len[i]] == str[i - Len[i]]) Len[i]++;//判断当前点是不是最长回文子串，不断的向右扩展
		if (Len[i] + i > mx) {//更新mx
			mx = Len[i] + i;
			id = i;//更新中间点
			maxx = max(maxx, Len[i]);//最长回文字串长度
		}
	}
	return (maxx - 1);
}

看上去不是很好理解，下面开始分步讲解

首先对数据进行初始化，再对i进行判断，分为两种情况
第一种情况，i>=mx,i在mx前面，直接让len[i]=1。
第二种情况，i<mx，这时候就又有两种情况了，对len[j]和mx-i进行比较：
（1）len[j]<=mx-i说明i的最右端还在mx里面，如上图所示，只需要让len[i]=len[j]即可。
（2）len[j]>mx说明i的最右端大于mx了，如下图所示，所以我们需要对这两种情况再讨论一下，当Len[j] < mx-i的时候，表示Len[i]的长度可能不会超过mx-i，所以我们就从i的Len[2*id - i]也就是Len[mx-i]的地方开始匹配。当Len[j] > mx - i的时候，说明i位置的子串长度超过了mx，但mx以外的地方还没有遍历到，所以我们就从mx-i也就是mx的位置开始对i匹配。

如果len[i]+i>mx，就对mx进行更新，并且将中间点id更换成i，再通过比较更新最长回文串长度，返回最大值。

 Manacher模板代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int maxn = 1e6 + 5;
char s[maxn * 2], str[maxn * 2];
int Len[maxn * 2], len;
void getstr() {//重定义字符串
	int k = 0;
	str[k++] = '@';//开头加个特殊字符防止越界
	for (int i = 0; i < len; i++) {
		str[k++] = '#';
		str[k++] = s[i];
	}
	str[k++] = '#';
	len = k;
	str[k] = 0;//字符串尾设置为0，防止越界
}
int manacher() {
	int mx = 0, id;//mx为最右边，id为中心点
	int maxx = 0;
	for (int i = 1; i < len; i++) {
		if (mx > i) Len[i] = min(mx - i, Len[2 * id - i]);//判断当前点超没超过mx
		else Len[i] = 1;//超过了就让他等于1，之后再进行查找
		while (str[i + Len[i]] == str[i - Len[i]]) Len[i]++;//判断当前点是不是最长回文子串，不断的向右扩展
		if (Len[i] + i > mx) {//更新mx
			mx = Len[i] + i;
			id = i;//更新中间点
			maxx = max(maxx, Len[i]);//最长回文字串长度
		}
	}
	return (maxx - 1);
}
int main() {
	scanf("%s", s);
	len = strlen(s);
	getstr();
	printf("%d\n",manacher());
	return 0;
}