设p[i]为当前节点i最长回文串的右节点与i的距离,i为回文中心,而有些回文串是以空格为回文中心的,比如bb,那么为了使
每个回文串都是以字符为中心,把每个字符的前后都加上一个#(操作上就是后面加上#,最后在最前面加#)。
设maxp为最长的长度,p0为maxp的回文中心。那么更新p[i]的时候有两中情况。如果i<maxp时,分两种
一种是可以直接到达maxp的位置也就是p0+p[p0]-i的长度,还有一种找到i关于p0对称的位置2*p0-i,
因为那个位置的长度和i关于p0对称,所以p[2*p0-i]也是p[i]的解,但是这种情况在大于maxp时有不确定性,
所以取最小值。但为了保证答案的准确性和考虑在i>maxp时的情况,暴力验证左右是否满足回文。并时时更新maxp和p0。
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#define N 51000005
using namespace std;
int p[N],l,maxp,p0,ans;
char ss[N*2],a[N];
void manacher()
{
for(int i=1;i<l;i++)
{
if(i<maxp)
{
p[i]=min(p[2*p0-i],p[p0]+p0-i);
}else
{
p[i]=1;
}
while(ss[i-p[i]]==ss[i+p[i]])
{
p[i]++;
}
if(i+p[i]>maxp)
{
maxp=i+p[i];p0=i;
}
}
}
int main()
{
scanf("%s",a+1);
l=strlen(a+1);
for(int i=1;i<=l;i++)
{
ss[i*2]=a[i];
ss[i*2+1]='#';
}
l=l*2+2;
ss[0]='^',ss[1]='#',ss[l]='0';
manacher();
for(int i=1;i<=l;i++)
{
ans=max(ans,p[i]-1);
}
printf("%d\n",ans);
return 0;
}