Manacher算法(O(n)求最大回文字串的长度)(HDU3068)

最近看了一个讲Manacher算法非常好的帖子，地址：http://acm.uestc.edu.cn/bbs/read.php?tid=3258&page=1

非常容易懂，简直太神了。

自己写了下hdu3068，代码如下：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
char s[211111],c[111111];
int p[211111];			//p[i]代表以i为中心，能最多向左或向右扩展多少个字符（包括i）
int main()
{
    //freopen("1.txt","r",stdin);
    //freopen("2.txt","w",stdout);
    while (~scanf("%s",c+1)){
        int len=strlen(c+1);
        for (int i=1;i<=len*2;i++){
            if (i&1) s[i]='#';
            else s[i]=c[i/2];
        }
        s[0]='@';       //防止越界
        s[len*2+1]='!';
        int mx,id;		//mx为前面能匹配到的最右边的位置+1
        mx=0;		
        for (int i=1;i<=len*2;i++){
            if (mx>i) p[i]=min(p[id*2-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];
                id=i;
            }
        }
        int ans=0;
        for (int i=1;i<=len*2;i++)
            if (s[i]=='#') ans=max(ans,p[i]/2*2);
            else ans=max(ans,(p[i]-1)/2*2+1);
        printf("%d\n",ans);
    }
    return 0;
}