【BZOJ 1031】 [JSOI2007]字符加密Cipher|后缀数组

先把串复制一遍

然后对这个2n的串求后缀数组

然后输出sa值小于n的


开始我只是对n做后缀数组 wa不停 耽误了很长时间

如果只是对n求后缀数组的话 不能求出后面n的rank值 就不能倍增!!!


#include <cstring>
#include <cstdio>
#include <iostream>
using namespace std;
const int MAXN=200000+100;
char s[MAXN];
int n,sa[MAXN],rank[MAXN],tmp[MAXN],cnt[MAXN],sec[MAXN];
void get_sa()
{
	for(int i=1;i<=n;i++) cnt[s[i]]++;
	for(int i=1;i<=200;i++) cnt[i]+=cnt[i-1];
	for(int i=1;i<=n;i++) sa[cnt[s[i]]--]=i;
	int nn=rank[sa[1]]=1;
	for(int i=2;i<=n;i++)
	{
		if(s[sa[i]]!=s[sa[i-1]]) nn++;
		rank[sa[i]]=nn;
	}
	for(int i=0;i<=200;i++) cnt[i]=0;
	for(int k=1;k<=n&&nn<n;k<<=1)
	{
		int top=0;
		for(int i=n-k+1;i<=n;i++) sec[++top]=i;
		for(int i=1;i<=n;i++) if(sa[i]>k) sec[++top]=sa[i]-k;

		for(int i=0;i<=nn;i++) cnt[i]=0;
		for(int i=1;i<=n;i++) cnt[rank[i]]++;
		for(int i=1;i<=nn;i++) cnt[i]+=cnt[i-1];
		for(int i=n;i>=1;i--) sa[cnt[rank[sec[i]]]--]=sec[i];
		top=tmp[sa[1]]=1;
		for(int i=2;i<=n;i++)
		{
			if(rank[sa[i]]!=rank[sa[i-1]]||rank[sa[i]+k]!=rank[sa[i-1]+k]) top++;
			tmp[sa[i]]=top;
		}
		nn=top;
		for(int i=1;i<=n;i++) rank[i]=tmp[i];
	}

}
int main()
{
	scanf("%s",s+1);
	n=strlen(s+1);int t=n;
	for(int i=1;i<=n;i++) s[i+n]=s[i];
	n*=2;get_sa();n=t;
	for(int i=1;i<=n*2;i++) if(sa[i]<=n)printf("%c",s[sa[i]+n-1]);printf("\n");
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值