后缀数组模板

最新推荐文章于 2024-09-12 08:56:38 发布

Frods

最新推荐文章于 2024-09-12 08:56:38 发布

阅读量416

点赞数

分类专栏：后缀数组最长公共前缀 RMQ 文章标签： c语言

本文链接：https://blog.csdn.net/frods/article/details/62285317

版权

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

3 篇文章 0 订阅

订阅专栏

最长公共前缀

1 篇文章 0 订阅

订阅专栏

RMQ

1 篇文章 0 订阅

订阅专栏

#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;
}