bzoj 4698 Sandy的卡片 后缀数组

可以把原序列转化为差分序列，题意就变成了求这n个差分序列的最长公共子串。然后就可以连成一串，用后缀数组做。二分一个答案，询问是否有连续n个后缀LCP>=mid而且首位所属的串不相同。

#include<cstdio>
#include<cstring>
#include<iostream>
#define maxn 1100005
using namespace std;
void read(int &a)
{
    a=0;char c=getchar();
    while(c<'0'||c>'9')c=getchar();
    while(c>='0'&&c<='9')
    {
        a*=10;a+=c-'0';
        c=getchar();
    }
}
int id[maxn];
int s[maxn],mx,tot;
struct Suffix
{
    int sa[maxn],rank[maxn],height[maxn];
    int sec[maxn],tmp[maxn],cnt[maxn],n;
    void get_sa()
    {
        n=tot;
        for(int i=1;i<=n;i++) rank[i]=s[i];
        for(int i=1;i<=n;i++) cnt[rank[i]]++;
        for(int i=1;i<=mx;i++)cnt[i]+=cnt[i-1];
        for(int i=1;i<=n;i++) sa[cnt[rank[i]]--]=i;
        for(int k=1;k<=n;k<<=1)
        {
            int top=0;
            for(int i=n-k+1;i<=n;i++) sec[++top]=i;
            for(int i=1;i<=n;i++) if(sa[i]>k) sec[++top]=sa[i]-k;

            for(int i=0;i<=n;i++) cnt[i]=0;
            for(int i=1;i<=n;i++) cnt[rank[i]]++;
            for(int i=1;i<=n;i++) cnt[i]+=cnt[i-1];
            for(int i=n;i>=1;i--) sa[cnt[rank[sec[i]]]--]=sec[i];
            top=tmp[sa[1]]=1;
            for(int i=2;i<=n;i++)
            {
                if(rank[sa[i]]!=rank[sa[i-1]]||rank[sa[i]+k]!=rank[sa[i-1]+k]) top++;
                tmp[sa[i]]=top;
            }
            for(int i=1;i<=n;i++) rank[i]=tmp[i];
        }
        int k=0;
        for(int i=1;i<=n;i++)
        {
            if(k) k--;int j=sa[rank[i]-1];
            while(i+k<=n&&j+k<=n&&s[i+k]==s[j+k])k++;
            height[rank[i]]=k;
        }
    }
    void print()
    {
        for(int i=1;i<=n;i++)
            printf("%d ",s[i]);
        puts("");
        for(int i=1;i<=n;i++)
            printf("%d ",sa[i]);
        puts("");
        for(int i=1;i<=n;i++)
            printf("%d ",height[i]);
    }
    bool use[maxn];
    int S[maxn],top,c;
    bool check(int x,int p)
    {
        memset(use,0,sizeof(use));
        top=c=0;use[0]=1;
        for(int i=1;i<=n;i++)
        {
            if(height[i]<x)
            {while(top) use[S[top--]]=0;c=0;}
            if(!use[id[sa[i]]]&&id[sa[i]])
            {
                use[id[sa[i]]]=1;c++;
                S[++top]=id[sa[i]];
            }
            if(c==p) return true; 
        }
        return false;
    }
}SA;
int n,mn,inf=1e9;
int len[1005];
int a[1005][105];
int main()
{
    scanf("%d",&n);
    int l=0,r=inf;mn=inf;
    for(int i=1;i<=n;i++)
    {
        read(len[i]);r=min(r,len[i]+1);
        for(int j=1;j<=len[i];j++)
        {
            read(a[i][j]);
            if(j!=1)mx=max(mx,a[i][j]-a[i][j-1]);
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=2;j<=len[i];j++)
        {
            s[++tot]=a[i][j]-a[i][j-1];
            mn=min(mn,s[tot]);id[tot]=i;
        }
        s[++tot]=++mx;
        mn=min(mn,s[tot]);
    }
    for(int i=1;i<=tot;i++)
        s[i]=s[i]-mn+1;
    mx=mx-mn+1;
    SA.get_sa();//SA.print();
    while(r-l>1)
    {
        int mid=l+r>>1;
        if(SA.check(mid,n)) l=mid;
        else r=mid;
    }
    printf("%d",l+1);
    return 0;
}