poj1743后缀数组

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define maxn 50000
#define inf 0x3fffffff
int wa[maxn],wb[maxn],ss[maxn],sa[maxn],wv[maxn];
int cmp(int *r,int a,int b,int c){
	return r[a]==r[b]&&r[a+c]==r[b+c];
}
int n,a[maxn];
void da(int *r,int m){
	int i,j,p,*x=wa,*y=wb;
	for(i=0;i<n;i++)r[i]++;
	r[n++]=0;	
	for(i=0;i<m;i++)ss[i]=0;
	for(i=0;i<n;i++)ss[x[i]=r[i]]++;
	for(i=1;i<m;i++)ss[i]+=ss[i-1];
	for(i=n-1;i>=0;i--)sa[--ss[x[i]]]=i;
	for(j=1,p=1;p<n;j<<=1,m=p){
		for(p=0,i=n-j;i<n;i++)y[p++]=i;
        for(i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
        for(i=0;i<n;i++)wv[i]=x[y[i]];
        for(i=0;i<m;i++)ss[i]=0;
        for(i=0;i<n;i++)ss[wv[i]]++;
        for(i=1;i<m;i++)ss[i]+=ss[i-1];
        for(i=n-1;i>=0;i--)sa[--ss[wv[i]]]=y[i];
        for(swap(x,y),p=1,x[sa[0]]=0,i=1;i<n;i++)
        x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
	}
}
int he[maxn],rank[maxn];
void  make(int *r){
	int i,j,k=0;
	for(i=1;i<n;i++)rank[sa[i]]=i;
	for(i=0;i<n-1;he[rank[i++]]=k)for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
}
char in[1000];
bool check(int k){
	int mi=inf,ma=-1;
	for(int i=1;i<n;i++){	
		if(he[i]<k)mi=inf,ma=-1;
		mi=min(mi,sa[i]);
		ma=max(ma,sa[i]);
		if(ma>mi+k)return 1;
	}
	return 0;
}
int le,ri;
void ask(){
	le=0,ri=n;
	while(le+1!=ri){
		int mid=(le+ri)>>1;
		if(check(mid))le=mid;
		else ri=mid;
	}
  le++;
	if(le>=5)printf("%d\n",le);
	else printf("0\n");
}
int main(){
	while(scanf("%d",&n)&&n){
		int la,no;
		scanf("%d",&la);
		for(int i=1;i<n;i++)scanf("%d",&no),a[i-1]=(no-la+100),la=no;
		da(a,200);
		make(a);
		ask();
	}
	
}

poj1743后缀数组,布布扣,bubuko.com

poj1743后缀数组

原文：http://www.cnblogs.com/wangyucheng/p/3595207.html