以下是manacher算法的板子
#include<bits/stdc++.h>
#define N 30000005
#define R register int
using namespace std;
char s[N];
int p[N],ans;
int main()
{
char ch=getchar();s[0]='@';int len=1;s[1]='#';
while (ch>='a' && ch<='z') s[++len]=ch,s[++len]='#',ch=getchar();s[len+1]='*';
int mx=0,mid=0;
for (R i=1;i<=len;++i)
{
if (i<mx) p[i]=min(p[(mid<<1)-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],mid=i;ans=max(ans,p[i]-1);
}
printf("%d\n",ans);
return 0;
}
以字符串aabbccbbd为例
跑出来的回文半径如下:
∗a∗a∗b∗b∗c∗c∗b∗b∗d∗
1232123212721232121
也就是说当前这个回文中心的位置也会产生一个贡献,那代码第12行的mx−i不应该是mx−i+1吗。以及第16行的i+p[i]不应该是i+p[i]−1吗?(虽然跑出来的结果是一样的,因为往后也会更新到,但是这样的打法本质上是有问题的?)