【BZOJ 1031】[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。
后缀数组模板题,我现在也只会比这模板打模板题,后缀数组依旧晕着
1 #include<cstdio> 2 #include<cstring> 3 #define N 200001 4 using namespace std; 5 char ch[N]; 6 int a[N],n,k; 7 int v[N],sa[2][N],rk[2][N]; 8 void init(){ 9 scanf("%s",ch+1); 10 n=strlen(ch+1); 11 for (int i=1;i<=n;i++){ 12 a[i]=int(ch[i]); 13 ch[i+n]=ch[i]; 14 a[i+n]=a[i]; 15 } 16 n=n<<1; 17 } 18 19 void calsa(int sa[N],int rk[N],int SA[N],int RK[N]){ 20 for (int i=1;i<=n;i++)v[rk[sa[i]]]=i; 21 for (int i=n;i>=1;i--) 22 if (sa[i]>k) SA[v[rk[sa[i]-k]]--]=sa[i]-k; 23 for (int i=n-k+1;i<=n;i++) SA[v[rk[i]]--]=i; 24 for (int i=1;i<=n;i++) 25 RK[SA[i]]=RK[SA[i-1]]+(rk[SA[i-1]]!=rk[SA[i]]||rk[SA[i-1]+k]!=rk[SA[i]+k]); 26 27 } 28 29 void work(){ 30 int p=0,q=1; 31 for (int i=1;i<=n;i++)v[a[i]]++; 32 for (int i=1;i<=256;i++) v[i]+=v[i-1]; 33 for (int i=1;i<=n;i++)sa[p][v[a[i]]--]=i; 34 for (int i=1;i<=n;i++) 35 rk[p][sa[p][i]]=rk[p][sa[p][i-1]]+(a[sa[p][i-1]]!=a[sa[p][i]]); 36 k=1; 37 while (k<n){ 38 calsa(sa[p],rk[p],sa[q],rk[q]); 39 q^=1;p^=1;k=k<<1; 40 } 41 for (int i=1;i<=n;i++) if (sa[p][i]<=n/2) 42 printf("%c",ch[sa[p][i]+n/2-1]); 43 } 44 45 int main(){ 46 init(); 47 work(); 48 }