bzoj1031: [JSOI2007]字符加密Cipher（后缀数组）

1031: [JSOI2007]字符加密Cipher

题目：传送门

 

正解：

　　 后缀数组裸题。。。

　　 刚开始看的时候可能觉得一个个排序很麻烦，其实在后面接多一个相同的字符串之后再跑一边sa就ok

　　 因为对于每一个组成顺序不同的字符串来说，其实长度都是一样的，这就保证了他们在这个长度的范围内就可以比较出大小。

　　 所以我们所做的预处理是没有影响的。

　　 输出的时候再稍微判断一下就AC了...有点水

 

手撸模版：

　　 

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define N 1100000
using namespace std;
int Rank[N],Rsort[N],sa[N],y[N],wr[N];
char s[N];
int a[N];
int len,n;
bool cmp(int k1,int k2,int ln){return wr[k1]==wr[k2] && wr[k1+ln]==wr[k2+ln];}
void get_sa(int n,int m)
{
    memcpy(Rank,a,sizeof(Rank));
    
    memset(Rsort,0,sizeof(Rsort));
    for(int i=1;i<=n;i++)Rsort[Rank[i]]++;
    for(int i=1;i<=m;i++)Rsort[i]+=Rsort[i-1];
    for(int i=n;i>=1;i--)sa[Rsort[Rank[i]]--]=i;
    
    int ln=1,k,p=0;
    while(p<n)
    {
        k=0;
        for(int i=n-ln+1;i<=n;i++)y[++k]=i;
        for(int i=1;i<=n;i++)if(sa[i]>ln)y[++k]=sa[i]-ln;
        
        for(int i=1;i<=n;i++)wr[i]=Rank[y[i]];
        
        memset(Rsort,0,sizeof(Rsort));
        for(int i=1;i<=n;i++)Rsort[wr[i]]++;
        for(int i=1;i<=m;i++)Rsort[i]+=Rsort[i-1];
        for(int i=n;i>=1;i--)sa[Rsort[wr[i]]--]=y[i];
        
        for(int i=1;i<=n;i++)wr[i]=Rank[i];
        p=1;Rank[sa[1]]=1;
        for(int i=2;i<=n;i++)
        {
            if(!cmp(sa[i],sa[i-1],ln))p++;
            Rank[sa[i]]=p;
        }
        m=p;ln*=2;
    }
}
int main()
{
    scanf("%s",s+1);len=strlen(s+1);
    n=len*2;int j=0;
    for(int i=len+1;i<=n;i++)s[i]=s[++j];
    for(int i=1;i<=n;i++)a[i]=s[i];
    get_sa(n,1000);
    for(int i=1;i<=n;i++)
    {
        if(sa[i]>len)continue;
        else
        {
            printf("%c",a[sa[i]+len-1]);
        }
    }
    printf("\n");
    return 0;
}