[POJ 2774] 最長公共子串
阿新 • • 發佈:2019-01-05
題目大意:給定兩個字串,求它們的最長公共子串(廢話)。
字尾陣列大模板題。
兩串接起來,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;
}