3160 最长公共子串

时间限制: 2 s

空间限制: 128000 KB

题目等级 : 大师 Master

题目描述 Description

给出两个由小写字母组成的字符串，求它们的最长公共子串的长度。

输入描述 Input Description

读入两个字符串

输出描述 Output Description

输出最长公共子串的长度

样例输入 Sample Input

yeshowmuchiloveyoumydearmotherreallyicannotbelieveit

yeaphowmuchiloveyoumydearmother

样例输出 Sample Output

27

数据范围及提示 Data Size & Hint

单个字符串的长度不超过100000

分类标签 Tags

后缀树 后缀树组 树结构

/*

后缀数组求lcp裸题.

先用无关字符将两个串连接起来.

然后用后缀数组求一下rank.

这样答案显然就是两个相邻rank后缀的最大前缀.

这样的话.

求出height数组,

然后找一个i和i-1分别属于前后两个串就可以了.

*/

#include<iostream>

#include<cstring>

#include<cstdio>

#define MAXN 2000010

using namespace std;

int n,sa[MAXN],l1,l2,m=128,rank[MAXN],height[MAXN],t1[MAXN],c[MAXN],t2[MAXN],ans,tot,s[MAXN];

char ch1[MAXN],ch2[MAXN],ch[MAXN];

bool cmp(int *y,int a,int b,int k)

{

    int a1=y[a],b1=y[b];

    int a2=a+k>=n?-1:y[a+k];

    int b2=b+k>=n?-1:y[b+k];

    return a1==b1&&a2==b2;

}

void slovesa()

{

    int *x=t1,*y=t2;

    for(int i=0;i<m;i++) c[i]=0;

    for(int i=0;i<n;i++) c[x[i]=s[i]]++;

    for(int i=1;i<m;i++) c[i]+=c[i-1];

    for(int i=n-1;i>=0;i--) sa[--c[x[i]]]=i;

    for(int k=1,p=0;k<=n;k<<=1,m=p,p=0)

    {

        for(int i=n-k;i<n;i++) y[p++]=i;

        for(int i=0;i<n;i++) if(sa[i]>=k) y[p++]=sa[i]-k;

        for(int i=0;i<m;i++) c[i]=0;

        for(int i=0;i<n;i++) c[x[y[i]]]++;

        for(int i=1;i<m;i++) c[i]+=c[i-1];

        for(int i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];

        swap(x,y);p=1;x[sa[0]]=0;

        for(int i=1;i<n;i++)

        {

            if(cmp(y,sa[i-1],sa[i],k)) x[sa[i]]=p-1;

            else x[sa[i]]=p++;

        }

        if(p>=n) break;

    }

}

void sloveheight()

{

    int k=0;

    for(int i=0;i<n;i++) rank[sa[i]]=i;

    for(int i=0;i<n;height[rank[i++]]=k)

    {

        if(!rank[i]) continue;

        int j=sa[rank[i]-1];

        if(k) k--;

        while(j+k<n&&i+k<n&&s[j+k]==s[i+k]) k++;

    }

    height[0]=0;

}

void sloveans()

{

    for(int i=0;i<n;i++)

      if(height[i]>ans)

      {

        if(sa[i]<l1&&sa[i-1]>l1) ans=height[i];

        if(sa[i]>l1&&sa[i-1]<l1) ans=height[i];

      }

}

int main()

{

    scanf("%s",ch1);

    scanf("%s",ch2);

    l1=strlen(ch1);

    l2=strlen(ch2);

    for(int i=0;i<l1;i++) s[n++]=ch1[i];

    s[n++]='$';

    for(int i=0;i<l2;i++) s[n++]=ch2[i];

    slovesa(),sloveheight(),sloveans();

    printf("%d",ans);

    return 0;

}