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);
}

 

posted @ 2017-04-01 08:18  SiriusRen  阅读(384)  评论(0编辑  收藏  举报