SPOJ 1811 Longest Common Substring(字尾自動機)

題目連結：http://www.spoj.com/problems/LCS/

題意：求兩個串的最長公共substring。

思路：建立一個串的SAM，對於第二個串的每個字首，找到其在第一個串上匹配的最大長度。

const int KIND=26;
struct SAM
{
    SAM *son[KIND],*pre;
    int len,cnt;

    void init()
    {
        clr(son,NULL);
        pre=NULL;
        len=0;
    }
};

SAM sam[N],*head,*last,*b[N];
int d[N],cnt,f[N],len;
char s[N];

void initSAM()
{
    sam[0].init();
    head=last=&sam[0];
    cnt=1;
}


void insert(int x)
{
    SAM *p=&sam[cnt++],*u=last;
    p->len=last->len+1;
    last=p;
    for(;u&&!u->son[x];u=u->pre) u->son[x]=p;
    if(!u) p->pre=head;
    else if(u->son[x]->len==u->len+1) p->pre=u->son[x];
    else
    {
        SAM *r=&sam[cnt++],*q=u->son[x];
        *r=*q; r->len=u->len+1;
        p->pre=q->pre=r;
        for(;u&&u->son[x]==q;u=u->pre) u->son[x]=r;
    }
}

int main()
{
    RD(s);
    len=strlen(s);
    initSAM();
    int i;
    FOR0(i,len) insert(s[i]-'a');
    RD(s); len=strlen(s);
    int ans=0,x,L=0;
    SAM *p=head;
    FOR0(i,len)
    {
        x=s[i]-'a';
        if(p->son[x]) L++,p=p->son[x];
        else
        {
            while(p&&!p->son[x]) p=p->pre;
            if(p) L=p->len+1,p=p->son[x];
            else p=head,L=0;
        }
        ans=max(ans,L);
    }
    PR(ans);
    return 0;
}