POJ1743 Musical Theme [后缀自动机]

时间:2023-03-09 13:33:16
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+;
inline int read(){
char c=getchar();int x=,f=;
while(c<''||c>''){if(c=='-')f=-; c=getchar();}
while(c>=''&&c<=''){x=x*+c-''; 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[],par,val;
}t[N];
int sz,root,last;
inline int nw(int _){t[++sz].val=_;return sz;}
inline void iniSAM(){sz=;root=last=nw();}
void extend(int c,int i){
int p=last,np=nw(t[p].val+); d[np]=;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+) t[np].par=q;
else{
int nq=nw(t[p].val+);
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=;i<=n;i++) c[i]=;
for(int i=;i<=sz;i++) c[t[i].val]++;
for(int i=;i<=n;i++) c[i]+=c[i-];
for(int i=sz;i>=;i--) a[c[t[i].val]--]=i;
}
void solve(){
RadixSort();
int ans=;
for(int i=sz;i>=;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]>=&&min(t[u].val,r[u]-l[u]-)>=) ans=max(ans,min(t[u].val,r[u]-l[u]-));
if(t[u].val<) break;
}
if(ans<) puts("");
else printf("%d\n",ans+);
}
int main(){
freopen("in","r",stdin);
while((n=read())){
n--;
memset(t,,sizeof(t));
memset(d,,sizeof(d));memset(r,,sizeof(r));memset(l,,sizeof(l));
iniSAM();
int last=read(),x;
for(int i=;i<=n;i++){
x=read();
s[i]=x-last+;
last=x;
extend(s[i],i);
}
solve();
}
}