[ 后缀数组 ] [ SDOI2008 ] BZOJ4698 Sandy的卡片

最新推荐文章于 2020-02-27 11:11:05 发布

gjghfd

最新推荐文章于 2020-02-27 11:11:05 发布

阅读量195

点赞数

分类专栏：后缀数组

本文链接：https://blog.csdn.net/gjghfd/article/details/79809360

版权

后缀数组专栏收录该内容

2 篇文章 0 订阅

订阅专栏

先差分，然后就转化为求最长公共子串。
将所有串依次接在一起，之间用分隔符隔开。那么就是求分属 $n$ 个串的后缀的 $lcp$ 最大值。
求出后缀数组，二分答案，把 $height_i\geq mid$ 的区间找出来，判断是否有 n <script type="math/tex" id="MathJax-Element-8">n</script> 个串就可以了。

#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
typedef pair<int,int> abcd;
const int N=111010;
const int INF=1e9;
int k,n,m,ncnt,cnt,cur;
int a[N],c[N];
int t[N],tp[N],rk[N],sa[N],h[N];
int A[N];
int v[N];
abcd d[N];
bool Cmp(int p,int l) {
    return tp[sa[p]]==tp[sa[p-1]]&&tp[sa[p]+l]==tp[sa[p-1]+l];
}
void Rsort(int m) {
    for(int i=1;i<=m;i++) t[i]=0;
    for(int i=1;i<=n;i++) ++t[rk[tp[i]]];
    for(int i=2;i<=m;i++) t[i]+=t[i-1];
    for(int i=n;i;i--) sa[t[rk[tp[i]]]--]=tp[i];
}
void init() {
    for(int i=1;i<=n;i++) rk[i]=a[i],tp[i]=i;
    Rsort(m);
    int p=0;
    for(int j=1;;j<<=1,m=p,p=0) {
        for(int i=n-j+1;i<=n;i++) tp[++p]=i;
        for(int i=1;i<=n;i++)
            if(sa[i]>j) tp[++p]=sa[i]-j;
        Rsort(m);swap(tp,rk);rk[sa[1]]=p=1;
        for(int i=2;i<=n;i++) rk[sa[i]]=Cmp(i,j)?p:++p;
        if(p==n) break;
    }
    p=0;int j;
    for(int i=1;i<=n;h[rk[i++]]=p)
        for(p=p?p-1:0,j=sa[rk[i]-1];a[i+p]==a[j+p];p++);
}
bool Check(int x) {
    int cnt=0;
    for(int i=1;i<=n;i++) {
        if(h[i]<x) ++cur,cnt=0;
        if(c[sa[i]]&&v[c[sa[i]]]!=cur) v[c[sa[i]]]=cur,++cnt;
        if(cnt==ncnt) return 1;
    }
    return 0;
}
int main() {
    scanf("%d",&ncnt);
    for(int i=1;i<=ncnt;i++) {
        scanf("%d",&m);
        for(int j=1;j<=m;j++) scanf("%d",&A[j]);
        for(int j=1;j<=m;j++) d[++n].fi=A[j]-A[j-1],d[n].se=n,c[n]=i;
        d[++n].fi=INF+i;d[n].se=n;
    }
    sort(d+1,d+n+1);
    m=0;
    for(int i=1;i<=n;i++) {
        if(i==1||d[i].fi!=d[i-1].fi) ++m;
        a[d[i].se]=m;
    }
    init();
    int ss=0;
    int l=1,r=n;
    while(l<=r) {
        int Mid=l+r>>1;
        if(Check(Mid)) l=Mid+1;else r=Mid-1;
    }
    cout<<l<<endl;
    return 0;
}