[HZOI 2015] Glass Beads

题目描述】

给定长度为n(n<=300000)的循环同构的字符串,定义最小表示为该字符串的字典序最小的同构表示,请输出这个表示。

【输入格式】

第一行是串的长度,第二行是字符串。

【输出格式】

串的最小表示。

【样例输入】


10

helloworld


【样例输出】

dhelloworl

【题目来源】

HZOI2015 改编自poj1509

题解:

先将字符串*2

然后求出新字符串的后缀数组 找一个在后缀数组中下标最小的 位置小于原长的下标 那么从当前下标在字符串中位置开始提取长度为n的串就是原串的最小表示法

代码:

#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
const int maxn=600000+10;
int c[maxn],t[maxn],t2[maxn];
int sa[maxn];
char s[maxn];
int maxm=500;
int n;
inline bool cmp(int *y,int p,int q,int l){
	int o1=p+l>=n?-1:y[p+l];
	int o2=q+l>=n?-1:y[q+l];
	return o1==o2&&y[p]==y[q];
}
inline void build_sa(int m){
	int *x=t,*y=t2;
	for(int i=0;i<=m;i++)
		c[i]=0;
	for(int i=0;i<n;i++)
		c[x[i]=(s[i]-'a'+1)]++;
	for(int i=1;i<=m;i++)
		c[i]+=c[i-1];
	for(int i=n-1;i>=0;i--)
		sa[--c[x[i]]]=i;
	for(int k=1;k<=n;k<<=1){
		int p=0;
		for(int i=0;i<=m;i++)
			c[i]=0;
		for(int i=n-k;i<n;i++)
			y[p++]=i;
		for(int i=0;i<n;i++)
			if(sa[i]>=k)
				y[p++]=sa[i]-k;
		for(int i=0;i<n;i++)
			c[x[y[i]]]++;
		for(int i=1;i<=m;i++)
			c[i]+=c[i-1];
		for(int i=n-1;i>=0;i--)
			sa[--c[x[y[i]]]]=y[i];
		swap(x,y);
		x[sa[0]]=0;
		p=1;
		for(int i=1;i<n;i++)
			x[sa[i]]=cmp(y,sa[i],sa[i-1],k)?p-1:p++;
		if(p>=n)
			break;
		m=p;
	}
}
int main(){
	freopen("MinRepresentations.in","r",stdin);
	freopen("MinRepresentations.out","w",stdout);
	scanf("%d\n",&n);
	for(int i=0;i<n;i++){
		s[i]=getchar();
		while(s[i]<'a'||s[i]>'z') s[i]=getchar();
	}
	for(int i=0;i<n;i++)
		s[i+n]=s[i];
	n=n*2;
	build_sa(maxm);
	n/=2;
	int ans;
	for(int i=0;;i++)
		if(sa[i]<=n){
			ans=sa[i]%n;
			break;
		}
	for(int i=0;i<n;i++)
		printf("%c",s[ans+i]);
return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值