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

最新推荐文章于 2018-12-22 14:23:22 发布

wzq_QwQ

最新推荐文章于 2018-12-22 14:23:22 发布

阅读量1.6k

点赞数 1

分类专栏： Suffix Array BZOJ刷题录文章标签：数组加密 namespace 代码

本文链接：https://blog.csdn.net/wzq_QwQ/article/details/48181693

版权

BZOJ刷题录同时被 2 个专栏收录

225 篇文章 13 订阅

订阅专栏

Suffix Array

6 篇文章 0 订阅

订阅专栏

题意:链接

方法:后缀数组

解析:

对于本题来说，我们只需要将与原串一模一样的串接到原串的后面，然后再求一下SA，之后输出合法的前n个即可。

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 100010
using namespace std;
char s[N<<1];
int sumfix[310];
int rnk[N<<1];
int sa[N<<1];
int use[N<<1];
int que[N<<1];
int n;
int cmp(int *z,int a,int b,int l)
{
    return z[a]==z[b]&&z[a+l]==z[b+l];
}
void get_sa(int lim)
{
    int *x=use,*y=que,*t;
    for(int i=0;i<lim;i++)sumfix[i]=0;
    for(int i=0;i<n;i++)x[i]=s[i];
    for(int i=0;i<n;i++)sumfix[x[i]]++;
    for(int i=1;i<lim;i++)sumfix[i]+=sumfix[i-1];
    for(int i=n-1;i>=0;i--)sa[--sumfix[x[i]]]=i;
    int nowlen=1,head=0;
    for(nowlen=1,head=1;head<n;nowlen<<=1,lim=head)
    {
        head=0;
        for(int i=n-nowlen;i<n;i++)y[head++]=i;
        for(int i=0;i<n;i++)if(sa[i]>=nowlen)y[head++]=sa[i]-nowlen;
        for(int i=0;i<lim;i++)sumfix[i]=0;
        for(int i=0;i<n;i++)sumfix[x[y[i]]]++;
        for(int i=1;i<lim;i++)sumfix[i]+=sumfix[i-1];
        for(int i=n-1;i>=0;i--)sa[--sumfix[x[y[i]]]]=y[i];
        int i;
        for(t=x,x=y,y=t,head=1,x[sa[0]]=0,i=1;i<n;i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],nowlen)?head-1:head++;
    }
} 
int main()
{
    scanf("%s",s);
    n=strlen(s);
    for(int i=n;i<n*2;i++)s[i]=s[i-n];n=2*n+1;
    //s[n-1]=0;
    get_sa(300);
    for(int i=0;i<n;i++)
        if(sa[i]<(n>>1))printf("%c",s[sa[i]+(n>>1)-1]);
}