字符加密Cipher

1031: [JSOI2007]字符加密Cipher

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 8920  Solved: 4034
[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。

#pragma GCC optimize(2)
#include <bits/stdc++.h>

typedef long long ll;
using namespace std;
const int maxn = 2e5 + 100000;
char str[maxn];
int sa[maxn],rk[maxn],ton[maxn],n,mx,tp[maxn];

inline void rsort(){
    for(register int i=0;i<=mx;++i)ton[i]=0;
    for(register int i=1;i<=n;++i)++ton[rk[i]];
    for(register int i=1;i<=mx;++i)ton[i]=ton[i]+ton[i-1];
    for(register int i=n;i>=1;--i)sa[ton[rk[tp[i]]]--]=tp[i];
}

inline void solve(){
    for(register int i=1;i<=n;++i)rk[i]=str[i],tp[i]=i;
    rsort();
    for(register int w=1,p=0;w<=n&&p<n;w*=2,mx=p){
        p=0;
        for(register int i=n-w+1;i<=n;++i)tp[++p]=i;
        for(register int i=1;i<=n;++i){
            if(sa[i]>w){
                tp[++p]=sa[i]-w;
            }
        }
        rsort();
        swap(rk,tp);
        rk[sa[1]]=p=1;
        for(register int i=2;i<=n;++i){
            rk[sa[i]]=((tp[sa[i]]==tp[sa[i-1]])&&(tp[sa[i]+w]==tp[sa[i-1]+w]))?p:++p;
        }
        if(p>=n)break;
    }
    for(register int i=1;i<=n;++i){
        rk[sa[i]]=i;
    }
}

int main() {
#ifndef ONLINE_JUDGE
    freopen("1.txt", "r", stdin);
#endif
    scanf("%s",str+1);
//    for(register int i=1;i<=100000;++i){
//        str[i]=rand()%10+'#';
//        printf("%c",str[i]);
//    }
//    printf("
");
    int init=strlen(str+1);
    //strcat(str+1,str+1);
    for(register int i=init+1;i<=2*init;++i)str[i]=str[i-init];
    n=strlen(str+1);
    str[n+1]=0;
    ++n;
    for(register int i=1;i<=n;++i){
        if(str[i]>mx)mx=str[i];
    }
//    mx=600;
//    printf("%s",str+1);
    solve();
//    for(register int i=1;i<=n;++i){
//        printf("debug sa[%d] = %d rk[%d] = %d
",i,sa[i],i,rk[i]);
//    }
    for(register int i=1;i<=n;++i){
        if(sa[i]<=init){
            printf("%c",str[sa[i]+init-1]);
        }
    }
    return 0;
}