最长回文（马拉车算法模板解析）

最长回文

给出一个只由小写英文字符a,b,c…y,z组成的字符串S,求S中最长回文串的长度.
回文就是正反读都是一样的字符串,如aba, abba等
Input
输入有多组case,不超过120组,每组输入为一行小写英文字符a,b,c…y,z组成的字符串S
两组case之间由空行隔开(该空行不用处理)
字符串长度len <= 110000
Output
每一行一个整数x,对应一组case,表示该组case的字符串中所包含的最长回文长度.
Sample Input

aaaa

abab

Sample Output

4
3

解析：
什么是马拉车算法：
比较官方的介绍，Manacher（马拉车）算法是查找一个字符串的最长回文子串的线性算法。
我的理解就是从字符串的开始，到末尾结束，找到当前右边界最大的回文串子串，记录对称中心id和右边界mx，利用对称性，给下一个以i为对称中心的回文子串赋个初值p[i]，以减少时间，之后再以i为对称中心，在p[i]的基础上向两边遍历，找到以i为中心的最大回文子串的半径p[i]，之后更新id和mx。

字符串的重构：
由于回文串分为奇回文串和偶回文串，为避免这种差异，可将一个初始的字符串用一个不会出现的字符隔开，例如#a#b#a#，#a#b#b#a#，此时重构的回文串都是奇回文串，在重构的字符串前面加一个不同于间隔字符并且原字符串中不会出现的字符，使字符串从下标1开始。

int l=0,i;
a[l++]='$';//使字符串从1开始 
a[l++]='#';//重构数组，使其都变成奇数串 
for(i=0;i<len;i++)
{
	a[l++]=s[i];
	a[l++]='#';
}
a[l]=0;//清空末尾 

核心代码的讲解：
用p[i]数组表示以i为回文中心时最长回文串的半径
从重构的字符串开头开始遍历，首先找到当前右边界最大的回文串子串，记录对称中心id和右边界mx，假设左边界为mn，则对于下一个i（此时i在id右面）有两种情况，一种是对于以id为对称中心，i的对称点j=2*id-i，i+p[j]<mx，即以j为中心的回文串在mn和mx之间，此时利用以id为中心的回文串的对称性，可知p[i]=p[j]，此时p[j]<mx-i;另一种就是i+p[j]>mx，即以j为中心的最大回文串左边界小于mn，此时再利用以id为中心的最大回文串的对称性只能得到i到mx之间的这一段是回文子串，所以此时p[i]=mx-i,(p[j]>mx-i)。之后在以上初值的基础上再判断a[i+p[i]]与a[i-p[i]是否相等，若相等则p[i]++，一直判断，知道不相等为止，最后更新id和mx即可。

int mx=0,id=0;
for(i=0;i<l;i++)
{
	if(i<mx)
	p[i]=min(p[2*id-i],mx-i);//2*id-i是以id为对称中心i的对称点 
	else
	p[i]=1;//半径就是啊a[i]本身，所以为1 
	while(a[i+p[i]]==a[i-p[i]])//在原基础上向两边拓展 
	p[i]++;
	if(i+p[i]>mx)
	{
		mx=i+p[i];//mx为以id为中心的最长回文子串的右边界 
		id=i;
	}
}

完整代码如下：

#include"stdio.h"
#include"string.h"
#include"algorithm"
using namespace std;
char a[200010],s[100010];
int p[200010];//p[i]为以i为回文中心时最长回文串的半径 
void Manacher(int len)
{
	int l=0,i;
	a[l++]='$';//使字符串从1开始 
	a[l++]='#';//重构数组，使其都变成奇数串 
	for(i=0;i<len;i++)
	{
		a[l++]=s[i];
		a[l++]='#';
	}
	a[l]=0;//清空末尾 
	int mx=0,id=0;
	for(i=0;i<l;i++)
	{
		if(i<mx)
		p[i]=min(p[2*id-i],mx-i);//2*id-i是以id为对称中心i的对称点 
		else
		p[i]=1;//半径就是啊a[i]本身，所以为1 
		while(a[i+p[i]]==a[i-p[i]])//在原基础上向两边拓展 
		p[i]++;
		if(i+p[i]>mx)
		{
			mx=i+p[i];//mx为以id为中心的最长回文子串的右边界 
			id=i;
		}
	}
}
int main()
{
	while(~scanf("%s",s))
	{
		int len=strlen(s),i;
		Manacher(len);
		int ans=0;
		for(i=0;i<2*len+2;i++)//依次遍历，找最长的半径，半径减一就是原字符串中回文串的总长度
		ans=max(ans,p[i]-1);
		printf("%d\n",ans);
	}
	return 0;
}