后缀数组的使用心得——POJ2774 最长连续公共子串

对于这道题，将两个字符串直接合并成为一个字符串，分别记录连个字符串结束的位置。

首先，应用黑暗圣典的模板，我们可以顺利得到height,rank,sa三个数组。

之后直接扫描1-n所有的位置，选出来一个，符合“两者都在不同子串的最大长度即可”。

此时我们会发现，sa数组记录了每个子串开头的位置，可以用于判断。

#include<iostream>
#include<stdio.h>
#include<string>
#include<string.h>
#include<vector>
#include<algorithm>
using namespace std;

const long long MAXN=;
const long long INF=1e8+;
char s[MAXN];
int sa[MAXN],t[MAXN],t2[MAXN],c[MAXN],n;
int rankk[MAXN],height[MAXN];
void getHeight()
{
    int k=;
    for(int i=;i<n;++i)rankk[sa[i]]=i;
    for(int i=;i<n;++i)
    {
        if(k)k--;
        int j=sa[rankk[i]-];
        while(s[i+k]==s[j+k])k++;
        height[rankk[i]]=k;
    }
} 

void build_sa(int m)
{
    int i,*x=t,*y=t2;

    for(int i=;i<m;++i)c[i]=;
    for(int i=;i<n;++i)c[x[i]=s[i]]++;
    for(int i=;i<m;++i)c[i]+=c[i-];
    for(int i=n-;i>=;--i)sa[--c[x[i]]]=i;
    for(int k=;k<=n;k*=)
    {
        int p=;
        for(int i=n-k;i<n;++i)y[p++]=i;
        for(int i=;i<n;++i)if(sa[i]>=k)y[p++]=sa[i]-k;

        for(int i=;i<m;++i)c[i]=;
        for(int i=;i<n;++i)c[x[y[i]]]++;
        for(int i=;i<m;++i)c[i]+=c[i-];
        for(int i=n-;i>=;i--)sa[--c[x[y[i]]]]=y[i];
        swap(x,y);
        p=;x[sa[]]=;
        for(int i=;i<n;++i)
        {
            x[sa[i]]=y[sa[i-]]==y[sa[i]]&&y[sa[i-]+k]==y[sa[i]+k]?p-:p++;

        }
        if(p>=n)break;
        m=p;
    }
}

int main()
{
//    cin.sync_with_stdio(false);
    scanf("%s",s);
    int len=strlen(s);
    scanf("%s",s+len);
    int len2=strlen(s);
    n=len2+;
    build_sa();
    getHeight();
    int maxx=-;
    for(int i=;i<n;++i)
    {
//        cout<<height[i]<<ends<<s+sa[i]<<endl;
        if((sa[i]>=len&&sa[i-]<len)||(sa[i]<len&&sa[i-]>=len))maxx=max(maxx,height[i]);
    }
    cout<<maxx<<"\n";
//    cout<<s<<endl;
//    cout<<len<<ends<<len2<<endl;
//
    return ;
}