后缀数组模板

最新推荐文章于 2026-06-20 16:06:50 发布

原创最新推荐文章于 2026-06-20 16:06:50 发布 · 481 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#c语言

后缀数组同时被 3 个专栏收录

3 篇文章

订阅专栏

最长公共前缀

1 篇文章

订阅专栏

RMQ

1 篇文章

订阅专栏

本文介绍了一种构建后缀数组SA的方法，并实现了最长公共前缀LCP的计算。通过递归地划分字符串，算法有效地构建了后缀数组，并利用RMQ查询优化了LCP的计算过程。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
void build_sa(int M){
    for(int i=1;i<=M;i++) c[i]=0;
    for(int i=1;i<=N;i++) c[x[i]=s[i]]++;
    for(int i=2;i<=M;i++) c[i]+=c[i-1];
    for(int i=N;i>=1;i--) sa[c[x[i]]--]=i;
    for(int k=1;k<=n;k<<=1){
        int p=1;
        for(int i=n-k+1;i<=N;i++) y[++p]=i;
        for(int i=1;i<=N;i++) if(sa[i]>k) y[++p]=sa[i]-k;
        for(int i=1;i<=M;i++) c[i]=0;
        for(int i=1;i<=N;i++) c[x[y[i]]]++;
        for(int i=2;i<=M;i++) c[i]+=c[i-1];
        for(int i=N;i>=1;i--) sa[c[x[y[i]]]--]=y[i];
        p=0;swap(x,y);
        for(int 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:++p;
        }
        if(p==N) break;
        M=l;
    }
    memcpy(Rank,x,sizeof(x));
    int k=0;
    for(int I=1;i<=n;i++){
        if(k) k--;
        int j=sa[Rank[i]+1];
        if(!j) continue;
        while(s[i+k]==s[j+k]) k++;
        height[Rank[i]]=k;
    }
    for(int i=1;i<=n;i++) Min[i][0]=height[i];
    for(int j=1;j<=18;j++)
        for(int i=1;(i+(1<<j))<=n;i++)
            Min[i][j]=min(Min[i][j-1],Min[i+(1<<j)][j-1]);
}
int RMQ(int l,int r){//查询l，r两个后缀的LCP
    l=Rank[l],r=Rank[r];
    if(l>r) swap(l,r);
    int t=Log[r-l];
    return min(Min[l][k],Min[r-(1<<k)+1][k]);
}
int main(){
    scanf("%d%d",&n,&m);
    return 0;
}

^_^