BZOJ1031: [JSOI2007]字符加密Cipher

浏览数：28 / 时间：2015年06月09日

1031: [JSOI2007]字符加密Cipher

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 3146 Solved: 1226
[Submit][Status]

Description

喜欢钻研问题的JS 同学，最近又迷上了对加密方法的思考。一天，他突然想出了一种他认为是终极的加密办法：把需要加密的信息排成一圈，显然，它们有很多种不同的读法。例如下图，可以读作： JSOI07 SOI07J OI07JS I07JSO 07JSOI 7JSOI0 把它们按照字符串的大小排序： 07JSOI 7JSOI0 I07JSO JSOI07 OI07JS SOI07J 读出最后一列字符：I0O7SJ，就是加密后的字符串（其实这个加密手段实在很容易破解，鉴于这是突然想出来的，那就^^）。但是，如果想加密的字符串实在太长，你能写一个程序完成这个任务吗？

Input

输入文件包含一行，欲加密的字符串。注意字符串的内容不一定是字母、数字，也可以是符号等。

Output

输出一行，为加密后的字符串。

Sample Input

JSOI07

Sample Output

I0O7SJ

数据规模
对于40%的数据字符串的长度不超过10000。
对于100%的数据字符串的长度不超过100000。

HINT

Source

题解：

以前不知道怎么水过的。。。

把原串复制一遍接到后面，做后缀数组，然后扫一遍就可

没在结尾+辅助字符导致昨天晚上没A

代码：

 1 #include<cstdio>
 2 #include<cstdlib>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #include<iostream>
 7 #include<vector>
 8 #include<map>
 9 #include<set>
10 #include<queue>
11 #include<string>
12 #define inf 1000000000
13 #define maxn 250000+10000
14 #define maxm 500+100
15 #define eps 1e-10
16 #define ll long long
17 #define pa pair<int,int>
18 using namespace std;
19 string s,ss;
20 int n,a[maxn],sa[maxn],wa[maxn],wb[maxn],wr[maxn],wv[maxn];
21 inline bool cmp(int *r,int a,int b,int l)
22 {
23     return r[a]==r[b]&&r[a+l]==r[b+l];
24 }
25 void da(int *r,int *sa,int n,int m)
26 {
27     int i,j,p,*x=wa,*y=wb,*t;
28         for(i=0;i<m;i++)wr[i]=0;
29     for(i=0;i<n;i++)wr[x[i]=r[i]]++;
30     for(i=1;i<m;i++)wr[i]+=wr[i-1];
31     for(i=n-1;i>=0;i--)sa[--wr[x[i]]]=i;    
32     for(j=1,p=1;p<n;j<<=1,m=p)
33      { 
34          for(p=0,i=n-j;i<n;i++)y[p++]=i;
35          for(i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
36          for(i=0;i<m;i++)wr[i]=0;
37          for(i=0;i<n;i++)wr[x[y[i]]]++;
38          for(i=1;i<m;i++)wr[i]+=wr[i-1];
39          for(i=n-1;i>=0;i--)sa[--wr[x[y[i]]]]=y[i];
40          for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
41           x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
42      }
43 }
44 int main()
45 {
46     freopen("input.txt","r",stdin);
47     freopen("output.txt","w",stdout);
48     cin>>s;n=s.length();
49     for(int i=0;i<n;i++)a[i]=int(s[i]);
50     for(int i=0;i<n;i++)a[n+i]=int(s[i]);
51     a[2*n]=0;
52     da(a,sa,2*n+1,129);
53     ss="";
54     for(int i=0;i<=2*n;i++)
55      if(sa[i]<n)ss=ss+s[(sa[i]-1+n)%n];
56     cout<<ss<<endl; 
57     return 0;
58 }