「luogu2463」[SDOI2008]Sandy的卡片

获得原数列的差分数列，离散化后接在一起，中间用没有出现过的数分隔，于是问题就转化成求这个数列的最长重复子串，且每一段上都至少有一个子串。

求出后缀数组，二分答案即可。

#include<bits/stdc++.h>
#define R register
using namespace std;
const int N=1010,M=1000010,oo=2e9;
int n,m,tot,s[M],a[M],numpos[M];
int x[M<<1],y[M<<1],sa[M<<1],rank[M],h[M],height[M],cnt[M],bel[M];
int cmp(const int xx,const int yy){return a[xx]<a[yy];}
void build(){
    for(R int i=1;i<=n;i++) cnt[x[i]]++;
    for(R int i=2;i<=m;i++) cnt[i]+=cnt[i-1];
    for(R int i=n;i;i--) sa[cnt[x[i]]--]=i;
    for(R int i=1;i<=m;i++) cnt[i]=0;
    for(R int k=1;k<=n;k<<=1){
        int num=0;
        for(R int i=n;i>n-k;i--) y[++num]=i;
        for(R int i=1;i<=n;i++) if(sa[i]>k) y[++num]=sa[i]-k;
        for(R int i=1;i<=n;i++) cnt[x[i]]++;
        for(R int i=2;i<=m;i++) cnt[i]+=cnt[i-1];
        for(R int i=n;i;i--) sa[cnt[x[y[i]]]--]=y[i];
        for(R int i=1;i<=m;i++) cnt[i]=0;
        swap(x,y);
        x[sa[1]]=num=1;
        for(R int i=2;i<=n;i++){
            if(y[sa[i]]!=y[sa[i-1]]||y[sa[i]+k]!=y[sa[i-1]+k]) num++;
            x[sa[i]]=num;
        }
        if(num>=n) break;
        m=num;
    }
    for(R int i=1;i<=n;i++) rank[sa[i]]=i;
    int t=0;
    for(R int i=1;i<=n;i++){
        if(t) t--;
        while(rank[i]-1&&s[i+t]==s[sa[rank[i]-1]+t]) t++;
        h[i]=t,height[rank[i]]=t;
    }
}
bool vis[N];
bool check(int k){
    memset(vis,0,sizeof(vis));
    vis[0]=1;
    int minh=oo,cc=0;
    for(int i=1;i<=n;i++){
        minh=min(minh,height[i]);
        if(minh<k){
            memset(vis,0,sizeof(vis));
            vis[0]=1;
            cc=0,minh=oo;
        }else{
            if(!vis[bel[sa[i]]]) vis[bel[sa[i]]]=1,cc++;
            if(i>=1&&!vis[bel[sa[i-1]]]) vis[bel[sa[i-1]]]=1,cc++;
            if(cc>=tot) return 1;
        }
    }
    return 0;
}
int main(){
    int t1,t2,t3,div=1000000000;
    scanf("%d",&tot);
    for(R int i=1;i<=tot;i++){
        scanf("%d%d",&t1,&t2);
        for(int j=1;j<t1;j++) scanf("%d",&t3),a[++n]=t3-t2,x[n]=a[n],t2=t3,bel[n]=i,numpos[n]=n;
        a[++n]=++div,numpos[n]=n;
    }
    sort(numpos+1,numpos+n+1,cmp);
    int temp=0;
    x[numpos[1]]=s[numpos[1]]=++temp;
    for(int i=2;i<=n;i++){
        if(a[numpos[i]]!=a[numpos[i-1]]) temp++;
        s[numpos[i]]=x[numpos[i]]=temp;
    }
    m=temp;    
    build();
    int l=0,r=n,ans=0;
    while(l<=r){
        int mid=(l+r)>>1;
        if(check(mid)) l=mid+1,ans=mid;
        else r=mid-1;
    }
    printf("%d",ans+1);
    return 0;
}