【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]='\n'; 
    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+ - - 不过感谢灰哥

——It's a lonely path. Don't make it any lonelier than it has to be.
posted @ 2016-04-22 23:42  DaD3zZ  阅读(205)  评论(0编辑  收藏  举报