BZOJ1031: [JSOI2007]字符加密Cipher

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

代码：

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 1000000000
#define maxn 250000+10000
#define maxm 500+100
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
using namespace std;
string s,ss;
int n,a[maxn],sa[maxn],wa[maxn],wb[maxn],wr[maxn],wv[maxn];
inline bool cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void da(int *r,int *sa,int n,int m)
{
    int i,j,p,*x=wa,*y=wb,*t;
        for(i=0;i<m;i++)wr[i]=0;
    for(i=0;i<n;i++)wr[x[i]=r[i]]++;
    for(i=1;i<m;i++)wr[i]+=wr[i-1];
    for(i=n-1;i>=0;i--)sa[--wr[x[i]]]=i;    
    for(j=1,p=1;p<n;j<<=1,m=p)
     { 
         for(p=0,i=n-j;i<n;i++)y[p++]=i;
         for(i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
         for(i=0;i<m;i++)wr[i]=0;
         for(i=0;i<n;i++)wr[x[y[i]]]++;
         for(i=1;i<m;i++)wr[i]+=wr[i-1];
         for(i=n-1;i>=0;i--)sa[--wr[x[y[i]]]]=y[i];
         for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
          x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
     }
}
int main()
{
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);
    cin>>s;n=s.length();
    for(int i=0;i<n;i++)a[i]=int(s[i]);
    for(int i=0;i<n;i++)a[n+i]=int(s[i]);
    a[2*n]=0;
    da(a,sa,2*n+1,129);
    ss="";
    for(int i=0;i<=2*n;i++)
     if(sa[i]<n)ss=ss+s[(sa[i]-1+n)%n];
    cout<<ss<<endl; 
    return 0;
}