最长回文子串：Manacher算法

详解见：

http://blog.csdn.net/xingyeyongheng/article/details/9310555

解析：

首先预处理字符串，在字符串首尾和每两个字符之间插入特殊符号（例如’#’）使字符串变成奇数长度，为了进一步减少编码复杂度，还可以再字符串开始再加入另一个特殊字符，这样就不用特殊处理越界问题了。

然后利用计算字符串的P[i]数组来记录以字符S[i]为中心的最长回文子串向左/右扩张的长度即回文半径。

eg：

S：  #  1  #  2  #  2  #  1  #  2  #  3  #  2  #  1  #
P：  1   2  1  2  5   2  1  4   1  2  1  6   1  2   1  2  1
(p.s. 可以看出，P[i]-1正好是原字符串中以S[i]为中心形成回文串的总长度）

计算半径数组：

void manacher()//p[]为半径数组，str[]为处理过后字符串
{
int id, mx, i;
id = mx = 0;
for (i=1; i<=l; ++i)
{
if (mx > i) p[i] = min(p[2*id-i], mx-i+1);
else p[i] = 1;
while (str[i+p[i]]==str[i-p[i]]) ++p[i];
if (p[i]+i-1>mx)
{
id = i;
mx = i + p[i] - 1;
}
}
}

最后遍历一次半径数组，值最大的就是最长回文子串。

练习：

http://acm.hdu.edu.cn/showproblem.php?pid=3068

分析&题解：

直接用Manacher算法求最长回文子串，不知为何之前自己敲的一直TLE，感觉差不多，后来找了个AC代码，以后直接套模板吧 (= =)

AC 代码：

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

const int nn = 220050;
char str[nn], str0[nn];
int len, p[nn];

void manacher()
{
int id, mx, i;
id = mx = 0;
for (i=1; i<=len; ++i)
{
if (mx > i) p[i] = min(p[2*id-i], mx-i+1);
else p[i] = 1;
while (str[i+p[i]]==str[i-p[i]]) ++p[i];
if (p[i]+i-1>mx)
{
id = i;
mx = i + p[i] - 1;
}
}
}

int main(){

while (scanf("%s", str0) != EOF)
{
len = strlen(str0);
str[0] = '$'; str[1] = '#';
for (int i=0; i<len; ++i)
{
str[2*i+2] = str0[i];
str[2*i+3] = '#';//预处理数组
}
len = 2*len+1;
str[len+1] = '\0';

manacher();
int ans = 0;
for (int i=1; i<len; ++i)
if (ans < p[i]-1) ans = p[i]-1;
printf("%d\n", ans );
}
return 0;
}