[bzoj3676]回文串[后缀数组+Manacher]

后缀数组+Manacher

 #include <iostream>
 #include <cstdio>
 #include <cstdlib>
 #include <cstring>
 #include <cmath>
 #include <ctime>
 #include <algorithm>

 using namespace std;

 long long   Ans;
 int A[],B[],U[];
 int SA[],*Rank,H[],Tmp[];
 int lg2[],ST[][],p[];
 char    str[];

 void    Get_H(const int n)
 {
     int i,j,k=;
     for(i=;i<n;H[Rank[i++]]=k)
         for(k?k--:k,j=SA[Rank[i]-];str[i+k]==str[j+k];++k);
     for(i=;i<=n;++i)lg2[i]=lg2[i>>]+;
     for(i=;i<=n;++i)ST[i][]=H[i];
     for(j=;(<<j)<=n;++j)
     {
         for(i=;i+(<<j)-<=n;++i)
         {
             ST[i][j]=min(ST[i][j-],ST[i+(<<(j-))][j-]);
         }
     }
     return ;
 }

 int Query(const int l,const int r)
 {
     int temp=lg2[r-l+];
     return min(ST[l][temp],ST[r-(<<temp)+][temp]);
 }

 bool    cmp(const int * s,const int a,const int b,const int l)
 {
     return s[a]==s[b] && s[a+l]==s[b+l];
 }

 int*    Get_SA(const int n,int  m)
 {
     int i,j,_p,*x=A,*y=B;
     for(i=;i<m;++i)U[i]=;
     for(i=;i<n;++i)U[x[i]=str[i]]++;
     for(i=;i<m;++i)U[i]+=U[i-];
     for(i=n-;i>=;--i)SA[--U[x[i]]]=i;

     for(j=,_p=;_p<n;m=_p,j<<=)
     {
         for(_p=,i=n-j;i<n;++i)y[_p++]=i;
         for(i=;i<n;++i)if(SA[i]>=j)y[_p++]=SA[i]-j;
         for(i=;i<n;++i)Tmp[i]=x[y[i]];
         for(i=;i<m;++i)U[i]=;
         for(i=;i<n;++i)U[Tmp[i]]++;
         for(i=;i<m;++i)U[i]+=U[i-];
         for(i=n-;i>=;--i)SA[--U[Tmp[i]]]=y[i];
         for(swap(x,y),_p=,x[SA[]]=,i=;i<n;++i)
             x[SA[i]]=cmp(y,SA[i-],SA[i],j)?_p-:_p++;
     }
     return x;
 }

 long long   Calc(int l,int r,const int n)
 {
     l=(l-)>>,r=(r-)>>;
     int pos=Rank[l],L,R,temp=;

     L=,R=pos;

     while(L<R-)
     {
         int mid=L+((R-L)>>);
         if(Query(mid+,pos)>=r-l+)R=mid;
         else    L=mid;
     }

     temp+=pos-R;
     L=pos,R=n+;

     while(L<R-)
     {
         int mid=L+((R-L)>>);
         if(Query(pos+,mid)>=r-l+)L=mid;
         else    R=mid;
     }
     temp+=L-pos;
     return (long long)temp*(r-l+);
 }

 void    Manacher(const int n)
 {
     int i,pos=;
     for(i=n-;i>=;--i)
     {
         str[i+i+]=str[i];
         str[i+i+]='#';
     }
     str[]='^';str[n<<|]='#';str[(n+)<<]='$';
     for(i=;i<=(n<<|);++i)
     {
         if(p[pos]+pos>i)
             p[i]=min(p[(pos<<)-i],p[pos]+pos-i);
         else    p[i]=;
         while(str[i-p[i]]==str[i+p[i]])
         { if(i+p[i]>p[pos]+pos) Ans=max(Ans,Calc(i-p[i],i+p[i],n)); p[i]++; }
         if(pos+p[pos]<i+p[i])pos=i;
     }

     return ;
 }

 int main()
 {
     int n;

     scanf("%s",str);

     n=strlen(str);
     Rank=Get_SA(n+,);
     Get_H(n);
     Manacher(n);

     printf("%lld",Ans);

     return ;
 }