poj3947最长回文串

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;

typedef long long LL;
const int maxn=1000005;
char s[maxn],t[maxn<<1];
int Len[maxn<<1];

int init(char *st)
{
int i,len=strlen(st);
t[0]='@';//字符串开头增加一个特殊字符，防止越界
for(i=1;i<=2*len;i+=2)
{
t[i]='#';
t[i+1]=st[i/2];
}
t[2*len+1]='#';
t[2*len+2]='$';//字符串结尾加一个字符，防止越界
t[2*len+3]=0;
return 2*len+1;//返回转换字符串的长度
}
//Manacher算法计算过程
int manacher(char *st,int len)
{
int p=0,ans=0,po=0;//mx即为当前计算回文串最右边字符的最大值
for(int i=1;i<=len;i++)
{
if(p>i)
Len[i]=min(p-i,Len[2*po-i]);//在Len[j]和mx-i中取个小
else
Len[i]=1;//如果i>=mx，要从头开始匹配
while(st[i-Len[i]]==st[i+Len[i]])
Len[i]++;
if(Len[i]+i>p)//若新计算的回文串右端点位置大于mx，要更新po和mx的值
{
p=Len[i]+i;
po=i;
}
ans=max(ans,Len[i]);
}
return ans-1;//返回Len[i]中的最大值-1即为原串的最长回文子串额长度
}
int main()
{
int k=1;
while(scanf("%s",s)!=-1)
{
if(strcmp(s,"END")==0)
break;
init(s);
int len=init(s);
int ans=manacher(t,len);
printf("Case %d: %d\n",k++,ans);
}
return 0;
}