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;
}