省选专练之后缀自动机SPOJ1811LCS - Longest Common Substring

陈老师引入SAM的例题

求两个串的LCS

直接建一个串的SAM，并且把另外一个串放进去跳。

这个方法很经典！

几乎满足了所有两个串的公共串问题的所有解法

有巨子说这个跳是均摊Log的，我也不知道QwQ对不对，反正SPOJ那么严苛的时间是过得了的

#include<bits/stdc++.h>
using namespace std;
const int N=4e5+1000;
char S[N];
struct SAM{
    int LCS;
    int ans;
    int now;
    struct Node{
        int pre;
        int step;
        int vis[26];
    }SA[N];
    int last,cnt;
    inline void Insert(char C){
        int p=last;
        int np=++cnt;
        last=np;
        SA[np].step=SA[p].step+1;
        for(;!SA[p].vis[C-'a'];p=SA[p].pre)SA[p].vis[C-'a']=np;
        if(!p)SA[np].pre=1;
        else{
            int q=SA[p].vis[C-'a'];
            if(SA[q].step==SA[p].step+1){
                SA[np].pre=q;
            }
            else{
                int nq=++cnt;
                SA[nq].step=SA[p].step+1;
                memcpy(SA[nq].vis,SA[q].vis,sizeof(SA[q].vis));
                SA[nq].pre=SA[q].pre;
                SA[np].pre=SA[q].pre=nq;
                for(;SA[p].vis[C-'a']==q;p=SA[p].pre)SA[p].vis[C-'a']=nq;
            }
        }
    }
    void GetAns(char C){
        if(SA[now].vis[C-'a']){

            LCS++;
            now=SA[now].vis[C-'a'];
        }
        else{
            for(;now&&!SA[now].vis[C-'a'];now=SA[now].pre);
//			cout<<now<<" "<<SA[now].step<<'\n';
            if(!now)LCS=0,now=1;
            else LCS=SA[now].step+1,now=SA[now].vis[C-'a'];
//			cout<<C<<" "<<now<<" "<<LCS<<'\n';
        }
        ans=max(ans,LCS);//			
//		cout<<now<<" "<<C<<" "<<LCS<<'\n';
    }
    void Solve(){
        LCS=0;
        now=1;
        ans=0;
        cnt=last=1;
        int len;
        scanf("%s",S+1);
        len=strlen(S+1);
        for(int i=1;i<=len;i++)Insert(S[i]);
        scanf("%s",S+1);
        len=strlen(S+1);
        for(int i=1;i<=len;i++)GetAns(S[i]);
        cout<<ans;
    }
}Solution;
int main(){
//	freopen("spoj1811.in","r",stdin);
    Solution.Solve();
    return 0;
}