BZOJ 4698 差分+后缀数组
思路:
对所有序列差分一下
公共串的长度+1就是答案了
二分 扫一遍height即可,..
//By SiriusRen #include <cstdio> #include <cstring> #include <algorithm> using namespace std; const int N=1333; int n,xx,ans,a[N],bb=257,top,s[N*N],pos[N*N],A[N*N],B[N*N],cntA[N*N],cntB[N*N],sa[N*N],rk[N*N],tsa[N*N],ht[N*N],vis[N]; void SA(){ for(int i=1;i<=top;i++)cntA[s[i]]++; for(int i=1;i<N;i++)cntA[i]+=cntA[i-1]; for(int i=top;i;i--)sa[cntA[s[i]]--]=i; rk[sa[1]]=1; for(int i=2;i<=top;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]); for(int l=1;rk[sa[top]]<top;l<<=1){ memset(cntA,0,sizeof(cntA)); memset(cntB,0,sizeof(cntB)); for(int i=1;i<=top;i++) cntA[A[i]=rk[i]]++, cntB[B[i]=i+l<=top?rk[i+l]:0]++; for(int i=1;i<=top;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1]; for(int i=top;i;i--)tsa[cntB[B[i]]--]=i; for(int i=top;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i]; rk[sa[1]]=1; for(int i=1;i<=top;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]); } for(int i=1,j=0;i<=top;i++){ j=j?j-1:0; while(s[i+j]==s[sa[rk[i]-1]+j])j++; ht[rk[i]]=j; } } bool check(int mid){ int num=0; for(int i=1;i<=top;i++){ if(ht[i]<mid)memset(vis,0,sizeof(vis)),num=0; else{ if(!vis[pos[sa[i]]])num++,vis[pos[sa[i]]]=1; if(!vis[pos[sa[i-1]]])num++,vis[pos[sa[i-1]]]=1; if(num==n)return 1; } }return 0; } int main(){ scanf("%d",&n); for(int i=1;i<=n;i++){ scanf("%d",&xx); for(int j=1;j<=xx;j++){ scanf("%d",&a[j]); if(j!=1)s[++top]=a[j]-a[j-1]; pos[top]=i; } s[++top]=++bb; } SA(); int l=0,r=1000; while(l<=r){ int mid=(l+r)>>1; if(check(mid))l=mid+1,ans=mid; else r=mid-1; }printf("%d\n",ans+1); }