JSOI2007——字符加密Cipher(后缀数组)

1031: [JSOI2007]字符加密Cipher

Time Limit: 10 Sec  Memory Limit: 162 MB

Submit: 1684  Solved: 640

[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
数据规模
对于40%的数据字符串的长度不超过10000。

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

解析：

        很水的后缀数组，练练代码。。。

代码：

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define N 100010
int n,m=256;
char s[N];
int ws[N],wv[N],wa[N],wb[N],sa[N];

void read()
{
    freopen("Cipher.in","r",stdin);
    freopen("Cipher.out","w",stdout);
}

int cmp(int *r,int a,int b,int l)
{
    return (r[a]==r[b]&&r[a+l]==r[b+l]);
}

void DA(char *s,int *sa,int n,int m)
{
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0;i<m;i++)ws[i]=0;
    for(i=0;i<n;i++)ws[x[i]=s[i]]++;
    for(i=1;i<m;i++)ws[i]+=ws[i-1];
    for(i=n-1;i>=0;i--)sa[--ws[x[i]]]=i;
    for(j=1,p=1;p<n;j*=2,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<n;i++)wv[i]=x[y[i]];
        for(i=0;i<m;i++)ws[i]=0;
        for(i=0;i<n;i++)ws[wv[i]]++;
        for(i=1;i<m;i++)ws[i]+=ws[i-1];
        for(i=n-1;i>=0;i--)sa[--ws[wv[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()
{
    read();
    scanf("%s",s);
    n=strlen(s);
    for(int i=0;i<n;i++)
    s[i+n]=s[i];
    s[2*n]=0;
    DA(s,sa,n*2+1,m);
    for(int i=0;i<=2*n;i++)
    if(sa[i]<n)printf("%c",s[sa[i]+n-1]);
}