poj_1743_Musical Theme(后缀数组)

题目链接：poj_1743_Musical Theme

题意：

给你一串数字，让你找最长的变化相同不重叠的子串，至少长度为5

题解：

处理数据后用后缀数组加二分答案，然后用height数组check答案，运用height数组求相同不重叠的子串经典运用

 #include<cstdio>
 #include<algorithm>
 #define F(i,a,b) for(int i=a;i<=b;i++)
 using namespace std;

 namespace suffixarray{
     #define FN(n) for(int i=0;i<n;i++)
     const int N =2E4+;
     int rnk[N],sa[N],height[N],c[N],s[N];
     void getsa(int n,int m,int *x=rnk,int *y=height){
         FN(m)c[i]=;FN(n)c[x[i]=s[i]]++;FN(m)c[i+]+=c[i];
         for(int i=n-;i>=;i--)sa[--c[x[i]]]=i;
         for(int k=,p;p=,k<=n;k=p>=n?N:k<<,m=p){
             for(int i=n-k;i<n;i++)y[p++]=i;
             FN(n)if(sa[i]>=k)y[p++]=sa[i]-k;
             FN(m)c[i]=;FN(n)c[x[y[i]]]++;FN(m)c[i+]+=c[i];
             for(int i=n-;i>=;i--)sa[--c[x[y[i]]]]=y[i];
             swap(x,y),p=,x[sa[]]=;
             for(int i=;i<n;i++)
             x[sa[i]]=y[sa[i-]]==y[sa[i]]&&y[sa[i-]+k]==y[sa[i]+k]?p-:p++;
         }
         FN(n)rnk[sa[i]]=i;
         for(int i=,j,k=;i<n-;height[rnk[i++]]=k)
         for(k=k?k-:k,j=sa[rnk[i]-];s[i+k]==s[j+k];k++);
     }
 }

 using namespace suffixarray;
 int n;
 inline void upd(int &a,int b){if(a>b)a=b;}
 inline void upu(int &a,int b){if(a<b)a=b;}

 inline bool check(int x)
 {
     int l=N,r=;
     F(i,,n)
     {
         if(height[i]>=x)
         {
             upd(l,sa[i]),upu(r,sa[i]);
             if(r-l>=x)return ;
         }else l=r=sa[i];
     }
     return ;
 }

 int main()
 {
     while(scanf("%d",&n),n)
     {
         F(i,,n-)scanf("%d",s+i);
         if(n<){puts("");continue;}
         F(i,,n-)s[i]=s[i+]-s[i]+;
         s[--n]=,getsa(n+,);
         int l=,r=n,mid;
         while(l<=r)mid=(l+r)>>,check(mid)?l=mid+:r=mid-;
         printf("%d\n",l<?:l);
     }
     return ;
 }