SP1811 【LCS - Longest Common Substring】
阿新 • • 發佈:2019-01-04
sca bst 中一 printf long long turn cst 最長 substr
\(SAM\)上匹配
我們就是需要找到兩個串的最長公共子串
先對其中一個串建出\(SAM\),之後我們把另一個串放到上面跑
如果當前在\(SAM\)的狀態是\(now\),下一個字符是\(c\),匹配出的的長度為\(L\)
如果\(now\)有\(c\)這個轉移,我們就轉移過去,\(L\)++
如果沒有我們就跳\(link\),知道跳到有這個轉移為止,同時把\(L\)搞成新狀態的\(len\)
這樣做就好了
代碼
#include<iostream> #include<cstring> #include<cstdio> #include<algorithm> #define LL long long #define maxn 500005 #define max(a,b) ((a)>(b)?(a):(b)) char S[maxn>>1],T[maxn>>1]; int lst=1,n,m,len[maxn],pre[maxn],son[maxn][26]; int now=1,L,ans,cnt=1; inline void ins(int c) { int f=lst,p=++cnt; lst=p; len[p]=len[f]+1; while(f&&!son[f][c]) {son[f][c]=p;f=pre[f];} if(!f) {pre[p]=1;return;} int x=son[f][c]; if(len[f]+1==len[x]) {pre[p]=x;return;} int y=++cnt; len[y]=len[f]+1;pre[y]=pre[x];pre[x]=pre[p]=y; for(int i=0;i<26;i++) son[y][i]=son[x][i]; while(f&&son[f][c]==x) {son[f][c]=y;f=pre[f];} } inline void q(int c) { if(son[now][c]) {now=son[now][c];L++;ans=max(ans,L);return;} while(now&&!son[now][c]) now=pre[now]; if(!now) {now=1,L=0;return;} L=len[now]+1;now=son[now][c];ans=max(ans,L); } int main() { scanf("%s",S+1);n=strlen(S+1); for(int i=1;i<=n;i++) ins((int)(S[i]-‘a‘)); scanf("%s",T+1);n=strlen(T+1); for(int i=1;i<=n;i++) q((int)(T[i]-‘a‘)); printf("%d\n",ans); return 0; }
SP1811 【LCS - Longest Common Substring】