bzoj1031 [JSOI2007]字符加密Cipher
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。
正解:后缀数组。
直接后缀排序求出sa数组,然后输出即可。然而后缀数组好难啊,只能背代码。。
1 #include <algorithm> 2 #include <iostream> 3 #include <iomanip> 4 #include <cstring> 5 #include <cstdlib> 6 #include <cstdio> 7 #include <cmath> 8 #include <string> 9 #define N (500010) 10 #define il inline 11 #define RG register 12 #define File(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout) 13 14 using namespace std; 15 16 int m=128,n,nn,wa[N],wb[N],sw[N],sa[N],rak[N]; 17 char s[N]; 18 19 int cmp(int *r,int a,int b,int l){ return r[a]==r[b]&&r[a+l]==r[b+l]; } 20 21 void SA(){ 22 RG int i,j,p;int *x=wa,*y=wb,*k; 23 for (i=0; i<n; i++) sw[x[i]=s[i]]++; 24 for (i=1; i<m; i++) sw[i]+=sw[i-1]; 25 for (i=n-1; i>=0; i--) sa[--sw[x[i]]]=i; 26 for (j=1; j<=n; j<<=1,m=p){ 27 for (p=0,i=n-j; i<n; i++) y[p++]=i; 28 for (i=0; i<n; i++) if (sa[i]>=j) y[p++]=sa[i]-j; 29 for (i=0; i<m; i++) sw[i]=0; 30 for (i=0; i<n; i++) sw[x[y[i]]]++; 31 for (i=1; i<m; i++) sw[i]+=sw[i-1]; 32 for (i=n-1; i>=0; i--) sa[--sw[x[y[i]]]]=y[i]; 33 for (k=x,x=y,y=k,p=1,x[sa[0]]=0,i=1; i<n; i++) 34 x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++; 35 if (p>=n) break; 36 } 37 return; 38 } 39 40 int main(){ 41 File("js"); scanf("%s",s); n=strlen(s),nn=n; 42 for (RG int i=n;i<(n<<1);++i) s[i]=s[i-n]; n<<=1; SA(); 43 for (RG int i=0;i<n;++i) if (sa[i]+nn<n) printf("%c",s[sa[i]+nn-1]); 44 return 0; 45 }