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

1031: [JSOI2007]字符加密Cipher

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 5485 Solved: 2266
[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
HINT
对于100%的数据字符串的长度不超过100000。
Source

【题解】后缀数组倍增法模板题

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define INF 300
using namespace std;
char s[500010];
int a[500010],sa[500010],c[500010],x[500010],y[500010];
int len;
inline void hzf()
{
    int n=INF,i,k;
    for(i=0;i<n;++i) c[i]=0;
    for(i=0;i<len;++i) c[x[i]=s[i]]++;
    for(i=1;i<n;++i) c[i]+=c[i-1];
    for(i=len-1;i>=0;--i) sa[--c[x[i]]]=i;
    for(k=1;k<=len;k<<=1)
     {
        int p=0;
        for(i=len-k;i<len;++i) y[p++]=i;
        for(i=0;i<len;++i)
         if(sa[i]>=k)  y[p++]=sa[i]-k;

        for(i=0;i<n;++i) c[i]=0;
        for(i=0;i<len;++i) c[x[y[i]]]++;
        for(i=1;i<n;++i) c[i]+=c[i-1];
        for(i=len-1;i>=0;--i) sa[--c[x[y[i]]]]=y[i];

        swap(x,y);
        p=1; x[sa[0]]=0;
        for(i=1;i<len;++i)
         x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?p-1:p++;
        if(p>=len) break;
        n=p;
     }
}
int main()
{
    int i;
    gets(s);
    len=strlen(s);
    for(i=0;i<len;++i) s[i+len]=s[i];
    len*=2;
    hzf();
    for(i=0;i<len;++i)
     if(sa[i]<len/2)
      putchar(s[sa[i]+len/2-1]);
    printf("\n");
    return 0;
}