求助
查看原帖
求助
128870
chen_qian楼主2021/3/9 11:56
#include<bits/stdc++.h>
#define N 1000006
using namespace std;
char s[N];
int px[N],cnt[N],id[N],sa[N],rk[N],oldrk[N<<1];
bool cmp(int x,int y,int w){
    return oldrk[x]==oldrk[y]&&oldrk[x+w]==oldrk[y+w];
}
int main(){
	scanf("%s",s+1);
	int n=strlen(s+1);
	int m=300;
	for(int i=1;i<=n;i++) cnt[rk[i]=s[i]]++;
	for(int i=1;i<=m;i++) cnt[i]+=cnt[i-1]; 
	for(int i=n;i>=1;i--) sa[cnt[rk[i]]--]=i;
	for(int w=1;w<n;w<<=1){
		int p=0;
		for(int i=n-w+1;i<=n;i++) id[++p]=i;
		for(int i=1;i<=n;i++)
			if(sa[i]>w) id[++p]=sa[i]-w;
		memset(cnt,0,sizeof(cnt));
		for(int i=1;i<=n;i++) cnt[rk[id[i]]]++;
		for(int i=1;i<=m;i++) cnt[i]+=cnt[i-1];
		for(int i=n;i>=1;i--) sa[cnt[rk[id[i]]--]=id[i];
		memcpy(oldrk,rk,sizeof(rk));
		for(int i=1;i<=n;i++)
			rk[sa[i]]=cmp(sa[i],sa[i-1],w)? p:++p;
	}
	for(int i=1;i<=n;i++) printf("%d ",sa[i]);
	return 0;
}

RT,对比 OIwiki 上的代码,感觉没有问题。

但在 luogu 上只能过一个点,求助

2021/3/9 11:56
加载中...