SPOJ 1811 Longest Common Substring(后缀自动机)

偏执的太偏执、 2022-03-29 08:42 160阅读 0赞

题目链接：[http://www.spoj.com/problems/LCS/][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;
    }

[http_www.spoj.com_problems_LCS]: http://www.spoj.com/problems/LCS/