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

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。

分析：自从开学以后就没有打过后缀数组这种鬼东西，连模板都忘得差不多了。以后还是要多复习才行。

这题其实就是把原串接到自己后面，把sa求出来，那么就可以得到每一个位置作为开头得到的字符串的排名了。

代码：

#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#define N 200005using namespace std;int len,b[N],c[N],d[N],sa[N],rank[N],height[N];char s[N],ans[N];void getsa(int n,int m){for (int i=0;i<=m;i++) b[i]=0;for (int i=1;i<=n;i++) b[s[i]]++;for (int i=1;i<=m;i++) b[i]+=b[i-1];for (int i=n;i>=1;i--) c[b[s[i]]--]=i;int t=0;for (int i=1;i<=n;i++){if (s[c[i]]!=s[c[i-1]]) t++;rank[c[i]]=t;}int j=1;while (j<=n){for (int i=0;i<=n;i++) b[i]=0;for (int i=1;i<=n;i++) b[rank[i+j]]++;for (int i=1;i<=n;i++) b[i]+=b[i-1];for (int i=n;i>=1;i--) c[b[rank[i+j]]--]=i;for (int i=0;i<=n;i++) b[i]=0;for (int i=1;i<=n;i++) b[rank[i]]++;for (int i=1;i<=n;i++) b[i]+=b[i-1];for (int i=n;i>=1;i--) d[b[rank[c[i]]]--]=c[i];int t=0;for (int i=1;i<=n;i++){if (rank[d[i]]!=rank[d[i-1]]||rank[d[i]]==rank[d[i-1]]&&rank[d[i]+j]!=rank[d[i-1]+j]) t++;c[d[i]]=t;}for (int i=1;i<=n;i++) rank[i]=c[i];if (t==n) break;j*=2;}for (int i=1;i<=n;i++)sa[rank[i]]=i;}void getheight(int n){int k=0;for (int i=1;i<=n;i++){if (k) k--;int j=sa[rank[i]-1];while (i+k<=n&&j+k<=n&&s[i+k]==s[j+k]) k++;height[rank[i]]=k;}}int main(){scanf("%s",s);int len=strlen(s),m=0;for (int i=len;i>=1;i--){m=max(m,(int)s[i]);s[i+len]=s[i]=s[i-1];}getsa(len*2,m);getheight(len*2);int j=0;for (int i=1;i<=len*2;i++)if (sa[i]<=len) ans[++j]=s[sa[i]+len-1];for (int i=1;i<=len;i++)printf("%c",ans[i]);return 0;}