BZOJ 3676: [Apio2014]回文串 后缀自动机 Manacher 倍增

http://www.lydsy.com/JudgeOnline/problem.php?id=3676

过程很艰难了，第一次提交Manacher忘了更新p数组，超时，第二次是倍增的第0维直接在自动机里完成，但是忽略了增加新点时fa变动的情况，还是肉眼查错最管用。

得到的教训是既然倍增就在倍增的函数里完成，自动机就在自动机里完成，不要随便乱搞赋值。

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<map>
 using namespace std;
 const int maxn=;
 char ch[maxn]={};
 char ch1[maxn*]={};
 int siz,siz1;
 struct sam{
     int sig[];
     int f,len,d[],num;
 }t[maxn*];
 int tot=,la=;
 int loc[maxn]={},cnt[maxn*]={},b[maxn*]={},p[maxn*]={};
 long long ans=;
 void add(int z){
     int x=++tot,i=la;
     t[x].len=t[la].len+;
     for(;i&&!t[i].sig[z];i=t[i].f)
         t[i].sig[z]=x;
     if(!i)t[x].f=;
     else{
         int p=t[i].sig[z];
         if(t[p].len==t[i].len+)t[x].f=p;
         else{
             int y=++tot;
             t[y]=t[p];t[y].len=t[i].len+;
             t[p].f=t[x].f=y;
             for(;i&&t[i].sig[z]==p;i=t[i].f)
                 t[i].sig[z]=y;
         }
     }
     la=x;
 }
 void pre(){
     int z;
     for(int i=;i<=tot;i++)cnt[t[i].len]++;
     for(int i=;i<=siz;i++)cnt[i]+=cnt[i-];
     for(int i=tot;i;i--)b[cnt[t[i].len]--]=i;
     for(int i=;i<=siz;i++)t[loc[i]].num=;
     for(int i=tot;i;i--){
         z=b[i];t[t[z].f].num+=t[z].num;
     }
     for(int i=;i<=tot;i++){
         z=b[i];t[z].d[]=t[z].f;
         for(int j=;j<=;j++)
             t[z].d[j]=t[t[z].d[j-]].d[j-];
     }
 }
 void getit(int l,int r){
     r=r/;l=(l+)/;
     if(r<l)return;
     int z=loc[r];
     for(int i=;i>-;i--){
         if(t[t[z].d[i]].len>=r-l+){
             z=t[z].d[i];
         }
     }
     long long sum=(long long)t[z].num*(r-l+);
     if(ans<sum)ans=sum;
 }
 void Manacher(){
     for(int i=;i<=siz;i++)ch1[i*]=ch[i],ch1[i*-]='#';
     ch1[siz*+]='#';ch1[siz*+]='\0';ch1[]='$';
     siz1=siz*+;
     int mx=,id=;
     for(int i=;i<=siz1;i++){
         if(mx>=i)p[i]=min(mx-i,p[id*-i]);
         else p[i]=;
         while(ch1[i+p[i]]==ch1[i-p[i]]){++p[i];getit(i-p[i]+,i+p[i]-);}
         if(mx<i+p[i]){mx=i+p[i];id=i;}
     }
 }
 int main(){
     memset(t,,sizeof(t));
     scanf("%s",ch+);siz=strlen(ch+);
     for(int i=;i<=siz;i++){loc[i]=tot+;add(ch[i]-'a');}
     pre();
     Manacher();
     printf("%lld\n",ans);
     return ;
 }