【BZOJ-1031】字符加密Cipher 后缀数组

1031: [JSOI2007]字符加密Cipher

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit:
5504  Solved: 2277
[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

Solution

非常裸的后缀数组

既然是处理环，那就复制一下，裸上SA即可

Code

#include<cstdio>
#include<cstring>
using namespace std;
#define maxn 2000010
int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
char S[maxn]; int SA[maxn];
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void DA(char *r,int *sa,int n,int m)
{
    int p,*x=wa,*y=wb,*t;
    for(int i=; i<m; i++) ws[i]=;
    for(int i=; i<n; i++) ws[x[i]=r[i]]++;
    for(int i=; i<m; i++) ws[i]+=ws[i-];
    for(int i=n-; i>=; i--) sa[--ws[x[i]]]=i;
    p=;for(int j=;p<n;j*=,m=p)
    {
        p=; for(int i=n-j;i<n;i++) y[p++]=i;
        for(int i=; i<n; i++) if(sa[i]>=j) y[p++]=sa[i]-j;
        for(int i=; i<n; i++) wv[i]=x[y[i]];
        for(int i=; i<m; i++) ws[i]=;
        for(int i=; i<n; i++) ws[wv[i]]++;
        for(int i=; i<m; i++) ws[i]+=ws[i-];
        for(int i=n-; i>=; i--) sa[--ws[wv[i]]]=y[i];
        t=x,x=y,y=t,p=,x[sa[]]=;
        for(int i=; i<n;i++)
            x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
    }
}
int main()
{
    scanf("%s",S);
    int len=strlen(S);
    for (int i=len; i<len*; i++) S[i]=S[i-len];
    len=len*+; S[len]='\n';
    DA(S,SA,len,);
    for (int i=; i<len; i++)
        if (SA[i]<(len>>))
            putchar(S[SA[i]+(len>>)-]);
    return ;
}

看完论文，自己搞了搞，似乎出了点意外..找ZYF2000神犇帮改了一下，1A1000ms+，不想换姿势，于是自己结合别人的订正了一波，500ms+ - - 不过感谢灰哥