SPOJ 1811 SAM 初探

思路：

一个串建SAM

另一个串在SAM上跑

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=;
int n;char a[N],b[N];
struct SAM{
    int ch[N][],dis[N],fa[N],root,tot,last;
    void init(){root=tot=last=;}
    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;
            }
        }
    }
}T;
int main(){
    T.init();
    scanf("%s%s",a+,b+);
    int n=strlen(a+),m=strlen(b+),ans=,len=;
    for(int i=;i<=n;i++)T.add(a[i]-'a');
    int p=T.root;
    for(int i=;i<=m;i++){
        int x=b[i]-'a';
        if(T.ch[p][x])len++,p=T.ch[p][x];
        else{
            for(;p&&!T.ch[p][x];p=T.fa[p]);
            if(!p)p=T.root,len=;
            else len=T.dis[p]+,p=T.ch[p][x];
        }ans=max(ans,len);
    }printf("%d\n",ans);
}