【hdu3068】最长回文——manacher算法

gdkoi时看了一下manacher，发现太弱看不懂……

如今4个月过去了，终于回想起这遗漏的算法，赶紧回来学了。

这道题是manacher裸题，其实只要理解了并不难。

网上有很多博客都写了该算法的原理，这里就不展开了。

说几个要注意的地方：

1.s数组要开大！（我才不会告诉你我第一次就是因为没开两倍RE了呢），因为我们要插入字符（一般是‘#’），所以至少要开2*len+2这么大的数组。

2.读进来的数组长为len，那么我们知道s[len]=''。请注意，这里在插入字符时也要从len开始（即s[len*2+1]='#',s[len*2+2]=''）,这样可以防止多组数据导致越界。

3.在从i开始向左右两边对比字符后，记得更新最右端的位置（即p[pos]+pos也就是改变对称轴pos）。

4.最后答案mx记得减一！

好了就注意这么多，具体细节照惯例看代码。

#include<cstdio>
#include<cstring>
#include<iostream>
const int maxn=220010;
using namespace std;
char s[maxn];
int len,pos;
int p[maxn];
int main()
{
    while(scanf("%s",s)!=EOF)
    {
        memset(p,0,sizeof(p));
        int mx=1;
        len=strlen(s);pos=0;
        for(int i=len;i>=0;i--)
        {
            s[i*2+2]=s[i];
            s[i*2+1]='#';
        }
        s[0]='$';
        for(int i=2;i<len*2+1;i++)
        {
            if(p[pos]+pos>i)p[i]=min(p[pos*2-i],p[pos]+pos-i);
            else p[i]=1;
            while(s[i+p[i]]==s[i-p[i]])p[i]++;
            if(pos+p[pos]<i+p[i])pos=i;
            if(mx<p[i])mx=p[i];
        }
        printf("%d
",mx-1);
    }
    return 0;
}