POJ 3974 回文串-Manacher

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
using namespace std; 
const int MAXN=1000000+5;
typedef long long int LL;
#define INF 0x3f3f3f3f
char str[MAXN],dstr[MAXN*3]; 
int lenstr,lendstr,p[MAXN*3],ans;
void manacher(){
    memset(p,0,sizeof(p));
    int id=0,mx=0;
    for(int i=1;i<lendstr;i++){
        if(mx>i){
            p[i]=min(p[2*id-i],mx-i);
        }
        else{
            p[i]=1;
        }
        while(dstr[i-p[i]]==dstr[i+p[i]]){
            p[i]++;
        }
        if(p[i]+i>mx){
            mx=p[i]+i;
            id=i;
        }
    }
}
void init(){
    dstr[0]=‘$‘;
    dstr[1]=‘#‘;
    for(int i=0;i<lenstr;i++){
        dstr[i*2+2]=str[i];
        dstr[i*2+3]=‘#‘;
    }
    lendstr=lenstr*2+2;
    dstr[lendstr]=‘*‘;
}
int main()
{
    int Ca=1;
    while(~scanf("%s",str)){
        if(strcmp(str,"END")==0){
            break;
        }
        lenstr=strlen(str); 
        init();
        manacher();
        ans=0;
        for(int i=0;i<lendstr;i++){
            ans=max(ans,p[i]);
        }
        printf("Case %d: %d\n",Ca++,ans-1);
    }
    return 0;
}