【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=20000+100;
const int INF=2147000000;
int s[maxn],a[maxn],n;
struct state{
    int len,link;
    int next[180];
}st[2*maxn];
int last,cur,sz;
int cnt[2*maxn],l[2*maxn],r[2*maxn],c[2*maxn];
void init(){
    sz=1;
    last=cur=0;
    st[0].len=0;
    st[0].link=-1;
    memset(st[0].next,0,sizeof(st[0].next));
}

void build_sam(int c,int pos){
    cur=sz++;
    st[cur].len=st[last].len+1;
    cnt[cur]=1;
    l[cur]=r[cur]=pos;
    memset(st[cur].next,0,sizeof(st[cur].next));
    int p;
    for(p=last;p!=-1&&st[p].next[c]==0;p=st[p].link)
        st[p].next[c]=cur;
    if(p==-1)
        st[cur].link=0;
    else{
        int q=st[p].next[c];
        if(st[q].len==st[p].len+1)
            st[cur].link=q;
        else{
            int clone=sz++;
            cnt[clone]=r[clone]=0;
            l[clone]=127;
            st[clone].len=st[p].len+1;
            //printf("%d ",st[clone].len);
            st[clone].link=st[q].link;
            memcpy(st[clone].next,st[q].next,sizeof(st[clone].next));
            for(;p!=-1&&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=0;
int main(){
    while(scanf("%d",&n)!=EOF&&n){
        for(int i=1;i<=n;i++){
            scanf("%d",&s[i]);
            a[i]=s[i]-s[i-1];
        }
        init();
        for(int i=1;i<=n;i++){
            build_sam(a[i]+88,i);
            //printf("%d ",a[i]);
        }
//        for(int i=0;i<sz;i++)
//            printf("%d ",st[i].len);
//        printf("
");
        for(int i=0;i<sz;i++)
            c[i]=i;
        sort(c,c+sz,cmp);
        ans=0;
        for(int i=0;i<sz;i++){
            int o=c[i];
            if(st[o].link!=-1){
                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
",cnt[o],st[o].len,l[o],r[o]);
            if(cnt[o]>=2&&min(st[o].len,r[o]-l[o])>=4){
                //printf("%d %d
",l[o],r[o]);
                ans=max(ans,min(st[o].len,r[o]-l[o]));
            }
        }
        if(ans<4)
            printf("0
");
        else
        printf("%d
",ans+1);
    }
return 0;
}