●BZOJ 4698 Sdoi2008 Sandy的卡片

题链:

http://www.lydsy.com/JudgeOnline/problem.php?id=4698

题解:

后缀数组,二分
这个题还是比较套路的。
首先依据题意,把各个串差分以后,用分割符号串在一起,
然后用倍增算法求出 sa[] rank[] height[]
接着就是二分答案,
并对 height[]分组(每一组内的LCP>+X),
判断是否在某一组内出现了所有串的子串 。

代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#define MAXN 1005000
#define INF 0x3f3f3f3f
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
int A[MAXN],S[MAXN];
int sa[MAXN],rak[MAXN],hei[MAXN],bel[MAXN];
void build(int N,int M){
	static int cc[MAXN],ta[MAXN],tb[MAXN],*x,*y,h,p;
	x=ta; y=tb; h=0; A[N]=-1;
	for(int i=0;i<M;i++) cc[i]=0;
	for(int i=0;i<N;i++) cc[x[i]=A[i]]++;
	for(int i=1;i<M;i++) cc[i]+=cc[i-1];
	for(int i=N-1;i>=0;i--) sa[--cc[x[i]]]=i;
	for(int k=1;p=0,k<N;k<<=1){
		for(int i=N-k;i<N;i++) y[p++]=i;
		for(int i=0;i<N;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
		for(int i=0;i<M;i++) cc[i]=0;
		for(int i=0;i<N;i++) cc[x[y[i]]]++;
		for(int i=1;i<M;i++) cc[i]+=cc[i-1];
		for(int i=N-1;i>=0;i--) sa[--cc[x[y[i]]]]=y[i];
		swap(x,y); y[N]=-1; x[sa[0]]=0; M=1;
		for(int i=1;i<N;i++)
			x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?M-1:M++;
		if(M>=N) break;
	}
	for(int i=0;i<N;i++) rak[sa[i]]=i;
	for(int i=0,j;i<N;i++){
		if(h) h--;
		if(rak[i]){
			j=sa[rak[i]-1];
			while(A[i+h]==A[j+h]) h++;
		}
		hei[rak[i]]=h;
	}
}
bool check(int X,int N,int K){
	static int vis[MAXN],tim,now; 
	tim++; now=1;
	for(int i=0;i<N;i++){
		if(hei[i]<X) tim++,vis[bel[sa[i]]]=tim,now=1;
		else{
			if(vis[bel[sa[i]]]==tim) continue;
			vis[bel[sa[i]]]=tim;
			now++;
		}
		if(now==K) return 1;
	}
	return 0;
}
int binary_search(int r,int N,int K){
	int l=1,ans=0,mid;
	while(l<=r){
		mid=(l+r)>>1;
		if(check(mid,N,K)) ans=mid,l=mid+1;
		else r=mid-1;
	}
	return ans;
}
int main()
{
	int N,M,len=0,r=INF;
	scanf("%d",&N);
	for(int i=0;i<N;i++){
		scanf("%d",&M); r=min(r,M-1);
		for(int j=0;j<M;j++){
			scanf("%d",&S[j]);
			if(j) S[j-1]=S[j]-S[j-1],A[len+j-1]=S[j-1],bel[len+j-1]=i;
		}
		A[len+M-1]=2000+i;
		len=len+M;
	}
	build(len,3005);
	int ans=binary_search(r,len,N)+1;
	printf("%d",ans);
	return 0;
}



Do not go gentle into that good night.
Rage, rage against the dying of the light.
————Dylan Thomas
posted @ 2017-12-07 19:34  *ZJ  阅读(111)  评论(0编辑  收藏  举报