來自蒟蒻XXJ的做題記錄

今天蒟蒻XXJ學習了SAM =w=

于是找了一道板子寫

題目大意：最長公共子串

首先對A串搞一個自動機

然后用B串在自動機上面走

突然發(fā)現(xiàn)了這樣一件事情的我內(nèi)心是崩潰的：

1.trans(p,c)表示的應該是可以直接把c接在p狀態(tài)尾部達到的狀態(tài)。

2.在SAM上面走的時候其實和AC自動機是差不多的 它的fa數(shù)組和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; }

轉(zhuǎn)載于:https://www.cnblogs.com/Xiaojianxiang/p/6596599.html

總結

以上是生活随笔為你收集整理的SPOJ1811的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。