bzoj1031[JSOI2007]字符加密Cipher（拆环变链，利用sa数组）

挺裸的一道后缀数组。。注意是环，把他在后面复制一遍再做就好了。。

#include <cstdio>
#include <cstring>
#define N 200005
char s[N];
int n,rank[N<<1],st[256],rank1[N],count[N],sa[N],tmp[N];
int main(){
//	freopen("a.in","r",stdin);
	scanf("%s",s+1);n=strlen(s+1);
	for(int i=1;i<=n;++i) s[n+i]=s[i];n<<=1;
	memset(st,0,sizeof(st));
	memset(rank,0,sizeof(rank));
	memset(rank1,0,sizeof(rank1));
	for(int i=1;i<=n;++i) st[s[i]]=1;
	for(int i=1;i<=255;++i) st[i]+=st[i-1];
	for(int i=n;i>=1;--i) rank[i]=st[s[i]];
	int k=1;
	for(int p=1;k!=n;p<<=1){
		memset(count,0,sizeof(count));
		for(int i=1;i<=n;++i) count[rank[i+p]]++;
		for(int i=1;i<=n;++i) count[i]+=count[i-1];
		for(int i=n;i>=1;--i) tmp[count[rank[i+p]]--]=i;
		memset(count,0,sizeof(count));
		for(int i=1;i<=n;++i) count[rank[tmp[i]]]++;
		for(int i=1;i<=n;++i) count[i]+=count[i-1];
		for(int i=n;i>=1;--i) sa[count[rank[tmp[i]]]--]=tmp[i];
		memcpy(rank1,rank,sizeof(rank1));
		rank[sa[1]]=k=1;
		for(int i=2;i<=n;++i){
			if(rank1[sa[i]]!=rank1[sa[i-1]]||rank1[sa[i]+p]!=rank1[sa[i-1]+p]) ++k;
			rank[sa[i]]=k;
		}
	}
	for(int i=1;i<=n;++i){
		if(sa[i]<=(n>>1)) printf("%c",s[sa[i]+n/2-1]);
	}
	return 0;
}