1031. [JSOI2007]字符加密【后缀数组】

Description

Input

Output

Sample Input

Sample Output

HINT

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

相同两个字符串连接然后裸后缀数组

#include<iostream>
#include<cstring>
#include<cstdio>
#define MAXN (200000+10)
using namespace std;
int wa[MAXN],wb[MAXN],wt[MAXN];
int SA[MAXN],Rank[MAXN],Height[MAXN];
int n,m;
char r[MAXN];

inline bool cmp(int *y,int a,int b,int k)
{
    int arank1=y[a];
    int brank1=y[b];
    int arank2=a+k>=n?-1:y[a+k];
    int brank2=b+k>=n?-1:y[b+k];
    return arank1==brank1 && arank2==brank2;
}

inline void Build_SA()
{
    int *x=wa,*y=wb;
    for (int i=0;i<m;++i) wt[i]=0;
    for (int i=0;i<n;++i) wt[x[i]=r[i]]++;
    for (int i=1;i<m;++i) wt[i]+=wt[i-1];
    for (int i=n-1;i>=0;--i) SA[--wt[x[i]]]=i;
    
    for (int j=1;j<=n;j<<=1)
    {
        int p=0;
        for (int i=n-j;i<n;++i) y[p++]=i;
        for (int i=0;i<n;++i) if (SA[i]>=j) y[p++]=SA[i]-j;
        
        for (int i=0;i<m;++i) wt[i]=0;
        for (int i=0;i<n;++i) wt[x[y[i]]]++;
        
        for (int i=1;i<m;++i) wt[i]+=wt[i-1];
        for (int i=n-1;i>=0;--i) SA[--wt[x[y[i]]]]=y[i];
        swap(x,y);m=1;
        x[SA[0]]=0;
        for (int i=1;i<n;++i)
            x[SA[i]]=cmp(y,SA[i-1],SA[i],j)?m-1:m++;
        if (m>=n) break;
    }
}

int main()
{
    scanf("%s",&r);
    n=strlen(r);
    for (int i=n;i<=n*2-1;++i) r[i]=r[i-n];
    n*=2;m=130;
    Build_SA();
    for (int i=0;i<n;i++)
        if (SA[i]<n/2)
            printf("%c",r[SA[i]+n/2-1]);
}