BZOJ_1031_[JSOI2007]_字符串加密_(后缀数组)

最新推荐文章于 2018-12-22 14:23:22 发布

aquz77529

最新推荐文章于 2018-12-22 14:23:22 发布

阅读量127

点赞数

文章标签： php

原文链接：http://www.cnblogs.com/Sunnie69/p/5615552.html

版权

描述

http://www.lydsy.com/JudgeOnline/problem.php?id=1031

长度为n的字符串形成环,枚举开头位置,得到n个长度为n的串,将这些串排序,然后按照顺序输出每一个串的最后一个字符.

分析

将原串复制一遍接在自己的后面,这样n个串都在总串里面出现了.

对于这n个串中的两个a,b,如果a<b,那么以a开头的后缀就<以b开头的后缀.

所以我们求一下后缀数组就好了.

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 
 4 const int maxn=200000+5;
 5 int n;
 6 int t1[maxn],t2[maxn],sa[maxn],c[maxn];
 7 char s[maxn];
 8 
 9 void build_sa(char *s,int n,int m){
10     int i,k,p,*x=t1,*y=t2,*t;
11     for(i=0;i<m;i++) c[i]=0;
12     for(i=0;i<n;i++) c[x[i]=s[i]-'\0']++;
13     for(i=1;i<m;i++) c[i]+=c[i-1];
14     for(i=0;i<n;i++) sa[--c[x[i]]]=i;//
15     for(k=1,p=1;p<n;k<<=1,m=p){
16         for(p=0,i=n-k;i<n;i++) y[p++]=i;//
17         for(i=0;i<n;i++)if(sa[i]>=k) y[p++]=sa[i]-k;
18         for(i=0;i<m;i++) c[i]=0;
19         for(i=0;i<n;i++) c[x[i]]++;//
20         for(i=1;i<m;i++) c[i]+=c[i-1];
21         for(i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
22         for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
23             x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?p-1:p++;
24     }
25 }
26 int main(){
27     scanf("%s",s+1); n=strlen(s+1);
28     for(int i=1;i<=n;i++) s[n+i]=s[i];
29     s[0]='\0';
30     build_sa(s,n*2+1,maxn);
31     for(int i=1;i<=2*n;i++) if(sa[i]<=n) printf("%c",s[sa[i]+n-1]);
32     return 0;
33 }