[BZOJ1031][JSOI2007]字符加密Cipher 后缀数组

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1031

求字符串排名的问题，直接上后缀数组。

由于这个字符串是环形的，我们把它倍增一下就好了。然后发现其实就是求sa数组裸题，后缀长度大于等于len的后缀就可以对应原来的一种字符串，直接输出就好了……

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m;
char s[200010];
int rk[200010],sa[200010],tmp[200010],c[200010];
bool inline cmp(int *r,int a,int b,int l){
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void Getsa(){
    m=300;
    for(int i=0;i<=m;i++) c[i]=0;
    for(int i=1;i<=n;i++) c[rk[i]=s[i]]++;
    for(int i=1;i<=m;i++) c[i]+=c[i-1];
    for(int i=n;i>=1;i--) sa[c[rk[i]]--]=i;
    int p;
    for(int i=1;i<=n;i<<=1){
        p=0;
        for(int j=n-i+1;j<=n;j++) tmp[++p]=j;
        for(int j=1;j<=n;j++) if(sa[j]>i) tmp[++p]=sa[j]-i;
        for(int j=0;j<=m;j++) c[j]=0;
        for(int j=1;j<=n;j++) c[rk[tmp[j]]]++;
        for(int j=1;j<=m;j++) c[j]+=c[j-1];
        for(int j=n;j>=1;j--) sa[c[rk[tmp[j]]]--]=tmp[j];
        swap(rk,tmp);
        rk[sa[1]]=p=1;
        for(int j=2;j<=n;j++) rk[sa[j]]=cmp(tmp,sa[j],sa[j-1],i)?p:++p;
        if(p>=n) break;
        m=p;
    }
}
int main(){
    scanf("%s",s+1);
    int len=strlen(s+1);
    for(int i=1;i<len;i++) s[len+i]=s[i];
    n=(len<<1)-1;
    Getsa();
    for(int i=1;i<=n;i++)
        if(sa[i]+len<=n+1)
            printf("%c",s[sa[i]+len-1]);
    return 0;
}