比赛 清华集训2017模板练习 评测结果 AAAAAAAATT
题目名称 后缀排序 最终得分 80
用户昵称 再见 运行时间 1.827 s
代码语言 C++ 内存使用 16.89 MiB
提交时间 2017-07-17 17:29:35
显示代码纯文本
#include <cstdio>
#include <cstdlib>
#include <cstring>

int n,wa[2000010],wb[2000010],sw[1000010],sa[1000010];
char str[1000010],top,out[110];

bool cmp(int *a,int i,int j,int l){
	return a[i]==a[j]&&a[i+l]==a[j+l];
}

void init(int n,int m){
	int *x=wa,*y=wb,*t;
	for(int i=1;i<=n;i++) sw[x[i]=str[i]]++;
	for(int i=1;i<=m;i++) sw[i]+=sw[i-1];
	for(int i=n;i>=1;i--) sa[sw[x[i]]--]=i;
	for(int j=1,p=0;p<n;j<<=1,m=p){
		p=0; for(int i=n-j+1;i<=n;i++) y[++p]=i;
		for(int i=1;i<=n;i++) if(sa[i]>j) y[++p]=sa[i]-j;
		for(int i=1;i<=m;i++) sw[i]=0;
		for(int i=1;i<=n;i++) sw[x[y[i]]]++;
		for(int i=1;i<=m;i++) sw[i]+=sw[i-1];
		for(int i=n;i>=1;i--) sa[sw[x[y[i]]]--]=y[i];
		int i=2; for(p=1,t=x,x=y,y=t,x[sa[1]]=1;i<=n;i++)
			x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p:++p;
	}
}

void output(int x){
	if(!x) out[++top]=0;
	while(x) out[++top]=x%10,x/=10;
	while(top) putchar(out[top--]+'0');
	putchar(' ');
}

int main(){
	freopen("sais.in","r",stdin);
	freopen("sais.out","w",stdout);
	scanf("%s",str+1); n=strlen(str+1);
	init(n,256);
	for(int i=1;i<=n;i++)
		output(sa[i]-1);
	return 0;
}