poj 1743 - 后缀自动机或后缀数组

后缀自动机:利用right集合的包含关系，利用一开始的原串的前缀串的最右端，去更新其他点right集合右端点包含它的子串出现的右端点位置，然后用最大位置减去最小位置和l【i】的关系看是否有重复

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int mx = 4e4 + 10;
int n,m,a[mx],tot,last,ma[mx],mi[mx];
int fa[mx],id[mx][190],len[mx];
int v[mx],pos[mx];
void init()
{
	last = tot = 1;
	fa[1] = len[1] = 0;
	memset(id[1],0,sizeof(id[1]));
	memset(v,0,sizeof(v));
}
void suf_auto(int ch)
{
	int np = ++tot,p = last;
	memset(id[np],0,sizeof(id[np]));
	len[np] = len[last] + 1,last = np;
	while(p&&!id[p][ch]) id[p][ch] = np,p = fa[p];
	if(!p) fa[np] = 1;
	else{
		int q = id[p][ch];
		if(len[p]+1==len[q]) fa[np] = q;
		else{
			int nq = ++tot;
			copy(id[q],id[q]+sizeof(id[q]),id[nq]);
			fa[nq] = fa[q];
			len[nq] = len[p] + 1;
			fa[q] = fa[np] = nq;
			while(p&&id[p][ch]==q) id[p][ch] = nq,p = fa[p];
		}
	} 
}
int solve()
{
	int p = 1,ans = 0;
	memset(mi,inf,sizeof(mi));
	memset(ma,0,sizeof(ma));
	for(int i=1;i<n;i++){//前缀串的最右位置 
		 p = id[p][a[i]];
		 ma[p] = mi[p] = i;
	}
	for(int i=tot;i>=1;i--){
		int k = pos[i];//fa包含k的right集合 
		ma[fa[k]] = max(ma[fa[k]],ma[k]);
		mi[fa[k]] = min(mi[fa[k]],mi[k]); 
	}
	for(int i=1;i<=tot;i++){
		if(ma[i]-mi[i]>len[i])
		ans = max(ans,len[i]);
	}
	if(ans<4) return 0;
	return ans+1;
}
int main()
{
	while(scanf("%d",&n)&&n)
	{
		init();
		for(int i=1;i<=n;i++) scanf("%d",a+i);
		for(int i=1;i<n;i++) a[i] = a[i+1] - a[i] + 88;
		for(int i=1;i<n;i++) suf_auto(a[i]);
		for(int i=1;i<=tot;i++) v[len[i]]++;
		for(int i=1;i<=n;i++) v[i] += v[i-1];
		for(int i=1;i<=tot;i++) pos[v[len[i]]--] = i;
		printf("%d\n",solve());
	}
	return 0;
}

后缀数组:算完排名之后，二分非重叠最长长度即可。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<set>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int mx = 2e4 + 10;
int n,m,a[mx];
int c[mx],wa[mx],wb[mx],sa[mx],he[mx],ra[mx];
void sort_sa(int m,int n)
{
	//如果末尾不加一个0作为排名为0的点，那么就有一个非空串就被当成0，跟空串排名一样，就会发生错误 
	int *x = wa,*y = wb, p  = 0;
	for(int i=0;i<m;i++) c[i] = 0;
	for(int i=0;i<n;i++) c[x[i]=a[i]]++;
	for(int i=1;i<m;i++) c[i] += c[i-1];
	for(int i=n-1;i>=0;i--) sa[--c[x[i]]] = i;
	for(int k=1;p<n;k<<=1,m=p){
		p = 0;
		for(int i=n-k;i<n;i++) y[p++] = i;//个位排位p的位置 
		for(int i=0;i<n;i++) if(sa[i]>=k) y[p++] = sa[i]-k;//考虑后一半的排名，那么位置肯定是这个排名的前k个位置 
		for(int i=0;i<m;i++) c[i] = 0;
		for(int i=0;i<n;i++) c[x[y[i]]]++;
		for(int i=1;i<m;i++) c[i] += c[i-1];
		for(int i=n-1;i>=0;i--) sa[--c[x[y[i]]]] = y[i];//个位数大的优先原则 
		swap(x,y),x[sa[0]] = 0,p = 1;//把y看做十位数，因为x已经比过了，没有必要保留了 
		for(int i=1;i<n;i++)
		if(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]) x[sa[i]] = p-1;//两半排名都一样，排名就一样 
		else x[sa[i]] = p++; 
	}
}
void get_he(int n)
{
	int k = 0;
	for(int i=0;i<n;i++) ra[sa[i]] = i;
	for(int i=0;i<n-1;i++){
		if(k) k--;//如果i-1向后延伸有k个那么i肯定至少是k-1个 
		int j = sa[ra[i]-1];
		while(a[i+k]==a[j+k]) k++;
		he[ra[i]] = k; 
	}
}
int main()
{
	while(scanf("%d",&n)&&n)
	{
		for(int i=0;i<n;i++) scanf("%d",a+i);
		for(int i=0;i<n-1;i++) a[i] = a[i+1] - a[i] + 88;
		a[n-1] = 0;
		sort_sa(190,n);
		get_he(n); 
		int l = 0,r = n;
		while(l<r){
			int mid = (l+r+1)>>1,flag = 0;
			int ma = sa[1],mi = sa[1];
			for(int i=2;i<n;i++){
				if(he[i]>=mid){
					mi = min(sa[i],mi);
					ma = max(sa[i],ma);
				}else{
					if(ma-mi>mid) flag = 1;
					ma = sa[i],mi = sa[i];
				}
			}
			if(ma-mi>mid||flag) l = mid;
			else r = mid - 1;
		}
		if(l<4) puts("0");
		else printf("%d\n",l+1);
	}
	return 0;
}