SPOJ LCS --后缀自动机

【题目大意】

        构造第一个串的后缀自动机,用第二个串与其逐位匹配。

        记录当前匹配到的位置p和位数len,如果下一个字符c仍能匹配成功,则更新匹配位置p=p.go[c]、匹配位数len++,否则沿后缀链接跳转至第一个能够接受c的位置q,令len=q.step。


【代码】

/************************
    ID:Ciocio
    LANG:C++
    DATE:2014-1-28
    TASK:LCS
************************/
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>

using namespace std;

#define MAXN 250010
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define rrep(i,a,b) for(int i=a;i>=b;--i)

char A[MAXN],B[MAXN];
struct SAM{
    int tot,last;
    int go[MAXN<<1][26],pre[MAXN<<1],step[MAXN<<1];

    void new_node(int s){
        step[++tot]=s;
        pre[tot]=0;
        memset(go[tot],0,sizeof go[tot]);
    }

    void build(char* s){
        tot=0;last=1;
        new_node(0);
        int n=strlen(s);
        rep(i,0,n-1){
            new_node(step[last]+1);
            int c=s[i]-'a';
            int p=last,np=tot,q,nq;
            for(;p&&!go[p][c];p=pre[p]) go[p][c]=np;
            if(!p)
                pre[np]=1;
            else{
                q=go[p][c];
                if(step[q]==step[p]+1) pre[np]=q;
                else{
                    new_node(step[p]+1);
                    nq=tot;
                    memcpy(go[nq],go[q],sizeof go[q]);
                    pre[nq]=pre[q];
                    pre[np]=
                    pre[q]=nq;
                    for(;p&&go[p][c]==q;p=pre[p]) go[p][c]=nq;
                }
            }
            last=np;
        }
    }

    int match(char* s){
        int p=1,len=0,res=0;
        int n=strlen(s);
        rep(i,0,n-1){
            int c=s[i]-'a';
            while(pre[p]&&!go[p][c]){
                p=pre[p];
                len=step[p];
            }
            if(go[p][c]){
                p=go[p][c];
                len++;
            }
            res=max(res,len);
        }
        return res;
    }
};

SAM run;
int main(){
    scanf("%s%s",A,B);
    run.build(A);
    cout<<run.match(B)<<endl;
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值