最长回文子串（Manacher算法）

问题：求最长回文子串，例如abba，最长回文子串为4；abaaaaa，最长回文子串为5。

问题详见：POJ3974

这里介绍的是Manacher算法，其时间复杂度为O(n)。

算法思想：

对于字符串“abba”的两边插入“#”，编程“#a#b#b#a#”。同时为了更好地处理越界问题，我们可以在字符串的开始处加入另一个特殊的字符，如“$”，

则字符串变为“$#a#b#b#a#”，则：

ｓ［ｉ］＝＂ #      ａ     ＃      ｂ      ＃      ｂ      ＃      ａ      ＃＂

ｐ［ｉ］＝＂１     ２      １      ２      ５       2       1        4      １＂

观察，可得max（t [ i ] - 1 ）正是原字符串中最长回文子串的总长度，如上面的最长回文子串的长度为4。

void manacher()
{
	//Manacher算法借助了两个辅助变量id和mx 
	// id表示最大回文子串中心的位置
	//mx为id+p[id]，也就是最大回文子串的边界 
	int id=0,mx=0;
	for(int i=1; t[i]!='\0' ;++i)	{
		//mx>i时，P[i]>= min(p[2*id-i],mx-i)
		p[i]= mx>i ? min(p[2*id-i],mx-i) : 1;
		//计算回文子串的长度	
		while(t[i+p[i]]==t[i-p[i]])
			p[i]++;
		//当 i+p[i] > mx时，
		//更新最大回文子串中心的位置和边界 
		if(i+p[i] > mx){
			mx=i+p[i];
			id=i;
		}
	}
}

完整代码为：

#include <iostream>

using namespace std;

const int maxn=1000010;
char s[maxn],t[2*maxn];
int p[2*maxn];

void manacher()
{
	//Manacher算法借助了两个辅助变量id和mx 
	// id表示最大回文子串中心的位置
	//mx为id+p[id]，也就是最大回文子串的边界 
	int id=0,mx=0;
	for(int i=1; t[i]!='\0' ;++i)	{
		//mx>i时，P[i]>= min(p[2*id-i],mx-i)
		p[i]= mx>i ? min(p[2*id-i],mx-i) : 1;
		//计算回文子串的长度	
		while(t[i+p[i]]==t[i-p[i]])
			p[i]++;
		//当 i+p[i] > mx时，
		//更新最大回文子串中心的位置和边界 
		if(i+p[i] > mx){
			mx=i+p[i];
			id=i;
		}
	}
}

int main()
{
	//添加这个的话，使得cin和cout的效率和scanf差不多 
	//不加的话，在遇到特别大的输入时，会影响效率 
	ios::sync_with_stdio(false);
	int _case=1;
	while(cin>>s && s[0]!='E')	{
		int k=0;
		//第一位置为"$" 
		t[k++]='$';
		//将输入的字符串每隔一位插入#
		//如abbd生成的新字符串为"#a#b#b#d#" 
		for(int i=0;s[i]!='\0';++i){
			t[k++] = '#';
			t[k++] = s[i];	
		}	
		t[k++]='#';
	  	//最后一位为'\0' ，此时该字符串为 "$#a#b#b#d#" 
		t[k]='\0';
		manacher();
		int len=0;
		//遍历字符串，得出字符串的最长回文子串的长度 
		for(int i=1;t[i]!='\0';++i){
			len=max(p[i],len);
		}
		cout << "Case" << _case++ << ":" << len - 1 << endl;
	}
	return 0;
}

综上所述，Manacher算法使用了辅助变量id和mx，可以在每次循环中直接对p [ i ] 快速赋值，这样使得计算回文字串的过程中，不必每次都从1开始比较，

从而减少了比较次数，使得达到线性时间复杂度O(n).