BZOJ 4698: Sdoi2008 Sandy的卡片（SA+二分）

Description

Sandy和Sue的热衷于收集干脆面中的卡片。然而，Sue收集卡片是因为卡片上漂亮的人物形象，而Sandy则是为了积
攒卡片兑换超炫的人物模型。每一张卡片都由一些数字进行标记，第i张卡片的序列长度为Mi，要想兑换人物模型
，首先必须要集够N张卡片，对于这N张卡片，如果他们都有一个相同的子串长度为k，则可以兑换一个等级为k的人
物模型。相同的定义为：两个子串长度相同且一个串的全部元素加上一个数就会变成另一个串。Sandy的卡片数远
远小于要求的N，于是Sue决定在Sandy的生日将自己的卡片送给Sandy，在Sue的帮助下，Sandy终于集够了N张卡片
，但是，Sandy并不清楚他可以兑换到哪个等级的人物模型，现在，请你帮助Sandy和Sue，看看他们最高能够得到
哪个等级的人物模型。
Input

第一行为一个数N，表示可以兑换人物模型最少需要的卡片数，即Sandy现在有的卡片数
第i+1行到第i+N行每行第一个数为第i张卡片序列的长度Mi，之后j+1到j+1+Mi个数，用空格分隔，分别表示序列中
的第j个数
n<=1000,M<=1000,2<=Mi<=101
Output

一个数k，表示可以获得的最高等级。

---

题解：
这道题只是在求多个串的最长公共子串问题上稍微变化了一点，我们只需要把数列转化成两两的差值就能确定了，因为两个数列通过加法可以相等的话，那么这个数列的差值数列一定是相同的，然后就可以直接二分答案AC啦

---

AC代码：

#pragma GCC optimize(2)
#include<bits/stdc++.h>
#include<ext/rope>
using namespace std;
using namespace __gnu_cxx;
#define LL long long
#define pii pair<int,int>
#define mp(a,b) make_pair(a,b)
const int MAXN = 2e5+1000;
const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
const int base = 2000;
int s[MAXN],a[MAXN],bel[MAXN],ls;
struct SA{
    int rk[MAXN],sa[MAXN],h[MAXN],y[MAXN],c[MAXN];
    int vis[1005];
    inline void get(int n,int m=5500){
        for(int i=0;i<=m;i++) c[i]=0;
        for(int i=1;i<=n;i++) ++c[rk[i]=s[i]];
        for(int i=1;i<=m;i++) c[i]+=c[i-1];
        for(int i=n;i>=1;i--) sa[c[rk[i]]--]=i;
        for(int k=1;k<=n;k<<=1){
            int num=0;
            for(int i=n-k+1;i<=n;i++) y[++num]=i;
            for(int i=1;i<=n;i++) if(sa[i]>k) y[++num]=sa[i]-k;
            for(int i=1;i<=m;i++) c[i]=0;
            for(int i=1;i<=n;i++) ++c[rk[i]];
            for(int i=1;i<=m;i++) c[i]+=c[i-1];
            for(int i=n;i>=1;i--) sa[c[rk[y[i]]]--]=y[i];
            swap(rk,y);
            rk[sa[1]]=num=1;
            for(int i=2;i<=n;i++)
                rk[sa[i]]=(y[sa[i]]==y[sa[i-1]] && y[sa[i]+k]==y[sa[i-1]+k] ? num:++num);
            if(num==n) break;
            m = num;
        }
        for(int i=1,k=0;i<=n;i++){
            if(k) --k; else k=0;
            int j=sa[rk[i]-1];
            while(s[i+k]==s[j+k]) ++k;
            h[rk[i]]=k;
        }
    }
    inline bool check(int mid,int n){
        int cnt=0;
        for(int i=1;i<=ls;i++){
            if(h[i]>=mid){
                if(bel[sa[i-1]] && !vis[bel[sa[i-1]]]) vis[bel[sa[i-1]]]=1,cnt++;
                if(bel[sa[i]] && !vis[bel[sa[i]]]) vis[bel[sa[i]]]=1,cnt++;
                if(cnt==n) return true;
            }else{
                cnt=0; memset(vis,0,sizeof(vis));
            }
        }
        return false;
    }
}sa;
int main(){
#ifndef ONLINE_JUDGE
    freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
#endif // ONLINE_JUDGE
    int n; scanf("%d",&n);
    for(int i=1;i<=n;i++){
        int m; scanf("%d",&m);
        for(int k=1;k<=m;k++) scanf("%d",&a[k]);
        for(int k=2;k<=m;k++) s[++ls]=a[k]-a[k-1]+base,bel[ls]=i;
        s[++ls]=2*base+i;//可能为负数，加一个base
    }
    sa.get(ls);
    int l=0,r=10000,ans=0;
    while(l<=r){
        int mid=(l+r)>>1;
        if(sa.check(mid,n)) l=mid+1,ans=mid;
        else r=mid-1;
    }
    printf("%d\n",ans+1);
    return 0;
}