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("");
}