POJ 1743 后缀数组不重叠最长重复子串

时间:2022-02-05 13:26:50
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define maxn 30000
using namespace std;
int n;
int jilu[maxn];
int wa[maxn],wb[maxn],wv[maxn],ws[maxn],r[maxn],sa[maxn],rankk[maxn],height[maxn];//r数组是将原始字符串进行整数化
int cmp(int *r,int a,int b,int l){
return r[a]==r[b]&&r[a+l]==r[b+l];
}
void da(int *r,int *sa,int n,int m){//待排序字符串长度是n,从0开始,最大是m-1
int i,j,p,*x=wa,*y=wb,*t;
for(i=;i<m;i++)ws[i]=;
for(i=;i<n;i++)ws[x[i]=r[i]]++;
for(i=;i<m;i++)ws[i]+=ws[i-];
for(i=n-;i>=;i--)sa[--ws[x[i]]]=i;
for(j=,p=;p<n;j*=,m=p){
for(p=,i=n-j;i<n;i++)y[p++]=i;
for(i=;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
for(i=;i<n;i++)wv[i]=x[y[i]];
for(i=;i<m;i++)ws[i]=;
for(i=;i<n;i++)ws[wv[i]]++;
for(i=;i<m;i++)ws[i]+=ws[i-];
for(i=n-;i>=;i--)sa[--ws[wv[i]]]=y[i];
for(t=x,x=y,y=t,p=,x[sa[]]=,i=;i<n;i++)
x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
}
}
void calheight(int *r,int *sa,int n){
int i,j,k=;
for(i=;i<=n;i++) rankk[sa[i]]=i;
for(i=;i<n;height[rankk[i++]]=k)
for(k?k--:,j=sa[rankk[i]-];r[i+k]==r[j+k];k++);
return;
}
bool ok(int mid){
int mmax=-,mmin=;
for(int i=;i<n;i++){
if(height[i]<mid){
if(mmax-mmin>=mid)return ;
mmax=sa[i];
mmin=sa[i];
}
mmax=max(sa[i],mmax);
mmin=min(sa[i],mmin);
}
if(mmax-mmin>=mid)return ;
else return ;
}
int main()
{
while(scanf("%d",&n)!=EOF){
if(n==)break;
for(int i=;i<n;i++){
scanf("%d",&jilu[i]);
if(i>)r[i-]=jilu[i]-jilu[i-]+;
}
r[n-]=;
da(r,sa,n,);
calheight(r,sa,n-);
int l=,r=n;
while(l<=r){
int mid=(l+r)>>;
if(ok(mid))l=mid+;
else r=mid-;
}
r++;
if(r<)r=;
printf("%d\n",r);
}
}