poj1743 Musical Theme（后缀数组+二分答案）

都做个差值，就解决了转调问题，然后就是求最长不重叠重复子串。二分答案即可。

判断不重叠：按h分组后，记录一组中的最小id--mn和最大id--mx，if(mx-mn>=x)即存在不重叠的。

#include <cstdio>
#include <cstring>
#define N 20005
inline int min(int x,int y){return x<y?x:y;}
inline int max(int x,int y){return x>y?x:y;}
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
	return x*f;
}
int n,a[N],aa[N],h[N];
int rank[N<<1],rank1[N],sa[N],count[N],st[200],tmp[N];
inline bool jud(int x){
	//连续的h[i]大于x，应该都满足，所以分在一组看有没有不重叠的 即可 
	//注意好h[i]是和前一个比，所以要一直更新mx，mn 
	int mn=sa[1],mx=sa[1];
	for(int i=2;i<=n;++i){
		if(h[i]>=x){
			mn=min(mn,sa[i]);mx=max(mx,sa[i]);
			if(mx-mn>=x) return true;
		}
		else{ mn=mx=sa[i];}
	}return false;
}
int main(){
//	freopen("a.in","r",stdin);
	while(1){
	scanf("%d",&n);if(n==0) break;
	for(int i=1;i<=n;++i) aa[i]=read();--n;
	for(int i=1;i<=n;++i) a[i]=aa[i+1]-aa[i]+88;
	memset(st,0,sizeof(st));
	memset(rank,0,sizeof(rank));
	for(int i=1;i<=n;++i) st[a[i]]=1;
	for(int i=1;i<=199;++i) st[i]+=st[i-1];
	for(int i=n;i>=1;--i) rank[i]=st[a[i]];
	int k=0;
	for(int p=1;k!=n;p<<=1){
		memset(count,0,sizeof(count));
		for(int i=1;i<=n;++i) count[rank[i+p]]++;
		for(int i=1;i<=n;++i) count[i]+=count[i-1];
		for(int i=n;i>=1;--i) tmp[count[rank[i+p]]--]=i;
		memset(count,0,sizeof(count));
		for(int i=1;i<=n;++i) count[rank[tmp[i]]]++;
		for(int i=1;i<=n;++i) count[i]+=count[i-1];
		for(int i=n;i>=1;--i) sa[count[rank[tmp[i]]]--]=tmp[i];
		memcpy(rank1,rank,sizeof(rank1));
		rank[sa[1]]=k=1;
		for(int i=2;i<=n;++i){
			if(rank1[sa[i]]!=rank1[sa[i-1]]||rank1[sa[i]+p]!=rank1[sa[i-1]+p]) ++k;
			rank[sa[i]]=k;
		}
	}k=0;
	for(int i=1;i<=n;++i){
		if(rank[i]==1){h[1]=0;continue;}
		if(i==1||h[rank[i-1]]<=1) k=0;
		if(k) --k;
		while(a[i+k]==a[sa[rank[i]-1]+k]) ++k;
		h[rank[i]]=k;
	}
	int l=4,r=n>>1,ans=0;
	while(l<=r){
		int mid=(l+r)>>1;
		if(jud(mid)){ans=mid;l=mid+1;}
		else r=mid-1;
	}
	if(ans<4) puts("0");
	else if((ans<<1)==n) printf("%d\n",ans);
	else printf("%d\n",ans+1);
	}
	return 0;
}