[Poj1743] [后缀数组论文例题] Musical Theme [后缀数组不可重叠最长重复子串]

利用后缀数组，先对读入整数处理str[i]=str[i+1]-str[i]+90这样可以避免负数，计算Height数组，二分答案，如果某处H<lim则将H数组分开，最终分成若干块，判断每块中是否存在SA[i]-SA[j]>lim（注意不是>=因为查分后如果相等那么两个公共串连接部分的元素是公用的，不符合题意），注意特判n=1的情况，因为查分后的结果是个空串，会导致RE。

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

 using namespace std;

 int    A[],B[],U[],SA[],H[];
 int    str[];

 void    Calc_H(const int n,const int * Rank)
 {
     int    i,j,k=;
     for(i=;i<n;H[Rank[i++]]=k)
         for(k?k--:,j=SA[Rank[i]-];str[i+k]==str[j+k];++k);
     return ;
 }

 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;
     memset(U,,sizeof(U));
     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;j<<=,m=p)
     {
         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<m;++i)U[i]=;
         for(i=;i<n;++i)U[x[y[i]]]++;
         for(i=;i<m;++i)U[i]+=U[i-];
         for(i=n-;i>=;--i)SA[--U[x[y[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++;
     }
     Calc_H(n,x);
     return x;
 }

 bool    Check(const int lim,const int n)
 {
     int    Min=SA[],Max=SA[];
     for(int i=;i<=n;++i)
     {
         if(i> && H[i]<lim)
         {
             if(Max-Min>=lim+)return true;
             Min=0x7fffffff;
             Max=-;
         }
         if(SA[i]<Min)Min=SA[i];
         if(SA[i]>Max)Max=SA[i];
     }
     if(Max-Min>lim)return true;
     return false;
 }

 void    Init()
 {
     memset(A,,sizeof(A));
     memset(B,,sizeof(B));
     memset(U,,sizeof(U));
     memset(H,,sizeof(H));
     memset(SA,,sizeof(SA));
     return ;
 }

 int main()
 {
     int    n,i,l,r;

     while(~scanf("%d",&n) && n)
     {
         Init();
         for(i=;i<n;++i)scanf("%d",&str[i]);
         if(n==){printf("0\n");continue;}
         n--;
         for(i=;i<n;++i)str[i]=str[i+]-str[i]+;
         str[n]=;

         Get_SA(n+,);

         l=,r=n;
         while(l<r-)
         {
             int    mid=l+((r-l)>>);
             if(Check(mid,n))l=mid;
             else        r=mid;
         }
         printf("%d\n",l>= ? l+:);
     }

     return ;
 }