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;

}