首页 > 代码库 > BZOJ 1031: [JSOI2007]字符加密Cipher 后缀数组

BZOJ 1031: [JSOI2007]字符加密Cipher 后缀数组

1031: [JSOI2007]字符加密Cipher

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 6014  Solved: 2503
[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
 
题解
    题目要求就是要形成环的串了,那么把这个字符串接起来跑SA就好了
#include <iostream>#include <cstdio>#include <cmath>#include <cstring>#include <algorithm>#include<vector>using namespace std;#pragma comment(linker, "/STACK:102400000,102400000")#define ls i<<1#define rs ls | 1#define mid ((ll+rr)>>1)#define pii pair<int,int>#define MP make_pairtypedef long long LL;const long long INF = 1e18+1LL;const double Pi = acos(-1.0);const int N = 3e5+10, M = 2e5+20, mod = 1e9+7, inf = 2e9;///heght[i] 表示 Suffix(sa[i-1])和Suffix(sa[i]) 的最长公共前缀:///rank[i] 表示 开头为i的后缀的等级:///sa[i] 表示 排名为i的后缀 的开头位置:int *rank,r[N],sa[N],height[N],wa[N],wb[N],wm[N];bool cmp(int *r,int a,int b,int l) {    return r[a] == r[b] && r[a+l] == r[b+l];}void SA(int *r,int *sa,int n,int m) {    int *x=wa,*y=wb,*t;    for(int i=0;i<m;++i)wm[i]=0;    for(int i=0;i<n;++i)wm[x[i]=r[i]]++;    for(int i=1;i<m;++i)wm[i]+=wm[i-1];    for(int i=n-1;i>=0;--i)sa[--wm[x[i]]]=i;    for(int i=0,j=1,p=0;p<n;j=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<m;++i)wm[i]=0;        for(i=0;i<n;++i)wm[x[y[i]]]++;        for(i=1;i<m;++i)wm[i]+=wm[i-1];        for(i=n-1;i>=0;--i)sa[--wm[x[y[i]]]]=y[i];        for(t=x,x=y,y=t,i=p=1,x[sa[0]]=0;i<n;++i) {            x[sa[i]]=cmp(y,sa[i],sa[i-1],j)?p-1:p++;        }    }    rank=x;}void Height(int *r,int *sa,int n) {    for(int i=0,j=0,k=0;i<n;height[rank[i++]]=k)    for(k?--k:0,j=sa[rank[i]-1];r[i+k] == r[j+k];++k);}char s[N],ans[N];int cnt = 0;int main() {        scanf("%s",s);        int n = strlen(s);        int mtp = n;        for(int i = 0; i < n; ++i) r[i] = s[i] + 1;        for(int i = n; i < n + n; ++i) r[i] = s[i-n] + 1;        r[n+n] = 0;        n = 2 * n;        SA(r,sa,n+1,1000);        Height(r,sa,n);        for(int i = 1; i <= n; ++i) {            if(sa[i] < mtp) {                ans[++cnt] = s[(sa[i]-1+mtp)%mtp];            }        }        for(int i = 1; i <= cnt; ++i) cout<<ans[i];        cout<<endl;        return 0;}

 

BZOJ 1031: [JSOI2007]字符加密Cipher 后缀数组