洛谷P4287 [SHOI2011]双倍回文（回文自动机）

传送门

听说有大佬用manacher$O(n)$过此题……太强啦……

说一下PAM的做法吧。（看了题解之后发现）蛮简单的

我们肯定要先建出回文自动机的

然后如果是枚举每一个节点暴跳fail指针肯定得T

那么我们对于每一个节点记录一个$trans[i]$，表示小于等于它长度一半的节点

这个可以在建自动机的时候顺便求出来，具体看代码

然后对每一个节点判断长度是否模4为0且$trans[i]$的长度是它的一半就好了

//minamoto
#include<cstdio>
#include<cstring>
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
const int N=500005;
int fail[N],ch[N][26],cnt[N],len[N],trans[N];
int n,m,tot,last,p,q,ans;
char s[N];
inline int newnode(int x){
    len[++tot]=x;return tot;
}
inline int getfail(int x,int n){
    while(s[n-1-len[x]]!=s[n]) x=fail[x];return x;
}
void build(){
    for(int i=1;s[i];++i){
        int x=s[i]-'a';
        p=getfail(last,i);
        if(!ch[p][x]){
            q=newnode(len[p]+2);
            fail[q]=ch[getfail(fail[p],i)][x];
            ch[p][x]=q;
            if(len[q]<=2) trans[q]=fail[q];
            else{
                int tmp=trans[p];
                while(s[i-1-len[tmp]]!=s[i]||(len[tmp]+2)*2>len[q]) tmp=fail[tmp];
                trans[q]=ch[tmp][x];
            }
        }
        cnt[last=ch[p][x]]++;
    }
}
int main(){
//    freopen("testdata.in","r",stdin);
    scanf("%d",&n);
    scanf("%s",s+1);
    s[0]=-1,fail[0]=1,last=0;
    len[0]=0,len[1]=-1,tot=1;
    build();
    for(int i=tot;i>=2;--i) cnt[fail[i]]+=cnt[i];
    for(int i=2;i<=tot;++i)
    if((len[trans[i]]<<1)==len[i]&&len[i]%4==0) cmax(ans,len[i]);
    printf("%d
",ans);
    return 0;
}