【BZOJ】1031 [JSOI2007]字符加密Cipher

Description

喜欢钻研问题的JS同学，最近又迷上了对加密方法的思考。一天，他突然想出了一种他认为是终极的加密办法：把需要加密的信息排成一圈，显然，它们有很多种不同的读法。例如下图，可以读作：
【图不见啦~~~】
JSOI07 SOI07J OI07JS I07JSO 07JSOI 7JSOI0把它们按照字符串的大小排序：07JSOI 7JSOI0 I07JSO JSOI07 OI07JS SOI07J读出最后一列字符：I0O7SJ，就是加密后的字符串（其实这个加密手段实在很容易破解，鉴于这是突然想出来的，那就^^）。但是，如果想加密的字符串实在太长，你能写一个程序完成这个任务吗？

Input

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

Output

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

Solution

这题是一个字符的同构排序（实际上拆环后就是后缀排序）
后缀排序用了一个阉割版的后缀数组做了

#include<stdio.h>
#include<cstring>
#define N 200005
#define M 400005

char s[M];
int sum[N],sa[N],rank[M],trank[N],tsa[N];

int main()
{
    scanf("%s",s+1);
    int n=strlen(s+1);
    for (int i=1;i<=n;i++) s[i+n]=s[i];n<<=1;
    for (int i=1;i<=n;i++) sum[s[i]]++;
    for (int i=1;i<=256;i++) sum[i]+=sum[i-1];
    for (int i=1;i<=n;i++) sa[sum[s[i]]--]=i;
    for (int i=1;i<=n;i++) rank[sa[i]]=rank[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);
    for (int j=1;j<n;j<<=1)
    {
        memset(sum,0,sizeof(sum));
        for (int i=1;i<=n;i++) sum[rank[i+j]]++;
        for (int i=1;i<=n;i++) sum[i]+=sum[i-1];
        for (int i=n;i;i--) tsa[sum[rank[i+j]]--]=i;
        memset(sum,0,sizeof(sum));
        for (int i=1;i<=n;i++) sum[rank[i]]++;
        for (int i=1;i<=n;i++) sum[i]+=sum[i-1];
        for (int i=n;i;i--) sa[sum[rank[tsa[i]]]--]=tsa[i];
        for (int i=1;i<=n;i++) trank[sa[i]]=trank[sa[i-1]]+(rank[sa[i]]!=rank[sa[i-1]] || rank[sa[i]+j]!=rank[sa[i-1]+j]);
        for (int i=1;i<=n;i++) rank[i]=trank[i];
    }
    for (int i=1;i<=n;i++) if (sa[i]<=(n>>1)) putchar(s[sa[i]+(n>>1)-1]);
}

邪教……邪教……（只慢了不到一倍……）

#include<stdio.h>
#include<algorithm>
#include<cstring>
using namespace std;
#define N 200005

char s[N];
int n,rk[N<<1];
#define mp make_pair
#define fi first
#define se second
pair<pair<int,int>,int> p[N];

int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);
    for (int i=1;i<=n;i++) s[n+i]=s[i];
    n<<=1;
    for (int i=1;i<=n;i++) rk[i]=s[i];
    for (int j=1;j<n;j<<=1)
    {
        for (int i=1;i<=n;i++) p[i]=mp(mp(rk[i],rk[i+j]),i);
        sort(p+1,p+n+1);
        for (int i=1;i<=n;i++) rk[p[i].se]=rk[p[i-1].se]+(p[i].fi!=p[i-1].fi);
    }
    for (int i=1;i<=n;i++) if ((p[i].se<<1)<=n) putchar(s[p[i].se+(n>>1)-1]);
}

我已经斜了（第一份1.3s，这份不到1.7s）

#include<stdio.h>
#include<algorithm>
#include<cstring>
using namespace std;
#define N 200005

char s[N];
typedef unsigned int uint;
uint n;
unsigned long long p[N],rk[N<<1];

int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);
    for (uint i=1;i<=n;i++) s[n+i]=s[i];
    n<<=1;
    for (uint i=1;i<=n;i++) rk[i]=s[i];
    for (uint j=1;j<n;j<<=1)
    {
        for (uint i=1;i<=n;i++) p[i]=(rk[i]<<36)+(rk[i+j]<<18)+i;
        sort(p+1,p+n+1);
        for (uint i=1;i<=n;i++) rk[p[i]&0x3ffff]=rk[p[i-1]&0x3ffff]+((p[i]>>18)!=(p[i-1]>>18));
    }
    for (uint i=1;i<=n;i++) if (((p[i]&0x3ffff)<<1)<=n) putchar(s[(p[i]&0x3ffff)+(n>>1)-1]);
}