Manacher模板

hdu 3068即一道Manacher的模板题，使用时记得处理每组新数据之前清空辅助数组b[i]（虽然好像这题不memset也不会错）。。。
后来研究发现，对于构造新串的操作，如果用指针一格格往后加入字符，最后指到哪儿长度就取就取多长的话可以不每次清空b[i]数组；如果做n次对b[(i<<1)+1]和b[(i<<1)+2]的修改最后用strlen()取新串长度的话就必须清空，不信可以试试下面这组数据（坏笑ing）。。。
INPUT
aaaaaaaaaaaaaaaaaaaaaa
a
aa

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#include<climits>using namespace std;const int maxn=110004;int n,res,p[maxn<<1];char a[maxn],b[maxn<<1];inline void init() {    memset(p,0,sizeof(p));    n=strlen(a);    for (register int i=0;i<n;i++)        b[(i<<1)+1]='#',b[(i<<1)+2]=a[i];    n=(n<<1)+1,res=0;    b[0]='$',b[n]='#';}inline void manacher() {    int id,mx=0;    for (register int i=1;i<=n;i++) {        p[i] = mx>i ? min(mx-i,p[(id<<1)-i]) : 1;        while (b[i-p[i]]==b[i+p[i]]) p[i]++;        if (mx<p[i]+i) mx=p[i]+i,id=i;        res=max(res,p[i]-1);     }    printf("%d\n",res);}int main() {    while (scanf("%s",a)!=EOF) {        init();        manacher();        memset(b,'\0',sizeof(b));    }    return 0;}