poj 2774 Long Long Message-优快云博客

本文链接：https://blog.youkuaiyun.com/threeh20/article/details/78508481

本文介绍了一种使用后缀数组寻找两个字符串间最长公共子串的方法，并提供了详细的C++实现代码。通过构建后缀数组和高度数组，算法能够高效地找到最长公共子串的长度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

http://poj.org/problem?id=2774

后缀数组，求最长公共子串。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=222222;
char s[N];
int height[N],sa[N],rank[N],source[N],cnt[N],wa[N],wb[N];
int cmp(int *x,int a,int b,int l)
{
    return (x[a]==x[b])&&(x[a+l]==x[b+l]);
}
void suffix(int *s,int len)
{
    int *y=wb,*t,*r=wa;
    int i,j,k,n,p,m;
    n = len;
    m=128;
    for(i=0;i<m;i++)
        cnt[i]=0;
    for(i=0;i<n;i++)
        cnt[r[i]=s[i]]++;
    for(i=1;i<m;i++)
        cnt[i]=cnt[i]+cnt[i-1];
    for(i=n-1;i>=0;i--)
        sa[--cnt[r[i]]]=i;
    for(j=1,p=1;p<n;m=p,j*=2)
    {
        for(p=0,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++)
            cnt[i]=0;
    
        for(i=0;i<n;i++)
            cnt[r[y[i]]]++;
    
        for(i=1;i<m;i++)
            cnt[i]=cnt[i]+cnt[i-1];
    
        for(i=n-1;i>=0;i--)
            sa[--cnt[r[y[i]]]]=y[i];

        for(t=r,r=y,y=t,p=1,r[sa[0]]=0,i=1;i<n;i++)
            r[sa[i]]=cmp(y,sa[i],sa[i-1],j)?p-1:p++;
    
    }
    for(i=1;i<n;i++)
        rank[sa[i]]=i;
    k=0;
    for(i=0;i<n-1;i++)
    {
        j=sa[rank[i]-1];
        for(k?k--:0;s[i+k]==s[j+k];k++);
        height[rank[i]]=k;
    }
}
int main()
{
    int i,j,n,ans=0;
    scanf("%s",s);
    j=strlen(s);
    s[j]=1;
    scanf("%s",s+j+1);
    n=strlen(s);
    for(i=0;i<n;i++) source[i]=s[i];
    source[n]=0;
    suffix(source,n+1);
    int result=0;
    for(i=2;i<=n;i++)
    {
        if(result<height[i])
            if((sa[i]>j && sa[i-1]<j) || (sa[i]<j && sa[i-1]>j))
                result=height[i];
    }
    cout<<result<<endl;
    return 0;
}