后缀数组+二分--bzoj4698: Sdoi2008 Sandy的卡片

传送门

差分一下，去掉首项，然后二分一个长度看连续的一段大于这个长度的$height$是否能包含所有$n$个序列

代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define maxm 110005
#define maxn 1005
using namespace std;

inline int rd(){
	int x=0,f=1;char c=getchar();
	while(c<'0' || c>'9') f=c=='-'?-1:1,c=getchar();
	while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
	return x*f;
}

int n,a[maxn],s[maxm],m,len,id[maxm];
int tax[maxm],h[maxm],rk[maxm],sa[maxm],tp[maxm];
bool vis[maxn];

inline void rsort(){
	for(int i=1;i<=m;i++) tax[i]=0;
	for(int i=1;i<=len;i++) tax[rk[i]]++;
	for(int i=1;i<=m;i++) tax[i]+=tax[i-1];
	for(int i=len;i;i--) sa[tax[rk[tp[i]]]--]=tp[i];
}

inline void ssort(){
	for(int i=1;i<=len;i++) rk[i]=s[i],tp[i]=i;
	rsort();
	for(int w=1,p=0;w<=len&&p<len;w<<=1,m=p){
		p=0;
		for(int i=len-w+1;i<=len;i++) tp[++p]=i;
		for(int i=1;i<=len;i++)
			if(sa[i]>w) tp[++p]=sa[i]-w;
		rsort(); swap(rk,tp);
		rk[sa[1]]=p=1;
		for(int i=2;i<=len;i++)
			if(tp[sa[i]]==tp[sa[i-1]] && tp[sa[i]+w]==tp[sa[i-1]+w])
				rk[sa[i]]=p;
			else rk[sa[i]]=++p;
	}
}

inline void get_h(){
	int j,k=0;
	for(int i=1;i<=len;i++){
		if(k) --k;
		j=sa[rk[i]-1];
		while(s[j+k]==s[i+k]) ++k;
		h[rk[i]]=k;
	}
}

inline bool check(int x){
	int cnt=0; memset(vis,0,sizeof vis);
	for(int i=1;i<=len;i++){
		if(h[i]>=x){
			if(!vis[id[sa[i]]]) vis[id[sa[i]]]=1,cnt++;
		}
		else{
			memset(vis,0,sizeof vis); cnt=1; vis[id[sa[i]]]=1;
		}
		if(cnt>=n) return true;
	}
	return cnt>=n;
}

int main(){
	n=rd();
	for(int i=1;i<=n;i++){
		int mi=rd();
		for(int j=1;j<=mi;j++) a[j]=rd();
		for(int j=mi;j;j--) a[j]-=a[j-1];
		for(int j=2;j<=mi;j++) s[++len]=a[j],id[len]=i;
		s[++len]=i+2000;
	}
	for(int i=1;i<=len;i++) m=max(m,s[i]);
	ssort(); get_h();
	int l=0,r=101,mid,ans=0;
	while(l<=r){
		mid=(l+r)>>1;
		if(check(mid)) ans=mid,l=mid+1;
		else r=mid-1;
	}
	printf("%d\n",ans+1);
	return 0;
}