【codevs3160】 LCS 【后缀自动机】

题意

  给出两个字符串，求它们的最长公共子串。

分析

  后缀自动机的基础应用。

  比如两个字符串s1和s2，我们把s1建为SAM，然后根据s2跑，找出s2每个前缀的最长公共后缀。

  我们可以理解为，当向尾部增加一个字母的时候，就按照匹配路径来走，当在SAM中找不到这样的字符串的时候，就要减少头部的字母，就顺着parent树往祖先结点走。

  

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <map>

using namespace std;
const int maxn=1e5+100;
char S[maxn],T[maxn];
int n,m;
struct state{
    int len,link;
    map<char,int>next;
}st[4*maxn];
int last,cur,sz;
void init(){
    last=cur=0;
    sz=1;
    st[0].len=0;
    st[0].link=-1;
//    st[0].next.clear();
}
void build_sam(char c){
    cur=sz++;
    st[cur].len=st[last].len+1;
    int p;
    for(p=last;p!=-1&&!st[p].next.count(c);p=st[p].link){
        st[p].next[c]=cur;
    }
    if(p==-1)
        st[cur].link=0;
    else{
        int q=st[p].next[c];
        if(st[p].len+1==st[q].len){
            st[cur].link=q;
        }else{
            int clone=sz++;
            st[clone].len=st[p].len+1;
            st[clone].next=st[q].next;
            st[clone].link=st[q].link;
            for(;p!=-1&&st[p].next[c]==q;p=st[p].link){
                st[p].next[c]=clone;
            }
            st[cur].link=st[q].link=clone;
        }
    }
    last=cur;
}

int main(){
    scanf("%s%s",S,T);
    n=strlen(S);
    m=strlen(T);
    init();
    for(int i=0;i<n;i++){
        build_sam(S[i]);
    }
    int cur=0;
    int v=0,len=0,ans=0;
    for(int i=0;i<m;i++){
        if(st[v].next.count(T[i])){
            v=st[v].next[T[i]];
            len++;
        }else{
            while(v!=-1&&!st[v].next.count(T[i])){
               // printf("%d
",v);
                v=st[v].link;
                len=st[v].len;
            }
            if(v==-1)
                v=len=0;
            else{
                v=st[v].next[T[i]];
                len++;
            }
        }
        ans=max(ans,len);
    }
    printf("%d
",ans);


return 0;
}