HDU 3068 最长回文 manacher 算法，基本上是O（n）复杂度

下面有别人的比较详细的解题报告：

http://wenku.baidu.com/view/3031d2d3360cba1aa811da42.html

下面贴我的代码，注释在代码中：

#include <cstdio>
#include <cstring>
#define N 110005
char a[N];
char s[N<<1];
int p[N<<1];
int min(int a,int b)
{
    return a<b?a:b;
}
int main()
{
//    freopen("out.cpp","r",stdin);
    s[0] = '#';
    while(scanf("%s",a) != EOF)
    {
        int len = strlen(a);
        int k=1;
        //在原字符串两个字符之间插入字符'*',左右两边也各插一个
        //这样求回文时只需要求奇数个数的回文
        for(int i=0; i<len; ++i)
        {
            s[k++] = '*';
            s[k++] = a[i];
        }
        s[k++] = '*';
        int n = 2*len +1;
        int max =0; //当以某个字符为中心求回文时延伸到右边最远的位置的序号
        int d;//上面最远位置对应的中心位置
        int ans =0;//记录最长回文串的长度
        for(int i=1; i <= n; ++i)
        {
            if(max > i)
                p[i] = min(max-i,p[2*d-i]); //一条线段，以位置d为中心，左边是j
                //右边是i,这样有i+j = 2*d , j = 2*d-i;这里是核心，避免了重复计算，类似记忆化搜索
            else
                p[i] = 1;
            while( s[i-p[i]] == s[i+p[i]] )  ++p[i];//小技巧，给字符串的最前面添加一个特殊字符，避免数组越界，而不是去判断
            if(i+p[i] > max)
            {
                max = i+p[i];
                d = i;
            }
            if(p[i] > ans) ans = p[i];
        }
        printf("%d\n",ans-1);
    }
    return 0;
}