Code vs 3160 最长公共子串【SAM】

http://codevs.cn/problem/3160/

有篇blog相当不错：https://blog.csdn.net/clover_hxy/article/details/53758535

初学后缀自动机，先按套路，后再丰富；

将一个串录入SAM，然后用另一个串进行匹配。

摘：
ch[x][s]表示的是节点x（也可以说是状态）字符s的节点编号。返回该结点对应的子串加上某个字符后生成的合法子串在后缀自动机中所对应的位置(其实就和字母树一样),如果该指针不存在,就说明这样的子串是不存在的(即不是s的子串)

fa[x]其实就是parent指针，但是注意这里的fa[x]表示的并不是他的父亲是谁，因为上面也说过了一个节点因为代表了多个含义所以他也有多个父亲。所以parent指针指向的位置表示的是如果当前节点可以接某个后缀那么parent指针指向的节点也一定可以接。即上一个可以接受后缀的节点。

right(s)表示的子串s在原串中所有出现位置的右端点的集合。对于right满足right(s)是right(fa(s))的子集。right的求法，按照parent树中深度从大到小，依次将每个状态的right集合并入他fa状态的right集合。

l[x]表示的是从根节点到该节点的最大距离，对于一个状态s,他代表的串的长度区间就是(len[fa],len[x]].

#include <bits/stdc++.h>
using namespace std;
const int N=2e5+2333;
int n,m,ch[N][30],fa[N],a[N],l[N];
int np,nq,p,q,last,cnt,root;
char s[N],s1[N];
void extend(int x){
    int c=a[x];
    p=last;np=++cnt;last=np;
    l[np]=x;
    for(;p&&!ch[p][c];p=fa[p]){
        ch[p][c]=np;
    }
    if(!p) fa[np]=root;
    else{
        int q=ch[p][c];
        if(l[p]+1==l[q]) fa[np]=q;
        else{
            int nq=++cnt;l[nq]=l[p]+1;
            memcpy(ch[nq],ch[q],sizeof ch[nq]);
            fa[nq]=fa[q];
            fa[np]=fa[q]=nq;
            for(;ch[p][c]==q;p=fa[p]) ch[p][c]=nq;
        }
    }
}
int solve()
{
    int tmp=0,ans=0;
    for(int i=1;i<=m;i++){
        int c=s1[i]-'a';
        if(ch[p][c]) p=ch[p][c],tmp++;
        else{
            while(p&&!ch[p][c]) p=fa[p];
            if(!p) p=1,tmp=0;
            else{
                tmp=l[p]+1;
                p=ch[p][c];
            }
        }
        ans=max(ans,tmp);
    }
    return ans;
}
int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);last=root=++cnt;
    for(int i=1;i<=n;i++){a[i]=s[i]-'a';}
    for(int i=1;i<=n;i++) extend(i);
    scanf("%s",s1+1);
    m=strlen(s1+1);
    printf("%d",solve());
    return 0;
}