[JSOI2007]字符加密Cipher

bzoj 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<cstdio>
#include<cstring>
#include<algorithm>
#define N 200010  
using namespace std;
int n,len,p,q=1,k,a[N],rk[2][N],sa[2][N],v[N];
char s[N]; 
void mul(int *sa,int *rk,int *SA,int *RK)
{
    for(int i=1;i<=n;i++) v[rk[sa[i]]]=i;
    for(int i=n;i;i--) if(sa[i]>k) SA[v[rk[sa[i]-k]]--]=sa[i]-k;
    for(int i=n-k+1;i<=n;i++) SA[v[rk[i]]--]=i;
    for(int i=1;i<=n;i++) RK[SA[i]]=RK[SA[i-1]]+(rk[SA[i]]!=rk[SA[i-1]]||rk[SA[i]+k]!=rk[SA[i-1]+k]); 
} 
void presa()
{
    for(int i=1;i<=n;i++) v[a[i]]++;
    for(int i=1;i<=130;i++) v[i]+=v[i-1];
    for(int i=1;i<=n;i++) sa[p][v[a[i]]--]=i;
    for(int i=1;i<=n;i++) rk[p][sa[p][i]]=rk[p][sa[p][i-1]]+(a[sa[p][i]]!=a[sa[p][i-1]]);
    for(k=1;k<n;k<<=1,swap(p,q)) 
    {
        mul(sa[p],rk[p],sa[q],rk[q]);
        if(rk[q][sa[q][n]]==n) 
        {
            swap(p,q);
            return;
        }
    }
}
void solve()
{
    int h;
    for(int i=1;i<=n;i++) 
    {
        if(sa[p][i]>len) continue;
        h=sa[p][i]+len-1;
        putchar(s[h]);
    }
}
int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);
    len=n;
    for(int i=1;i<=len;i++) a[++n]=s[i],s[n]=s[i];
    for(int i=1;i<=n;i++) a[i]=s[i]-' '+1;
    presa();
    solve();
}

开始的错误做法:

直接做后缀数组,按顺序输出第h=(sa[p][i]-1+n-1)%n+1个字符

错误原因:例:zbb 

排序结果:bbz,bzb,zbb

正确的sa数组:2,3,1

错误方法得到的sa数组:3,2,1

posted @ 2017-03-21 10:51  TRTTG  阅读(324)  评论(0编辑  收藏  举报