BZOJ 2946 SA/SAM

思路：

1. 二分+后缀数组

2.SAM

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=;
int n,num,cntA[N],cntB[N],A[N],B[N],sa[N],rk[N],tsa[N],ht[N],rec[N],vis[];
char ch[],s[N];
void SA(){
    for(int i=;i<=n;i++)cntA[s[i]]++;
    for(int i=;i<=;i++)cntA[i]+=cntA[i-];
    for(int i=n;i;i--)sa[cntA[s[i]]--]=i;
    rk[sa[]]=;
    for(int i=;i<=n;i++)rk[sa[i]]=rk[sa[i-]]+(s[sa[i]]!=s[sa[i-]]);
    for(int l=;rk[sa[n]]<n;l<<=){
        memset(cntA,,sizeof(cntA));
        memset(cntB,,sizeof(cntB));
        for(int i=;i<=n;i++)
            cntA[A[i]=rk[i]]++,
            cntB[B[i]=i+l<=n?rk[i+l]:]++;
        for(int i=;i<=n;i++)cntA[i]+=cntA[i-],cntB[i]+=cntB[i-];
        for(int i=n;i;i--)tsa[cntB[B[i]]--]=i;
        for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i];
        rk[sa[]]=;
        for(int i=;i<=n;i++)
            rk[sa[i]]=rk[sa[i-]]+(A[sa[i]]!=A[sa[i-]]||B[sa[i]]!=B[sa[i-]]);
    }
    for(int i=,j=;i<=n;i++){
        j=j?j-:;
        while(s[i+j]==s[sa[rk[i]-]+j])j++;
        ht[rk[i]]=j;
    }
}
bool check(int mid){
    int cnt=;
    for(int i=;i<=n;i++){
        if(ht[i]<mid)memset(vis,,sizeof(vis)),cnt=;
        else{
            if(!vis[rec[sa[i]]])vis[rec[sa[i]]]=,cnt++;
            if(!vis[rec[sa[i-]]])vis[rec[sa[i-]]]=,cnt++;
        }
        if(cnt>=num)return ;
    }return ;
}
int main(){
    scanf("%d",&num);
    for(int i=;i<=num;i++){
        scanf("%s",ch+);
        int templen=strlen(ch+);
        for(int j=;j<=templen;j++)s[++n]=ch[j],rec[n]=i;
        s[++n]=i;
    }SA();
    int l=,r=n,ans=;
    while(l<=r){
        int mid=(l+r)>>;
        if(check(mid))ans=mid,l=mid+;
        else r=mid-;
    }printf("%d\n",ans);
}

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=;
int n,m;char s[N];
struct SAM{
    int ch[N][],fa[N],dis[N],root,last,tot,ans[N],v[N],q[N],maxx[N];
    void init(){root=last=tot=;}
    int newnode(int v){return dis[++tot]=v,tot;}
    void add(int x){
        int p=last,np=newnode(dis[p]+);last=np;
        for(;p&&!ch[p][x];p=fa[p])ch[p][x]=np;
        if(!p)fa[np]=root;
        else{
            int q=ch[p][x];
            if(dis[q]==dis[p]+)fa[np]=q;
            else{
                int nq=newnode(dis[p]+);
                memcpy(ch[nq],ch[q],sizeof(ch[q]));
                fa[nq]=fa[q],fa[np]=fa[q]=nq;
                for(;ch[p][x]==q;p=fa[p])ch[p][x]=nq;
            }
        }
    }
    void topsort(){
        for(int i=;i<=tot;i++)ans[i]=dis[i],v[dis[i]]++;
        for(int i=;i<=tot;i++)v[i]+=v[i-];
        for(int i=tot;i;i--)q[v[dis[i]]--]=i;
    }
    void solve(){
        memset(maxx,,sizeof(maxx));
        scanf("%s",s),n=strlen(s);
        int p=root,len=;
        for(int i=;i<n;i++){
            int x=s[i]-'a';
            if(ch[p][x])len++,p=ch[p][x];
            else{
                for(;p&&!ch[p][x];p=fa[p]);
                if(!p)p=root,len=;
                else len=dis[p]+,p=ch[p][x];
            }maxx[p]=max(maxx[p],len);
        }
        for(int i=tot;i;i--){
            int x=q[i];
            ans[x]=min(ans[x],maxx[x]);
            if(maxx[x]&&fa[x])maxx[fa[x]]=max(maxx[fa[x]],dis[fa[x]]);
        }
    }
}T;
int main(){
    T.init();
    scanf("%d%s",&m,s),n=strlen(s);
    for(int i=;i<n;i++)T.add(s[i]-'a');
    T.topsort();
    for(int i=;i<m;i++)T.solve();
    for(int i=;i<=T.tot;i++)T.ans[]=max(T.ans[],T.ans[i]);
    printf("%d\n",T.ans[]);
}