【poj1743】Musical Theme 【后缀自动机】

题意

给出一个n个数字的序列，找出相同变化趋势且不重叠的两个最长子串。

分析

这个题以前应该用后缀数组+二分做过。学了后缀自动机后可以用后缀自动机搞一下。

先差分，然后把查分后的数组建SAM。然后对于每个状态记录一个l[u],和r[u]，分别代表right集合中，最大的v和最小的v。（这里如果不明白可以去看clj的课件）。

然后对于每个状态，当这个状态cnt[u]>=2的时候，说明有两个以上的子串，然后min(st[u].len,r[u]-l[u])就是这个状态最长不重叠相同子串的长度。

 #include <cstring>
 #include <algorithm>
 #include <iostream>
 #include <cstdio>
 #include <map>
 using namespace std;
 const int maxn=+;
 const int INF=;
 int s[maxn],a[maxn],n;
 struct state{
     int len,link;
     int next[];
 }st[*maxn];
 int last,cur,sz;
 int cnt[*maxn],l[*maxn],r[*maxn],c[*maxn];
 void init(){
     sz=;
     last=cur=;
     st[].len=;
     st[].link=-;
     memset(st[].next,,sizeof(st[].next));
 }

 void build_sam(int c,int pos){
     cur=sz++;
     st[cur].len=st[last].len+;
     cnt[cur]=;
     l[cur]=r[cur]=pos;
     memset(st[cur].next,,sizeof(st[cur].next));
     int p;
     for(p=last;p!=-&&st[p].next[c]==;p=st[p].link)
         st[p].next[c]=cur;
     if(p==-)
         st[cur].link=;
     else{
         int q=st[p].next[c];
         if(st[q].len==st[p].len+)
             st[cur].link=q;
         else{
             int clone=sz++;
             cnt[clone]=r[clone]=;
             l[clone]=;
             st[clone].len=st[p].len+;
             //printf("%d ",st[clone].len);
             st[clone].link=st[q].link;
             memcpy(st[clone].next,st[q].next,sizeof(st[clone].next));
             for(;p!=-&&st[p].next[c]==q;p=st[p].link){
                 st[p].next[c]=clone;
             }
             st[cur].link=st[q].link=clone;
         }
     }
     last=cur;
 }
 int cmp(int a,int b){
     return st[a].len>st[b].len;
 }

 int ans=;
 int main(){
     while(scanf("%d",&n)!=EOF&&n){
         for(int i=;i<=n;i++){
             scanf("%d",&s[i]);
             a[i]=s[i]-s[i-];
         }
         init();
         for(int i=;i<=n;i++){
             build_sam(a[i]+,i);
             //printf("%d ",a[i]);
         }
 //        for(int i=0;i<sz;i++)
 //            printf("%d ",st[i].len);
 //        printf("\n");
         for(int i=;i<sz;i++)
             c[i]=i;
         sort(c,c+sz,cmp);
         ans=;
         for(int i=;i<sz;i++){
             int o=c[i];
             if(st[o].link!=-){
                 cnt[st[o].link]+=cnt[o];
                 l[st[o].link]=min(l[st[o].link],l[o]);
                 r[st[o].link]=max(r[st[o].link],r[o]);
             }
          //   printf("%d %d %d %d\n",cnt[o],st[o].len,l[o],r[o]);
             if(cnt[o]>=&&min(st[o].len,r[o]-l[o])>=){
                 //printf("%d %d\n",l[o],r[o]);
                 ans=max(ans,min(st[o].len,r[o]-l[o]));
             }
         }
         if(ans<)
             printf("0\n");
         else
         printf("%d\n",ans+);
     }
 return ;
 }