1. 程式人生 > >[POJ 2774] 最長公共子串

[POJ 2774] 最長公共子串

題目大意:給定兩個字串,求它們的最長公共子串(廢話)。

字尾陣列大模板題。
兩串接起來,height陣列求個最大的且分別在兩串中的就行了。
啥都不說貼程式碼。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <vector>
#include <queue>
#include <deque> #include <map> #include <set> #include <ctime> #define putarray(a,n) for (int i=1;i<=n-1;i++) printf("%d ",a[i]); printf("%d\n",a[n]); using namespace std; const int MAXN=100005; char A[MAXN],B[MAXN]; int ans,l1,l2; class Suffix_array { char str[MAXN<<1
]; int len,sa[MAXN<<1],rank[MAXN<<1],height[MAXN<<1],C[MAXN<<1]; void build_sa(int m) { int *x=rank,*y=height; //這裡rank,height只是拿來快取的。。。。 for (int i=0;i<=m;i++) C[i]=0; for (int i=1;i<=len;i++) x[i]=int(str[i]-96),C[x[i]]++; for
(int i=1;i<=m;i++) C[i]+=C[i-1]; for (int i=1;i<=len;i++) sa[C[x[i]]--]=i; for (int k=1;k<=len;k<<=1) { int p=0; for (int i=len-k+1;i<=len;i++) y[++p]=i; for (int i=1;i<=len;i++) if (sa[i]-k>0) y[++p]=sa[i]-k; for (int i=0;i<=m;i++) C[i]=0; for (int i=1;i<=len;i++) C[x[y[i]]]++; for (int i=1;i<=m;i++) C[i]+=C[i-1]; for (int i=len;i>=1;i--) sa[C[x[y[i]]]--]=y[i]; swap(x,y); p=1; x[sa[1]]=1; for (int i=2;i<=len;i++) if (y[sa[i]]==y[sa[i-1]] && y[sa[i]+k]==y[sa[i-1]+k]) x[sa[i]]=p;else x[sa[i]]=++p; if (p==len) break; m=p; } } void build_height() { int k=0; for (int i=1;i<=len;i++) rank[sa[i]]=i; for (int i=1;i<=len;i++) { if (k) k--; int j=sa[rank[i]-1]; while (str[i+k]==str[j+k]) k++; height[rank[i]]=k; } } public: void init() { int _max=0; l1=strlen(A),l2=strlen(B); len=l1+l2+1; for (int i=0;i<l1;i++) str[i+1]=A[i]; str[l1+1]=96; for (int i=0;i<l2;i++) str[i+l1+2]=B[i]; build_sa(26); build_height(); } void solve() { ans=0; for (int i=1;i<=l1;i++) if (sa[rank[i]-1]>l1+1) ans=max(ans,height[rank[i]]); for (int i=l1+1;i<=len;i++) if (sa[rank[i]-1]<=l1) ans=max(ans,height[rank[i]]); } }SA; int main() { //freopen("Longest common substring.in","r",stdin); //freopen("Longest common substring.out","w",stdout); scanf("%s%s",A,B); SA.init(); SA.solve(); printf("%d\n",ans); return 0; }