bzoj4698 [Sdoi2008] Sandy的卡片(后缀数组+二分答案)

差分，就像poj1743那样，多个串的最长公共子串，就像poj3450那样。然后就过了。据说kmp也可以写过这道题。。。各位可以试试。

#include <cstdio>
#include <cstring>
#define N 1100000
int n=0,m=4000,tot,belong[N],s[N],ans=0;
int sa[N],rank[N<<1],rank1[N],tmp[N],count[N],h[N];
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;
}
inline bool jud(int x){
	int cnt=1;bool vis[1100];
	for(int i=1;i<=n+1;++i){
		if(h[i]>=x){
			if(!vis[belong[sa[i]]]) ++cnt,vis[belong[sa[i]]]=true;
		}
		else{
			if(cnt==tot) return true;
			memset(vis,0,sizeof(vis));
			cnt=1;vis[belong[sa[i]]]=true;
		}
	}
	return false;
}
int main(){
//	freopen("a.in","r",stdin);
	tot=read();
	for(int k=1;k<=tot;++k){
		int len=read();int a[1100];
		for(int i=1;i<=len;++i) a[i]=read();
		for(int i=1;i<len;++i) s[++n]=a[i+1]-a[i]+2000,belong[n]=k;
		s[++n]=++m;belong[n]=k;
	}
	for(int i=1;i<=m;++i) count[i]=0;
	for(int i=1;i<=n;++i) count[s[i]]=1;
	for(int i=1;i<=m;++i) count[i]+=count[i-1];
	for(int i=1;i<=n;++i) rank[i]=count[s[i]];
	int k=0;
	for(int p=1;k!=n;p<<=1,m=k){
		for(int i=1;i<=m;++i) count[i]=0;
		for(int i=1;i<=n;++i) count[rank[i+p]]++;
		for(int i=1;i<=m;++i) count[i]+=count[i-1];
		for(int i=n;i>=1;--i) tmp[count[rank[i+p]]--]=i;
		for(int i=1;i<=m;++i) count[i]=0;
		for(int i=1;i<=n;++i) count[rank[tmp[i]]]++;
		for(int i=1;i<=m;++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(s[i+k]==s[sa[rank[i]-1]+k]) ++k;
		h[rank[i]]=k;
	}
	int l=1,r=1000;
	while(l<=r){
		int mid=(l+r)>>1;
		if(jud(mid)) l=mid+1,ans=mid;
		else r=mid-1;
	}
	if(l==1) puts("1");
	else printf("%d\n",ans+1);
	return 0;
}