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 }
posted @ 2018-03-30 22:37  Refun  阅读(153)  评论(0编辑  收藏  举报