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

解题思路：

sa裸题，将字符串倍长求sa

代码：

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 const int N=1000000;
 5 int rnk[N];
 6 int tmr[N];
 7 int has[N];
 8 int sa[N];
 9 int ln[N];
10 char a[N];
11 int n;
12 int len;
13 int cnt;
14 bool Same(int a,int b,int l)
15 {
16     if(a+l>len||b+l>len)
17         return false;
18     return (rnk[a+l]==rnk[b+l])&&(rnk[a]==rnk[b]);
19 }
20 int main()
21 {
22     scanf("%s",a+1);
23     n=strlen(a+1);
24     len=2*n;
25     for(int i=1;i<=n;i++)
26         ln[i]=ln[i+n]=a[i];
27     for(int i=1;i<=len;i++)
28         has[ln[i]]++;
29     for(int i=0;i<128;i++)
30         if(has[i])
31             tmr[i]=++cnt;
32     for(int i=1;i<128;i++)
33         has[i]+=has[i-1];
34     for(int i=1;i<=len;i++)
35     {
36         rnk[i]=tmr[ln[i]];
37         sa[has[ln[i]]--]=i;
38     }
39     for(int k=1;cnt!=len;k<<=1)
40     {
41         cnt=0;
42         for(int i=0;i<=len;i++)
43             has[i]=0;
44         for(int i=1;i<=len;i++)
45             has[rnk[i]]++;
46         for(int i=1;i<=len;i++)
47             has[i]+=has[i-1];
48         for(int i=len;i;i--)
49             if(sa[i]>k)
50                 tmr[sa[i]-k]=has[rnk[sa[i]-k]]--;
51         for(int i=1;i<=k;i++)
52             tmr[len-i+1]=has[rnk[len-i+1]]--;
53         for(int i=1;i<=len;i++)
54             sa[tmr[i]]=i;
55         for(int i=1;i<=len;i++)
56             if(Same(sa[i],sa[i-1],k))
57                 tmr[sa[i]]=cnt;
58             else
59                 tmr[sa[i]]=++cnt;
60         for(int i=1;i<=len;i++)
61             rnk[i]=tmr[i];
62     }
63     for(int i=1;i<=len;i++)
64         if(sa[i]<=n)
65             printf("%c",a[(sa[i]-2+n)%n+1]);
66     puts("");
67     return 0;
68 }