Poj 1743 后缀数组

最新推荐文章于 2020-04-03 15:50:00 发布

MrBird_to_fly

最新推荐文章于 2020-04-03 15:50:00 发布

阅读量504

点赞数

CC 4.0 BY-SA版权

分类专栏：字符串

本文链接：https://blog.youkuaiyun.com/MrBird_to_fly/article/details/59110028

字符串专栏收录该内容

11 篇文章

订阅专栏

本文详细介绍了后缀数组的构建过程及其实现代码，并通过一个具体示例进行讲解。作者针对初学者可能遇到的问题进行了说明，并对lrj的白书中缺乏的一些细节给出了补充解释。

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

第一次写后缀数组，熟练一下模板。纠结了很久正负1的问题，大概是这样的：字符串从0开始，长度为n，在后面加一个0，这个0在build_sa里是要算上的，算rank的时候也要算上的，但是显然这个最后0的rank是0，所以为了简写，直接在算height的时候略过了，但我觉得直接统一一下，所有的数组都算上，然后加个if特判一下，这样比较好记也不太容易搞混。然后吐槽下lrj的白书，缺少了很多细节的解释，对初学者实在是太不友好了。

这里主要是为了整理一个比较好的板子，关于这个题的本身倒不是特别重要了。

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=20500;
int s[maxn];
int sa[maxn], t[maxn], t2[maxn], c[maxn], n;
int Rank[maxn], height[maxn];
void build_sa(int m, int n){
    int i, *x=t, *y=t2;
    for(i=0;i<m;i++)c[i]=0;
    for(i=0;i<n;i++)c[x[i]=s[i]]++;
    for(i=1;i<m;i++)c[i]+=c[i-1];
    for(i=n-1;i>=0;i--)sa[--c[x[i]]]=i;
    for(int k=1;k<=n;k<<=1){
        int p=0;
        for(i=n-k;i<n;i++)y[p++]=i;
        for(i=0;i<n;i++)if(sa[i]>=k)y[p++]=sa[i]-k;
        for(i=0;i<m;i++)c[i]=0;
        for(i=0;i<n;i++)c[x[y[i]]]++;
        for(i=1;i<m;i++)c[i]+=c[i-1];
        for(i=n-1;i>=0;i--)sa[--c[x[y[i]]]]=y[i];
        swap(x, y);
        p=1;x[sa[0]]=0;
        for(i=1;i<n;i++)
            x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
        if(p>=n)break;
        m=p;
    }
}
void geth(int n){
    int i, j, k=0;
    for(i=0;i<n;i++)Rank[sa[i]]=i;
    for(i=0;i<n;i++){
        if(k)k--;
        if(Rank[i]==0)break;
        j=sa[Rank[i]-1];
        while(s[i+k]==s[j+k])k++;
        height[Rank[i]]=k;
    }
}
bool check(int mid, int n){
    int i, j;
    int st=0;
    int mi, ma;
    for(i=1;i<n;i++){
        if(height[i]<mid){
            mi=n+1, ma=-1;
            for(j=st;j<i;j++){
                mi=min(mi, sa[j]);
                ma=max(ma, sa[j]);
            }
            if(ma-mi>mid)return true;
            st=i;
        }
    }
    mi=n+1, ma=-1;
    for(i=st;i<n;i++){
        mi=min(mi, sa[i]);
        ma=max(ma, sa[i]);
    }
    if(ma-mi>mid)return true;
    return false;
}
int main()
{
    while(scanf("%d", &n), n){
        for(int i=0;i<n;i++)scanf("%d", &s[i]);
        for(int i=0;i<n-1;i++)s[i]=s[i+1]-s[i]+100;
        s[n-1]=0;
        build_sa(200, n);
        geth(n);
        int l=3, r=n-1, ans=0;
        while(l<=r){
            if(r-l<=1){
                if(check(r, n))ans=r;
                else ans=l;
                break;
            }
            int mid=(l+r)/2;
            if(check(mid, n))l=mid;
            else r=mid;
        }
        if(ans<4)printf("0\n");
        else printf("%d\n", ans+1);
    }
}