后缀数组入门题——2323后缀排序

最新推荐文章于 2024-06-14 10:28:18 发布

原创最新推荐文章于 2024-06-14 10:28:18 发布 · 788 阅读

1 ·

CC 4.0 BY-SA版权

后缀数组专栏收录该内容

12 篇文章

订阅专栏

本文介绍了一种构建后缀数组的有效方法，通过详细解析代码实现过程，帮助读者理解后缀排序算法及其应用。该方法首先初始化单字符排序，然后逐步扩大排序范围直至整个字符串长度。

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

题目大意：给定一个字符串，输出排名为1到n的后缀的起始位置

思路：按后缀数组的sa的定义，直接输出即可

code：

#include<cstdio>
#include<cstring>
using namespace std;
const int maxl=200000;
int rank[maxl+10],sa[maxl+10],n,sum[maxl+10],trank[maxl+10],tsa[maxl+10];
char s[maxl+10];

void sorts(int j){
memset(sum,0,sizeof(sum));
for (int i=1;i<=n;i++) sum[rank[i+j]]++;
for (int i=1;i<=n;i++) sum[i]+=sum[i-1];
for (int i=n;i>0;i--) tsa[sum[rank[i+j]]--]=i;

memset(sum,0,sizeof(sum));
for (int i=1;i<=n;i++) sum[rank[tsa[i]]]++;
for (int i=1;i<=n;i++) sum[i]+=sum[i-1];
for (int i=n;i>0;i--) sa[sum[rank[tsa[i]]]--]=tsa[i];
}

void getsa(){
for (int i=1;i<=n;i++) trank[i]=s[i];
for (int i=1;i<=n;i++) sum[trank[i]]++;
for (int i=1;i<=255;i++) sum[i]+=sum[i-1];
for (int i=n;i>0;i--) sa[sum[trank[i]]--]=i;
rank[sa[1]]=1;
for (int i=2,p=1;i<=n;i++){
if (trank[sa[i]]!=trank[sa[i-1]]) p++;
rank[sa[i]]=p;
}
for (int j=1;j<=n;j*=2){
sorts(j);
trank[sa[1]]=1;
for (int i=2,p=1;i<=n;i++){
if (rank[sa[i]]!=rank[sa[i-1]]||rank[sa[i-1]+j]!=rank[sa[i]+j]) p++;
trank[sa[i]]=p;
}
memcpy(rank,trank,sizeof(rank));
}
}

int main(){
scanf("%s",s+1);
n=strlen(s+1);
getsa();
for (int i=1;i<=n;i++)
printf("%d ",sa[i]);
puts("");for (;;);
return 0;
}