1031: [JSOI2007]字符加密Cipher

1031: [JSOI2007]字符加密Cipher

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

分析

把字符串复制一份到原字符串的后面,然后后缀排序,输出。

code

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 
 5 using namespace std;
 6 
 7 const int N = 200100;
 8 
 9 char s[N];
10 int sa[N],t1[N],t2[N],c[N];
11 int n,m = 130;
12 
13 void get_sa() {
14     int i,p,*x = t1,*y = t2;
15     for (i=0; i<m; ++i) c[i] = 0;
16     for (i=0; i<n; ++i) 
17     x[i] = s[i],c[x[i]]++;
18     for (i=1; i<m; ++i) c[i] += c[i-1];
19     for (i=n-1; i>=0; --i) sa[--c[x[i]]] = i;
20     for (int k=1; k<=n; k<<=1) {
21         p = 0;
22         for (i=n-k; i<n; ++i) y[p++] = i;
23         for (i=0; i<n; ++i) if(sa[i]>=k) y[p++]=sa[i]-k;
24         for (i=0; i<m; ++i) c[i] = 0;
25         for (i=0; i<n; ++i) c[ x[y[i]] ]++;
26         for (i=1; i<m; ++i) c[i] += c[i-1];
27         for (i=n-1; i>=0; --i) sa[--c[ x[y[i]] ]] = y[i];
28         swap(x,y);
29         p = 1;
30         x[sa[0]] = 0;
31         for (i=1; i<n; ++i) 
32             x[sa[i]] = (y[sa[i-1]]==y[sa[i]] && sa[i-1]+k<n && sa[i]+k<n &&    
33             y[sa[i-1]+k]==y[sa[i]+k]) ? p-1 : p++;
34         if (p >= n) break;
35         m = p;
36     }
37 }
38 
39 int main() {
40     scanf("%s",s);
41     n = strlen(s);
42     for (int i=0; i<n; ++i) {
43         s[i+n] = s[i];    
44     }
45     n = n * 2;
46     get_sa();
47     for (int i=0; i<n; ++i) {
48         if (sa[i] < (n>>1)) printf("%c",s[sa[i]+(n>>1)-1]);    
49     }
50     return 0;    
51 }
原文地址:https://www.cnblogs.com/mjtcn/p/8456191.html