題目鏈接：點擊查看

題目大意：給出兩個字符串，求出其最長的公共子串

題目分析：后綴數(shù)組水題，直接用一個不同的符號拼接一下兩個字符串，然后跑出height數(shù)組和sa數(shù)組，再遍歷一遍取最大值就是答案了，需要注意的一個地方就是，若想讓height對答案作出貢獻(xiàn)，必須滿足sa數(shù)組表示的兩個數(shù)在不同符號的兩側(cè)才行

代碼：

#include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cstring> #include<algorithm> #include<stack> #include<queue> #include<map> #include<set> #include<sstream> using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=2e5+100;char s1[N],s2[N];int sa[N]; //SA數(shù)組，表示將S的n個后綴從小到大排序后把排好序的 //的后綴的開頭位置順次放入SA中 int t1[N],t2[N],c[N];int rk[N],height[N],n;int s[N];void build_sa(int s[],int n,int m)//n為添加0后的總長 {int i,j,p,*x=t1,*y=t2;for(i=0;i<m;i++) c[i]=0;for(i=0;i<n;i++) c[x[i]=s[i]]++;for(i=1;i<m;i++) c[i]+=c[i-1];for(i=n-1;i>=0;i--) sa[--c[x[i]]]=i;for(j=1;j<=n;j<<=1) {p=0;for(i=n-j;i<n;i++) y[p++]=i;for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;for(i=0;i<m;i++) c[i]=0;for(i=0;i<n;i++) c[x[y[i]]]++;for(i=1;i<m;i++) c[i]+=c[i-1];for(i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];swap(x,y);p=1,x[sa[0]]=0;for(i=1;i<n;i++) x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+j]==y[sa[i]+j]?p-1:p++;if(p>=n) break;m=p;} }void get_height(int s[],int n)//n為添加0后的總長 {int i,j,k=0;for(i=0;i<=n;i++)rk[sa[i]]=i;for(i=0;i<n;i++) {if(k) k--;j=sa[rk[i]-1];while(s[i+k]==s[j+k]) k++;height[rk[i]]=k;} }int main() { // freopen("input.txt","r",stdin); // ios::sync_with_stdio(false);scanf("%s%s",s1,s2);for(int i=0;s1[i];i++)s[i]=s1[i]-'a'+1;int mark=strlen(s1);s[mark]=30;for(int i=0;s2[i];i++)s[i+mark+1]=s2[i]-'a'+1;n=strlen(s1)+strlen(s2)+1;s[n]=0;build_sa(s,n+1,31);get_height(s,n);int ans=0;for(int i=2;i<=n-1;i++)if(sa[i-1]<mark&&sa[i]>mark||sa[i-1]>mark&&sa[i]<mark)ans=max(ans,height[i]);printf("%d\n",ans);return 0; }

?

總結(jié)

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

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