后缀数组（bzoj 1031: [JSOI2007]字符加密Cipher）

后缀数组主要功能：

长度为n的字符串总共有n个后缀，求这n个后缀的字典序

实现方法：倍增+基数排序，过程就是下面那张表

求log(n)次rank数组，每次的rank数组都可以通过上次的rank数组得出

最后的rank就是答案，过程看起来简单，但代码实现挺难得！

Rank[]：rank[i]=p表示下标为i开头的字符串排第p个

Sa[]：Sa[i]=p表示排第i个的字符串是下标为p开头的字符串

cot[]：辅助数组，cot[i]=p表示当前rank数组中值小于等于i的数的个数

例如rank[]：{1,2,4,1,1,1,2,3}，那么cot[1]==4，cor[2]==6，cot[3]==7，cot[4]==8

（或初始字符串中ASCLL码<=i的字符个数）

Temp[]：Temp[i]=p意思就是第二个数（图里y数组）中第i的小的是下标p上的数

参考博文：

http://www.cnblogs.com/shanchuan04/p/5324009.html

http://blog.csdn.net/yxuanwkeith/article/details/50636898

http://blog.csdn.net/qq_34731703/article/details/52934284

代码中有注释

1031: [JSOI2007]字符加密Cipher

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 7202 Solved: 3120
[Submit][Status][Discuss]

Description

喜欢钻研问题的JS同学，最近又迷上了对加密方法的思考。一天，他突然想出了一种他认为是终极的加密办法：把需要加密的信息排成一圈，显然，它们有很多种不同的读法。例如下图，可以读作：

JSOI07 SOI07J OI07JS I07JSO 07JSOI 7JSOI0把它们按照字符串的大小排序：07JSOI 7JSOI0 I07JSO JSOI07

OI07JS SOI07J读出最后一列字符：I0O7SJ，就是加密后的字符串（其实这个加密手段实在很容易破解，鉴于这是突然想出来的，那就^^）。但是，如果想加密的字符串实在太长，你能写一个程序完成这个任务吗？

Input

　　输入文件包含一行，欲加密的字符串。注意字符串的内容不一定是字母、数字，也可以是符号等。

Output

　　输出一行，为加密后的字符串。

Sample Input

JSOI07

Sample Output

I0O7SJ

题解：裸的后缀数组，因为是个环不是个后缀所以要将字符串复制一遍接在后面

然后就是模板

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
char str[200100];
int a[200100], Rank[200100], Temp[200100], cot[200100], Sa[200100];
int main(void)
{int n, i, m, p, k;while(scanf("%s", str+1)!=EOF){n = strlen(str+1);for(i=1;i<=n;i++){a[i] = a[i+n] = str[i];str[i+n] = str[i];}n *= 2, m = 128;memset(cot, 0, sizeof(cot));for(i=1;i<=n;i++){Rank[i] = a[i];cot[a[i]]++;}for(i=1;i<=m;i++)cot[i] += cot[i-1];           //  cot[i]表示ASCLL码小于等于i的字符个数for(i=n;i>=1;i--)Sa[cot[a[i]]--] = i;/*puts(str+1);for(i=1;i<=n;i++)printf("%d ", Sa[i]);                       //  第一次排序，输入ABAAB，输出结果为：printf("\n");                                    //  ABAABABAABfor(i=1;i<=n;i++)                              //  1 3 4 6 8 9 2 5 7 10printf("%d ", Rank[i]);                       //  97 98 97 97 98 97 98 97 97 98printf("\n");*/for(k=1;k<=n;k*=2){p = 0;for(i=n-k+1;i<=n;i++)         //下标在[n+1, n+k]范围内的字符是空字符，所以字典序肯定是最小的Temp[++p] = i;for(i=1;i<=n;i++)         //前k个字符因为不参与排序，后面所有字符自动向前面移动k位（下标-k）{if(Sa[i]>k)Temp[++p] = Sa[i]-k;}/*printf("k=%d\nTemp: ", k);for(i=1;i<=n;i++)                     //  当k=1时，输出结果为：printf("%d ", Temp[i]);                //  10 2 3 5 7 8 1 4 6 9printf("\n");*/memset(cot, 0, sizeof(cot));for(i=1;i<=n;i++)                   /*这个时候，可以把rank当成新的a数组，之前的数组已经毫无意义*/cot[Rank[i]]++;for(i=1;i<=m;i++)cot[i] += cot[i-1];for(i=n;i>=1;i--)Sa[cot[Rank[Temp[i]]]--] = Temp[i];/* Temp[i]=p意思就是第二个数（图里y数组）中第i的大的是下标p上的数*/swap(Rank, Temp);p = 1;Rank[Sa[1]] = 1;         //求出当前的rank数组，在这里之前的rank数组就已经毫无意义了for(i=2;i<=n;i++){if(Temp[Sa[i-1]]==Temp[Sa[i]] && Temp[Sa[i-1]+k]==Temp[Sa[i]+k])Rank[Sa[i]] = p;elseRank[Sa[i]] = ++p;}m = p;/*printf("Rank: ");for(i=1;i<=n;i++)                 //  当k=1时，输出结果为：printf("%d ", Rank[i]);            //  2 4 1 2 4 2 4 1 2 3printf("\n");*/}n /= 2;for(i=1;i<=2*n;i++){if(Sa[i]>n)continue;printf("%c", str[Sa[i]+n-1]);}printf("\n");}
}