【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=0; i<m; i++) ws[i]=0;   
    for(int i=0; i<n; i++) ws[x[i]=r[i]]++; 
    for(int i=1; i<m; i++) ws[i]+=ws[i-1];   
    for(int i=n-1; i>=0; i--) sa[--ws[x[i]]]=i; 
    p=1;for(int j=1;p<n;j*=2,m=p)   
    {   
        p=0; for(int i=n-j;i<n;i++) y[p++]=i; 
        for(int i=0; i<n; i++) if(sa[i]>=j) y[p++]=sa[i]-j;   
        for(int i=0; i<n; i++) wv[i]=x[y[i]];   
        for(int i=0; i<m; i++) ws[i]=0;   
        for(int i=0; i<n; i++) ws[wv[i]]++;   
        for(int i=1; i<m; i++) ws[i]+=ws[i-1];   
        for(int i=n-1; i>=0; i--) sa[--ws[wv[i]]]=y[i]; 
        t=x,x=y,y=t,p=1,x[sa[0]]=0; 
        for(int i=1; i<n;i++)   
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;   
    }    
}   
int main() 
{ 
    scanf("%s",S); 
    int len=strlen(S); 
    for (int i=len; i<len*2; i++) S[i]=S[i-len]; 
    len=len*2+1; S[len]='
'; 
    DA(S,SA,len,300); 
    for (int i=0; i<len; i++) 
        if (SA[i]<(len>>1)) 
            putchar(S[SA[i]+(len>>1)-1]);  
    return 0; 
}

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