1. 程式人生 > >SP1811 【LCS - Longest Common Substring】

SP1811 【LCS - Longest Common Substring】

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】