BZOJ 2342 [SHOI2011]双倍回文 (回文自动机)

题目大意：略

先建出$PAM$

因为双倍回文串一定是4的倍数，所以找出$PAM$里所有$dep$能整除4的节点

看这个串是否存在一个回文后缀，长度恰好为它的一半，沿着$pre$链往上跳就行了

暴跳可能会$T$，所以倍增了跳

如果被卡空间，可以把trs数组当成倍增数组

#include <cmath>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N1 500100
#define S1 (N1<<1)
#define ll long long
#define uint unsigned int
#define rint register int 
#define dd double
#define il inline 
#define inf 0x3f3f3f3f
#define idx(X) (X-'a')
using namespace std;

int gint()
{
    int ret=0,fh=1;char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')fh=-1;c=getchar();}
    while(c>='0'&&c<='9'){ret=ret*10+c-'0';c=getchar();}
    return ret*fh;
}
int len;
namespace PAM{
int trs[N1][26],pre[N1],dep[N1],la,tot;
void init(){la=tot=1,pre[0]=pre[1]=1,dep[1]=-1;}
int chk(char *str,int i,int p){return str[i-dep[p]-1]!=str[i]?1:0;}
void insert(char *str,int i)
{
    int p=la,np,fp,c=idx(str[i]);
    while(chk(str,i,p)) p=pre[p];
    if(!trs[p][c])
    {
        np=++tot;
        dep[np]=dep[p]+2;
        fp=pre[p];
        while(chk(str,i,fp)) fp=pre[fp];
        pre[np]=trs[fp][c];
        trs[p][c]=np;
    }
    la=trs[p][c];
}
int multi()
{
    int ans=0,i,j,x;
    trs[1][0]=trs[1][1]=0;
    trs[0][0]=trs[0][1]=0;
    for(i=2;i<=tot;i++) trs[i][0]=i,trs[i][1]=pre[i];
    for(j=2;j<=20;j++)
        for(int i=2;i<=tot;i++)
        trs[i][j]=trs[ trs[i][j-1] ][j-1];
    for(i=2;i<=tot;i++)
        if(dep[i]%4==0)
        {
            x=i;
            for(j=20;j>=0;j--)
                if(dep[trs[x][j]]>=dep[i]/2) x=trs[x][j];
            if(dep[x]!=dep[i]/2) continue;
            ans=max(ans,dep[i]);
        }
    return ans;
}
};
char str[N1];

int main()
{
    //freopen("t2.in","r",stdin);
    //freopen("a.out","w",stdout);
    scanf("%d",&len);
    scanf("%s",str+1);
    PAM::init();
    for(int i=1;i<=len;i++) PAM::insert(str,i);
    printf("%d\n",PAM::multi());
    return 0;
}