HDU 3068 最长回文(manacher算法:回文字串)

#include <iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
const int MAXN=110000+1000;
char S[MAXN*2];//构造的新串
char T[MAXN];//原始串
int n;
int ans;
int p[MAXN*2];
void init()
{
    int i,j;
    S[0]='@',S[1]='#';
    for(i=0,j=2;T[i];i++,j+=2)
    {
        S[j]=T[i];
        S[j+1]='#';
    }
    S[j]=0;
}
int manacher()
{
    int ans=0;
    int id,i,mx=0;
    for(int i=1;S[i];i++)
    {
        if(mx>i) p[i]=min(p[2*id-i] , mx-i);
        else p[i]=1;
        while(S[i-p[i]] == S[i+p[i]])
            p[i]++;
        if(i+p[i]>mx)
        {
            mx=i+p[i];
            id=i;
        }
        ans = max(ans,p[i]-1);
    }
    return ans;
}
int main()
{
    while(scanf("%s",T)==1)
    {
        init();
        printf("%d\n",manacher());
    }
    return 0;
}