題目鏈接 POJ2774
SPOJ1811 LCS - Longest Common Substring
比后綴自動機慢好多(廢話→_→)。

\(Description\)

求兩個字符串最長公共子串

\(Solution\)

任何一個子串一定是某個后綴的前綴
可以將兩個字符串拼在一起，中間用一個從未出現(xiàn)過的字符隔開，這樣ht[]的最大值就是答案？
不一定，最大的ht[]可能是由同一個字符串得到的，判一下屬于哪個字符串即可

//3772K 516MS //SPOJ:26M 0.11s(N=5e5) #include <cstdio> #include <cstring> #include <algorithm> const int N=2e5+10;int n,l,sa[N],rk[N],ht[N],sa2[N],tm[N]; char s[N];void Get_SA() {int *x=rk,*y=sa2,m=30;for(int i=1; i<=n; ++i) ++tm[x[i]=s[i]-'a'+1];for(int i=1; i<=m; ++i) tm[i]+=tm[i-1];for(int i=n; i; --i) sa[tm[x[i]]--]=i;for(int p=0,k=1; k<n; m=p,p=0,k<<=1){for(int i=n-k+1; i<=n; ++i) y[++p]=i;for(int i=1; i<=n; ++i) if(sa[i]>k) y[++p]=sa[i]-k;for(int i=0; i<=m; ++i) tm[i]=0;for(int i=1; i<=n; ++i) ++tm[x[i]];for(int i=1; i<=m; ++i) tm[i]+=tm[i-1];for(int i=n; i; --i) sa[tm[x[y[i]]]--]=y[i];std::swap(x,y), p=x[sa[1]]=1;for(int i=2; i<=n; ++i)x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p:++p;if(p>=n) break;}for(int i=1; i<=n; ++i) rk[sa[i]]=i;ht[1]=0;for(int k=0,p,i=1; i<=n; ++i){if(rk[i]==1) continue;if(k) --k;p=sa[rk[i]-1];while(i+k<=n&&p+k<=n&&s[i+k]==s[p+k]) ++k;ht[rk[i]]=k;} }int main() {scanf("%s",s+1), l=strlen(s+1);s[l+1]='z'+1;scanf("%s",s+2+l), n=strlen(s+1);Get_SA();int res=0;for(int i=2; i<=n; ++i)if((sa[i]<=l)^(sa[i-1]<=l)) res=std::max(res,ht[i]);printf("%d",res);return 0; }

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

總結(jié)

以上是生活随笔為你收集整理的POJ.2774.Long Long Message/SPOJ.1811.LCS(后缀数组 倍增)的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

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