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

题目大意:给定一个字符串,输出排名为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;
}


posted @ 2014-12-03 11:49  orzpps  阅读(139)  评论(0编辑  收藏  举报