SPOJ1811

来自蒟蒻XXJ的做题记录

今天蒟蒻XXJ学习了SAM =w=

于是找了一道板子写

题目大意：最长公共子串

首先对A串搞一个自动机

然后用B串在自动机上面走

突然发现了这样一件事情的我内心是崩溃的：

1.trans(p,c)表示的应该是可以直接把c接在p状态尾部达到的状态。

2.在SAM上面走的时候其实和AC自动机是差不多的 它的fa数组和AC自动机中的fail指针很像。

B串在自动机上面走的时候 失配就利用fa向前找 看能不能当前配上的串里面找到一个后缀去配上这个字母进行trans

上代码，可以读读代码理解

#include<bits/stdc++.h>
#define mem(i,j) memset(i,j,sizeof(i))
#define mcy(i,j) memcpy(i,j,sizeof(i))
using namespace std;
const int MAXN=500010;
char a[MAXN],b[MAXN];
int len1,len2,ans;
struct SAM{
    int trans[MAXN][26],lth[MAXN],fa[MAXN];
    int cnt,last;
    SAM(){
        cnt=0;mem(trans,0);mem(fa,0);
        last=++cnt;
    }
    void extend(int c){
        int np=++cnt,p=last;last=cnt;
        for(;p&&!trans[p][c];p=fa[p]) trans[p][c]=np;
        if(!p) fa[np]=1;
        else{
            int q=trans[p][c];
            if(lth[q]==lth[p]+1) fa[np]=q;
            else{
                int nq=++cnt;lth[nq]=lth[p]+1;
                mcy(trans[nq],trans[q]);
                fa[nq]=fa[q];
                fa[np]=fa[q]=nq;
                for(;p&&trans[p][c]==q;p=fa[p]) trans[p][c]=nq;
            }
        }
    }
    void build(char str[]){
        int len=strlen(str);
        for(int i=0;i<len;i++) extend(str[i]-'a');
    }
}sam;
void input(){
    scanf("%s",a);
    scanf("%s",b);
    len1=strlen(a);len2=strlen(b);
    sam.build(a);
}
void xxj(){
    int len(0),p(1);
    for(int i=0;i<len2;i++){
        int c=b[i]-'a';
        if(sam.trans[p][c]){
            ++len;
            p=sam.trans[p][c];
        }
        else{
            while(p&&!sam.trans[p][c]) p=sam.fa[p];
            if(!p) p=1,len=0;
            else {len=sam.lth[p]+1;p=sam.trans[p][c];}
        }
        ans=max(ans,len);
    }
}
void output(){
    cout<<ans<<endl;
}

int main(){
    input();
    xxj();
    output();
    return 0;
}

转载于:https://www.cnblogs.com/Xiaojianxiang/p/6596599.html