1031. [JSOI2007]字符加密【后缀数组】
Description
喜欢钻研问题的JS同学,最近又迷上了对加密方法的思考。一天,他突然想出了一种他认为是终极的加密办法
:把需要加密的信息排成一圈,显然,它们有很多种不同的读法。例如下图,可以读作:
JSOI07 SOI07J OI07JS I07JSO 07JSOI 7JSOI0把它们按照字符串的大小排序:07JSOI 7JSOI0 I07JSO JSOI07
OI07JS SOI07J读出最后一列字符:I0O7SJ,就是加密后的字符串(其实这个加密手段实在很容易破解,鉴于这是
突然想出来的,那就^^)。但是,如果想加密的字符串实在太长,你能写一个程序完成这个任务吗?
Input
输入文件包含一行,欲加密的字符串。注意字符串的内容不一定是字母、数字,也可以是符号等。
Output
输出一行,为加密后的字符串。
Sample Input
JSOI07
Sample Output
I0O7SJ
HINT
对于100%的数据字符串的长度不超过100000。
相同两个字符串连接然后裸后缀数组
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #define MAXN (200000+10) 5 using namespace std; 6 int wa[MAXN],wb[MAXN],wt[MAXN]; 7 int SA[MAXN],Rank[MAXN],Height[MAXN]; 8 int n,m; 9 char r[MAXN]; 10 11 inline bool cmp(int *y,int a,int b,int k) 12 { 13 int arank1=y[a]; 14 int brank1=y[b]; 15 int arank2=a+k>=n?-1:y[a+k]; 16 int brank2=b+k>=n?-1:y[b+k]; 17 return arank1==brank1 && arank2==brank2; 18 } 19 20 inline void Build_SA() 21 { 22 int *x=wa,*y=wb; 23 for (int i=0; i<m; ++i) wt[i]=0; 24 for (int i=0; i<n; ++i) wt[x[i]=r[i]]++; 25 for (int i=1; i<m; ++i) wt[i]+=wt[i-1]; 26 for (int i=n-1; i>=0; --i) SA[--wt[x[i]]]=i; 27 28 for (int j=1; j<=n; j<<=1) 29 { 30 int p=0; 31 for (int i=n-j; i<n; ++i) y[p++]=i; 32 for (int i=0; i<n; ++i) if (SA[i]>=j) y[p++]=SA[i]-j; 33 34 for (int i=0; i<m; ++i) wt[i]=0; 35 for (int i=0; i<n; ++i) wt[x[y[i]]]++; 36 37 for (int i=1; i<m; ++i) wt[i]+=wt[i-1]; 38 for (int i=n-1; i>=0; --i) SA[--wt[x[y[i]]]]=y[i]; 39 swap(x,y); 40 m=1; 41 x[SA[0]]=0; 42 for (int i=1; i<n; ++i) 43 x[SA[i]]=cmp(y,SA[i-1],SA[i],j)?m-1:m++; 44 if (m>=n) break; 45 } 46 } 47 48 int main() 49 { 50 scanf("%s",&r); 51 n=strlen(r); 52 for (int i=n; i<=n*2-1; ++i) r[i]=r[i-n]; 53 n*=2; 54 m=130; 55 Build_SA(); 56 for (int i=0; i<n; i++) 57 if (SA[i]<n/2) 58 printf("%c",r[SA[i]+n/2-1]); 59 }