bzoj1031 [JSOI2007]字符加密Cipher

Description

---

　喜欢钻研问题的JS同学，最近又迷上了对加密方法的思考。一天，他突然想出了一种他认为是终极的加密办法 
：把需要加密的信息排成一圈，显然，它们有很多种不同的读法。例如下图，可以读作： 

JSOI07 SOI07J OI07JS I07JSO 07JSOI 7JSOI0把它们按照字符串的大小排序：07JSOI 7JSOI0 I07JSO JSOI07
OI07JS SOI07J读出最后一列字符：I0O7SJ，就是加密后的字符串（其实这个加密手段实在很容易破解，鉴于这是 
突然想出来的，那就^^）。但是，如果想加密的字符串实在太长，你能写一个程序完成这个任务吗？

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

Solution

---

本蒟蒻的第一题后缀数组SA，纪念一下 
把原串复制一份加在后面，跑出来的所有后缀中长度>=len的按顺序排序，每个取最后一位就ok 
膜来的板子a啊b啊c啊d啊乱七八糟的有点晕，但是想不到更好的替代 
辣鸡csdn，卡屎了

Code

---

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<vector>
#define rep(i,st,ed)for(int i=st;i<=ed;++i)
#define drp(i,st,ed)for(int i=st;i>=ed;--i)
#define fill(x,t) memset(x,t,sizeof(x))
#define copy(x,t) memcpy(x,t,sizeof(x))

constint N=200005;

int rank[N],sa[N],h[N];
int b[N],c[N],d[N];

char s[N],ans[N];

void get_sa(int n,int m){
    rep(i,0,m) b[i]=0;
    rep(i,1,n) b[s[i]]++;
    rep(i,1,m) b[i]+=b[i-1];
    drp(i,n,1) c[b[s[i]]--]=i;
int t=0;
    rep(i,1,n){
if(s[c[i]]!=s[c[i-1]]) t++;
        rank[c[i]]=t;
}
int j=1;
while(j<=n){
        rep(i,0,n) b[i]=0;
        rep(i,1,n) b[rank[i+j]]++;
        rep(i,1,n) b[i]+=b[i-1];
        drp(i,n,1) c[b[rank[i+j]]--]=i;
        rep(i,0,n) b[i]=0;
        rep(i,1,n) b[rank[i]]++;
        rep(i,1,n) b[i]+=b[i-1];
        drp(i,n,1) d[b[rank[c[i]]]--]=c[i];
int t=0;
        rep(i,1,n){
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;
}
        rep(i,1,n) rank[i]=c[i];
if(t==n)break;
        j*=2;
}
    rep(i,1,n)
        sa[rank[i]]=i;
}

void get_height(int n){
int k=0;
    rep(i,1,n){
if(k) k--;
int j=sa[rank[i]-1];
while(i+k<=n&&j+k<=n&&s[i+k]==s[j+k]) k++;
        h[rank[i]]=k;
}
}

int main(void){
scanf("%s", s);int len=strlen(s),m=0;
    drp(i,len,1){
        s[i+len]=s[i]=s[i-1];
        m=std:: max(m,(int)s[i]);
}
    get_sa(len*2,m);
    get_height(len*2);
int j=0;
    rep(i,1,len*2)if(sa[i]<=len) ans[++j]=s[sa[i]+len-1];
    rep(i,1,len)printf("%c",ans[i]);
return0;
}

 本文转载自：http://www.javaxxz.com/thread-359193-1-1.html