原文链接http://www.cnblogs.com/zhouzhendong/p/8982392.html
题目传送门 - SPOJ LCS
题意
求两个字符串的最长公共连续子串长度。
字符串长$\leq 250000$
题解
首先对于第一个字符串建一个$SAM$。
然后拿第二个串在$SAM$上面走一遍就好了。
具体地:
将第二个串的字符一个一个地按照顺序加入。
设当前状态为$now$,要加入字符$c$,当前匹配的字符串长度为$len$(答案自然是各种情况下$len$的最大值)。
如果在$SAM$上面,状态$now$有标号为$c$的转移,那么,$len=len+1$,$now$更新为转移后的结果。
否则,我们跳$now$的$fa$,直到得到一个新的$now$使得$now$有标号为$c$的转移,并使$len=Max(now)+1$,$now$更新为新的$now$再走$c$转移之后的状态。
代码
#include <bits/stdc++.h> using namespace std; const int N=500005; int n,last=1,size=1; char s[N]; struct SAM{ int Next[26],fa,Max; }t[N]; void expend(int c){ int p=last,np=++size,q,nq; t[np].Max=t[p].Max+1; for (;!t[p].Next[c];p=t[p].fa) t[p].Next[c]=np; q=t[p].Next[c]; if (t[q].Max==t[p].Max+1) t[np].fa=q; else { nq=++size; t[nq]=t[q],t[nq].Max=t[p].Max+1; t[q].fa=t[np].fa=nq; for (;t[p].Next[c]==q;p=t[p].fa) t[p].Next[c]=nq; } last=np; } int main(){ t[0].Max=-1; for (int i=0;i<26;i++) t[0].Next[i]=1; scanf("%s",s); n=strlen(s); for (int i=0;i<n;i++) expend(s[i]-'a'); int ans=0; scanf("%s",s); n=strlen(s); for (int i=0,now=1,len=0;i<n;i++){ int c=s[i]-'a'; if (t[now].Next[c]){ now=t[now].Next[c]; ans=max(ans,++len); continue; } while (!t[now].Next[c]) now=t[now].fa; ans=max(ans,len=t[now].Max+1); now=t[now].Next[c]; } printf("%d",ans); return 0; }