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

1031: [JSOI2007]字符加密Cipher

Time Limit: 20 Sec  Memory Limit: 256 MB

题目连接

http://acm.hust.edu.cn/problem/show/1590

Description

 

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

 

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

Input

输入文件包含一行，欲加密的字符串。注意字符串的内容不一定是字母、数字，也可以是符号等。

Output

输出一行，为加密后的字符串。

 

Sample Input

JSOI07

Sample Output

I0O7SJ

HINT

 

题意

 

题解：

后缀数组裸题，首先把字符串延长一半，然后跑一发后缀数组，他自己就可以排序

然后判断一下这个字符串的长度，然后就输出就好了

后缀数组代码见下，已经加了注释了

代码：

 

//qscqesze

#include <cstdio>

#include <cmath>

#include <cstring>

#include <ctime>

#include <iostream>

#include <algorithm>

#include <set>

#include <vector>

#include <sstream>

#include <queue>

#include <typeinfo>

#include <fstream>

#include <map>

#include <stack>

typedef long long ll;

using namespace std;

//freopen("D.in","r",stdin);

//freopen("D.out","w",stdout);

#define sspeed ios_base::sync_with_stdio(0);cin.tie(0)

#define maxn 1000001

#define mod 10007

#define eps 1e-9

int Num;

char CH[20];

//const int inf=0x7fffffff;   //нчоч╢С

const int inf=0x3f3f3f3f;

/*



inline void P(int x)

{

    Num=0;if(!x){putchar('0');puts("");return;}

    while(x>0)CH[++Num]=x%10,x/=10;

    while(Num)putchar(CH[Num--]+48);

    puts("");

}

*/

inline ll read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}

    return x*f;

}

inline void P(int x)

{

    Num=0;if(!x){putchar('0');puts("");return;}

    while(x>0)CH[++Num]=x%10,x/=10;

    while(Num)putchar(CH[Num--]+48);

    puts("");

}

//**************************************************************************************



char s[2*maxn];

int n;

int sa[maxn],ra[maxn],tmp[maxn],height[maxn];

int top[maxn+10];

void build_sa()

{

    int m = 128; //最大字符值

    int *x = ra, *y = tmp; // x 第一关键字 y 备份

    // 对单个字符进行基数排序

    for(int i = 0; i < m; i++) top[i] = 0;

    for(int i = 0; i < n; i++) top[x[i] = s[i]]++;

    for(int i = 1; i < m; i++) top[i] += top[i-1];

    for(int i = n-1; i >= 0; i--) sa[--top[x[i]]] = i;

    for(int k = 1; k <= n; k <<= 1) {

        int p = 0;

        //利用上一次的sa值算出第二关键字的排序

        for(int i = n - k; i < n; i++) y[p++] = i;

        for(int i = 0; i < n; i++) if(sa[i] >= k) y[p++] = sa[i] - k;

        //对第一关键字进行基数排序

        for(int i = 0; i < m; i++) top[i] = 0;

        for(int i = 0; i < n; i++) top[x[y[i]]]++;

        for(int i = 0; i < m; i++) top[i] += top[i-1];

        for(int i = n-1; i >= 0; i--) sa[--top[x[y[i]]]] = y[i];

        swap(x, y);

        p = 1; x[sa[0]] = 0;

        // 计算rank值

        for(int i = 1; i < n; i++)

            x[sa[i]] = (y[sa[i-1]] == y[sa[i]] && y[sa[i-1]+k] == y[sa[i]+k]) ? p-1 : p++;

        if(p >= n) break;

        m = p;

    }

}

void build_height(){

    //rep (i,0,n) rank[sa[i]] = i;

    int k = 0;

    for(int i=0;i<n;i++) ra[sa[i]] = i;

    for(int i=0;i<n;i++){

        if (k) k--;

        int j = sa[ra[i] - 1];

        while (s[i+k] == s[j+k]) k++;

        height[ra[i]] = k;

    }

}

int main()

{

    scanf("%s", s);

    int len = strlen(s);

    for(int i=len;i<2*len;i++) s[i] = s[i-len];

    n = 2*len;

    build_sa();

    for(int i=1;i<n;i++)

    {

        if (sa[i] < len)

        {

            putchar(s[sa[i]+len-1]);

        }

    }

    puts("");

}