Manacher求最长回文子串

经典例题：POJ 3974

爱吃代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<string>
#include<cmath>
#include<cstring>
#include<queue>
#define inf 0x3f3f3f3f
const double PI=acos(-1.0);
using namespace std;
char str[1000010],strnew[2000020];
int p[2000020],len;
void init()
{
    int i;
    len=strlen(str);
    strnew[0]='@';
    strnew[1]='#';
    for(i=0;i<len;i++)
    {
        strnew[i*2+2]=str[i];
        strnew[i*2+3]='#';
    }
    strnew[len*2+2]='$';
}
int Manacher()
{
    memset(p,0,sizeof(p));
    int i,mx=0,di=0,ans=0;
    for(i=1;i<len*2+2;i++)
    {
        if(i<mx)
            p[i]=min(p[2*di-i],mx-i);
        else
            p[i]=1;
        for(;strnew[i-p[i]]==strnew[i+p[i]];p[i]++);
        if(p[i]+i>mx)
        {
            mx=p[i]+i;
            di=i;
        }
        ans=max(ans,p[i]);
    }
    return ans;
}
int main()
{
    int ans=0,cnt=1;
    while(scanf("%s",str)!=EOF)
    {
        if(strcmp(str,"END")==0)
            break;
        init();
        ans=Manacher();
        printf("Case %d: %d\n",cnt++,--ans);
    }
    return 0;
}