bzoj 1031[JSOI2007]字符加密 - 后缀数组

　　　　　　　　　　　　　　　　

1031: [JSOI2007]字符加密Cipher

Time Limit: 10 Sec  Memory Limit: 162 MB

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。

 

应该是一道板子题吧

把字符串复制一遍搁到后边就行了

 

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#define LL long long
using namespace std;

const int MAXN = 3e5 + 10;
char str[MAXN];
int m;
int ra[MAXN], tp[MAXN];
int SA[MAXN];
int cur[MAXN];
int c[MAXN];
int len;
int p;
inline LL read()
{
    LL x = 0, w = 1; char ch = 0;
    while(ch < '0' || ch > '9') {
        if(ch == '-') {
            w = -1;
        }
        ch = getchar();
    }
    while(ch >= '0' && ch <= '9') {
        x = x * 10 + ch - '0';
        ch = getchar();
    }
    return x * w;
}

void sort(int m)
{
    for(int i = 0; i <= m; i++) {
        c[i] = 0;
    }
    for(int i = 1; i <= len; i++) {
        c[ra[tp[i]]]++;
    }
    for(int i = 1; i <= m; i++) {
        c[i] += c[i - 1];
    }
    for(int i = len; i >= 1; i--) {
        SA[c[ra[tp[i]]]--] = tp[i];
    }
}

void copy()
{
    for(int i = 1; i <= len; i++) {
        cur[i] = ra[i];
    }
}

void suffix()
{
    for(int i = 1; i <= len; i++) {
        ra[i] = str[i];
        tp[i] = i;
    } 
    sort(m = 128);
    for(int w = 1, p = 0; p < len; w += w, m = p) {
        p = 0;
        for(int j = len - w + 1; j <= len; j++) {
            tp[++p] = j;
        }
        for(int i = 1; i <= len; i++) {
            if(SA[i] > w) {
                tp[++p] = SA[i] - w;
            }
        }
        sort(m);
        copy();
        ra[SA[1]] = p = 1;
        for(int i = 2; i <= len; i++) {
            if(cur[SA[i]] == cur[SA[i - 1]] && cur[SA[i - 1] + w] == cur[SA[i] + w]) {
                ra[SA[i]] = p; 
            } else {
                ra[SA[i]] = ++p;
            }
        }
    }
}
int main()
{
    scanf("%s", str + 1);
    len = strlen(str + 1);
    for(int i = 1; i <= len; i++) {
        str[i + len] = str[i];
    }
    len += len;
    suffix();
    for(int i = 1; i <= len; i++) {
        if(SA[i] <= len / 2) 
        cout<<str[SA[i] + len / 2  - 1];
    }
    cout<<endl;
    return 0;
}