bzoj4698(洛谷P2463): Sdoi2008 Sandy的卡片（后缀数组）

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/sunshiness_s/article/details/81879841

题意

有 $n$ 个串，求最大的 $k$ ,使得这 $n$ 个串都可以找出长度为 $k$ 的子串且“相同”（这里“相同”意为:经过共同加、减一个数得到一样）

题解

首先这一定是一道找匹配的题目，但是之间匹配肯定得不到我们想要的，因此把每一位变成与前一位相减得到的，既然串经过共同加、减能变成一样，所以每一位的差肯定是一样的
然后就求rank、sa、height
最后求一步二分找下答案就可以了

代码

#include <cstdio>
#include <cstring>
#define N 1000010
int n,num,a[1010][110],len[N],b[N],c[N],cnt[N],tmp[N],height[N],rank[N<<1],rank1[N<<1],sa[N];
inline void getsa(){
    for(int i=1;i<=num;i++) cnt[b[i]]=1;
    for(int i=1;i<=num;i++) cnt[i]+=cnt[i-1];
    for(int i=1;i<=num;i++) rank[i]=cnt[b[i]]; 
    int k=1;for(int len=1;k!=num;len<<=1){
        memset(cnt,0,sizeof(cnt));
        for(int i=1;i<=num;i++) cnt[rank[i+len]]++;
        for(int i=1;i<=num;i++) cnt[i]+=cnt[i-1];
        for(int i=num;i>=1;i--) tmp[cnt[rank[i+len]]--]=i;

        memset(cnt,0,sizeof(cnt));
        for(int i=1;i<=num;i++) cnt[rank[tmp[i]]]++;
        for(int i=1;i<=num;i++) cnt[i]+=cnt[i-1];
        for(int i=num;i>=1;i--) sa[cnt[rank[tmp[i]]]--]=tmp[i];

        memcpy(rank1,rank,sizeof(rank));
        k=rank[sa[1]]=1;
        for(int i=2;i<=num;i++){
            if(rank1[sa[i]]!=rank1[sa[i-1]] || rank1[sa[i]+len]!=rank1[sa[i-1]+len]) k++;
            rank[sa[i]]=k;
        }
    }
}
inline void getheight(){
    int k=0;
    for(int i=1;i<=num;i++){
        if(rank[i]==1) k=0;
        else{
            if(--k<0) k=0;
            while(b[i+k]==b[sa[rank[i]-1]+k]) k++;
        }
        height[rank[i]]=k;
    } 
}
inline bool check1(int l,int r){
    if(r-l+1<n) return false;
    memset(cnt,0,sizeof(cnt));
    int ans=0;
    for(int i=l;i<=r;i++){
        if(!cnt[c[sa[i]]]){
            cnt[c[sa[i]]]=1;
            ans++;
        }
    }
    if(ans==n) return true;
    return false;
} 
inline bool check(int x){
    int l=1,r=1;
    while(l<=num){
        while(height[r+1]>=x) r++;
        if(check1(l,r)) return true;
        l=r+1;r=l;
    }
    return false;
}
int main(){
    freopen("a.in","r",stdin);
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d%d",&len[i],&a[i][1]);
        for(int j=2;j<=len[i];j++){
            scanf("%d",&a[i][j]);
            b[++num]=a[i][j]-a[i][j-1];
            c[num]=i;
        }
        b[++num]=0;c[num]=i;
    }
    getsa();getheight();
//  for(int i=1;i<=num;i++) printf("%d\n",height[i]);
    int l=1,r=num,ans=0;
    while(l<=r){
        int mid=l+r>>1; 
        if(check(mid)) ans=mid,l=mid+1;
        else r=mid-1;
    }
    printf("%d\n",ans+1);
    return 0;
}