【Manacher】最长回文子串

caioj任意门hz2016评测传送门可以的话来一下hz2016评测吧，有的题caioj没有的我也可以给到数据嘛。
#include<map>#include<queue>#include<cmath>#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>#define Maxchar 100000#define mes(x,y) memset(x,y,sizeof(x));#define mpy(x,y) memcpy(x,y,sizeof(x))#define INF 2147483647using namespace std;int p[2*Maxchar+1];char s[2*Maxchar+1],now[2*Maxchar+1];int Manacher(){    int len=strlen(s+1);    for(int i=1;i<=len;i++)now[2*i-1]='#',now[2*i]=s[i];    len=len*2+1;    now[len]='#';    int pos=0,R=0;//R表示当前访问到的所有回文子串，所能触及的最右一个字符的位置,pos表示R位置的回文串的中心    for(int i=1;i<=len;i++){        int j=2*pos-i;//j为i关于pos的对称点         if(i<=R)p[i]=min(p[j],R-i);/*如果i<=R的话，分两种情况第一种情况p[j]>R-i时，表示j对称点的最长回文子串已经越出R的界限了这时，因为我们不确定大于R时的情况，所以p[i]暂时等于R-i第二种情况p[j]<=R-i时，那就可以直接继承p[j]得p[i]=p[j]综合以上两种情况，我们可以用p[i]=min(p[j],R-i)来归纳*/        else p[i]=1;//不然就要暴力枚举        while(1<=i-p[i]&&i+p[i]<=len&&now[i-p[i]]==now[i+p[i]])p[i]++;//向左右枚举，更新p[i]，即是作为继承后的补充，又是作为暴力枚举的操作        if(i+p[i]-1>R){pos=i;R=i+p[i]-1;}//适时更新pos和R    }    int ans=0;    for(int i=1;i<=len;i++)ans=max(ans,p[i]-1);/*首先当前的最长回文子串长度为2*p[i]-1因为我们得到的p数组是在加了#号后的字符串上操作的，所以我们要对答案进行处理因为#号处于首尾和每个字符之间，所以我们就可以保证所得出的最长回文子串的首尾都为#这时我们可以得出不带#号的回文串的长度为(2*p[i]-1-1)/2=p[i]-1所以真正的最长回文子串就是p[i]-1ans记录最长的回文子串长度*/    return ans;}int main(){    while(scanf("%s",s+1)!=EOF){        printf("%d\n",Manacher());    }    return 0;}
Manacher就是kmp嘛，不就是把原来的那个向右匹配变成向左匹配，程序里面右注释，读注释去，很好理解的。

查看原文：http://hz2016.tk/blog/?p=23