SPOJ 1811 LCS 后缀自动机

题意：

给两个字符串A、B求他们的最长公共字串。

 

题解：

后缀自动机啊。你怎么这么恶心这么神啊。。。

必然还是我太弱了。。

第一次搞，参考了别人的代码。

将A建立成后缀自动机，后缀自动机的任意一个节点都表示若干个A的字串，让B在后缀自动机上匹配，不能匹配就沿着f指针转移就是了。

在这个过程中维护自动机的一个节点对应的right集合的最大匹配长度。取最大值就是答案。

 

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

#define N 500500

using namespace std;

struct SAM
{
    SAM *son[26],*f;
    int l;
}sam[N],*head,*last;

char str[N],b[N];
int cnt,ans,len;

inline void add(int x)
{
    SAM *p=&sam[++cnt],*bj=last;
    p->l=last->l+1; last=p;
    for(;bj&&!bj->son[x];bj=bj->f) bj->son[x]=p;
    if(!bj) p->f=head;
    else if(bj->l+1==bj->son[x]->l) p->f=bj->son[x];
    else
    {
        SAM *r=&sam[++cnt],*q=bj->son[x];
        *r=*q; r->l=bj->l+1; 
        p->f=q->f=r;
        for(;bj&&bj->son[x]==q;bj=bj->f) bj->son[x]=r;
    }
}

inline void read()
{
    scanf("%s%s",str,b);
    len=strlen(str);
    head=last=&sam[cnt=0];
    for(int i=0;i<len;i++) add(str[i]-'a');
    len=strlen(b);
    for(int i=0;i<len;i++) b[i]-='a';
}

inline void go()
{
    last=head;
    int res=0;
    for(int i=0;i<len;i++)
    {
        if(last->son[b[i]]) res++,last=last->son[b[i]];
        else
        {
            for(;last&&!last->son[b[i]];last=last->f);
            if(!last) last=head,res=0;
            else res=last->l+1,last=last->son[b[i]];
        }
        ans=max(ans,res);
    }
    printf("%d\n",ans);
}

int main()
{
    read(),go();
    return 0;
}