后缀数组入门题——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;
}