[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