最长回文子串：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;
}