Manacher算法学习

注意事项：
1.manacher只适用于长度为奇数的序列
2.为了回避情况1，我们用其他字符填充原来的字符串让他变成长度为奇数的串
例如：
原串:

aabbcab
acca

变换之后

#a#a#b#b#c#a#b#
#a#c#c#a#

设原串为s，变换后的串为sa

求s的最长回文

int manacher(){
    for(int i=0;i<=(110001<<1);i++)sa[i]=' ';
    n=strlen(s+1);
    l=0;
    for(int i=1;i<=n;i++){
        sa[++l]='#';
        sa[++l]=s[i];
    }
    sa[++l]='#';
    memset(r,0,sizeof(r));
    for(int i=2,j=1,k;i<=l;i+=k){
        while(i+j+1<=l&&i-j-1>0&&sa[i-j-1]==sa[i+j+1])++j;
        r[i]=j;
        for(k=1;k<=r[i]&&r[i-k]!=r[i]-k;k++)r[i+k]=min(r[i-k],r[i]-k);
        j=max(j-k,0);
    }
    int ans=1;
    for(int i=1;i<=l;i++)ans=max(ans,r[i]);
    printf("%d\n",ans);
    return 0;
}

求原串的回文串总数

int manacher(){
    for(int i=0;i<=(100001<<1);i++)sa[i]=' ';
    int n=strlen(s+1);
    l=0;
    for(int i=1;i<=n;i++){
        sa[++l]='#';
        sa[++l]=s[i];
    }
    sa[++l]='#';
    memset(r,0,sizeof(r));
    for(int i=1,maxn=1,now=1;i<=l;i++){
        if(maxn>i)r[i]=min(r[2*now-i],maxn-i);
        else r[i]=1;
        while(i+r[i]<=l&&i-r[i]>0&&sa[i+r[i]]==sa[i-r[i]])++r[i];
        if(i+r[i]>maxn){
            maxn=i+r[i];
            now=i;
        }
    }
    ll ans=0;
    for(int i=1;i<=l;i++)ans+=r[i]/2;
    printf("%lld\n",ans);
    return 0;
}