POJ1743 Musical Theme [后缀自动机]

题意:不重叠最长重复子串


后缀数组做法:http://www.cnblogs.com/candy99/p/6227659.html

后缀自动机的话,首先|Right|>=2 然后min(t[u].val,最右r-最左r-1)>=4就可以了

注意判断条件不能写错哦 并且需要两者取min

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=4e4+5;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1; c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0'; c=getchar();}
    return x*f;
}
int n,m,k;
int c[N],a[N],d[N],l[N],r[N];
int s[N];
struct State{
    int ch[180],par,val;
}t[N];
int sz,root,last;
inline int nw(int _){t[++sz].val=_;return sz;}
inline void iniSAM(){sz=0;root=last=nw(0);}
void extend(int c,int i){
    int p=last,np=nw(t[p].val+1); d[np]=1;l[np]=r[np]=i;
    for(;p&&!t[p].ch[c];p=t[p].par) t[p].ch[c]=np;
    if(!p) t[np].par=root;
    else{
        int q=t[p].ch[c];
        if(t[q].val==t[p].val+1) t[np].par=q;
        else{
            int nq=nw(t[p].val+1);
            memcpy(t[nq].ch,t[q].ch,sizeof(t[q].ch));
            t[nq].par=t[q].par;
            t[q].par=t[np].par=nq;
            for(;p&&t[p].ch[c]==q;p=t[p].par) t[p].ch[c]=nq;
        }
    }
    last=np;
}
void RadixSort(){
    for(int i=0;i<=n;i++) c[i]=0; 
    for(int i=1;i<=sz;i++) c[t[i].val]++;
    for(int i=1;i<=n;i++) c[i]+=c[i-1];
    for(int i=sz;i>=1;i--) a[c[t[i].val]--]=i; 
}
void solve(){
    RadixSort();
    int ans=0;
    for(int i=sz;i>=1;i--){
        int u=a[i],par=t[u].par;
        d[par]+=d[u];
        l[par]=min(l[par],l[u]);
        r[par]=max(r[par],r[u]);
        if(d[u]>=2&&min(t[u].val,r[u]-l[u]-1)>=4) ans=max(ans,min(t[u].val,r[u]-l[u]-1));
        if(t[u].val<4) break;
    }
    if(ans<4) puts("0");
    else printf("%d\n",ans+1);
}
int main(){
    freopen("in","r",stdin);
    while((n=read())){
        n--;
        memset(t,0,sizeof(t));
        memset(d,0,sizeof(d));memset(r,0,sizeof(r));memset(l,127,sizeof(l));
        iniSAM();
        int last=read(),x;
        for(int i=1;i<=n;i++){
            x=read();
            s[i]=x-last+88;
            last=x;
            extend(s[i],i);
        }
        solve();
    }
}

 

posted @ 2017-02-08 18:42  Candy?  阅读(215)  评论(0编辑  收藏  举报