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=,fh=;char c=getchar();
     while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
     while(c>=''&&c<=''){ret=ret*+c-'';c=getchar();}
     return ret*fh;
 }
 int len;
 namespace PAM{
 int trs[N1][],pre[N1],dep[N1],la,tot;
 void init(){la=tot=,pre[]=pre[]=,dep[]=-;}
 int chk(char *str,int i,int p){return str[i-dep[p]-]!=str[i]?:;}
 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]+;
         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=,i,j,x;
     trs[][]=trs[][]=;
     trs[][]=trs[][]=;
     for(i=;i<=tot;i++) trs[i][]=i,trs[i][]=pre[i];
     for(j=;j<=;j++)
         for(int i=;i<=tot;i++)
         trs[i][j]=trs[ trs[i][j-] ][j-];
     for(i=;i<=tot;i++)
         if(dep[i]%==)
         {
             x=i;
             for(j=;j>=;j--)
                 if(dep[trs[x][j]]>=dep[i]/) x=trs[x][j];
             if(dep[x]!=dep[i]/) 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+);
     PAM::init();
     for(int i=;i<=len;i++) PAM::insert(str,i);
     printf("%d\n",PAM::multi());
     return ;
 }