最平常的算法是O(n^3)暴力枚举的。O(n^2)的是从中开始枚举,对奇数串和偶数串分别处理。而manacher采取了增加字符的方式避开了奇偶的讨论,而且达到了理论的下界O(n)的时间复杂度。设mx为字符串能往右延伸的最大点,id为最大点时的对称点。设p数组为当前点能延伸距离,j点为i关于id的对称点,则j=2*id-i,(0<j<i)所以当处理到i时,j点是已经被处理过的了,i点为当前处理点(1)当mx>i时表明最大延伸距离大于当前点,即当前点能向右延伸。但是只能取p[j]和mx-i两者中小的个才合法,即p[i]=min(p[2*id-i],mx-i).(2)mx<=i时p[i]=1; 然后判断之后当前点能延伸的最大距离,即while(s[i+p[i]]==s[i-p[i]]) ++p[i]; 然后处理下mx和id的变化。具体代码如下
#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
using namespace std;
const int maxn=110000+10;
char s[maxn*2],str[maxn];
int p[maxn*2];
int manancher(int len)
{
int l=0;
s[l++]='$';
s[l++]='#';
for(int i=0;i<len;i++)
{
s[l++]=str[i];
s[l++]='#';
}
s[l]=0;
int mx=0,id=0;
for(int i=0;i<l;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(p[i]+i>mx)
{
mx=p[i]+i;
id=i;
}
}
int ans=0;
for(int i=0;i<l;i++)
{
ans=max(ans,p[i]-1);
}
return ans;
}
int main()
{
while(~scanf("%s",str))
{
int len=strlen(str);
printf("%d\n",manancher(len));
}
}