BZOJ1031: [JSOI2007]字符加密Cipher

本文详细介绍了一种基于字符串处理的后缀数组(SA)构建算法。通过将原字符串复制并进行特定处理,该算法能在O(n log n)的时间复杂度内完成后缀数组的构建。文中提供了一个完整的C++实现示例,并展示了如何通过枚举和条件判断来优化输出过程。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目链接

将原串复制一遍,后缀数组区间长度只枚举到n/2,即原串长,最后输出时枚举,若sa[i]小于n/2,输出s[sa[i]+n/2-1]。

【代码】

#include <cstdio>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cmath>
#define N 200005
#define M 10000005
#define INF 1e9
using namespace std;
typedef unsigned int ll;
typedef pair<ll,ll> pa;

int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n,m,nn;
char s[N];
int sa[N],t[N],t2[N],c[N];

void Build_Sa()
{
    int *x=t,*y=t2;
    for(int i=0;i<m;i++) c[i]=0;
    for(int i=0;i<n;i++) c[x[i]=s[i]]++;
    for(int i=1;i<m;i++) c[i]+=c[i-1];
    for(int i=n-1;i>=0;i--) sa[--c[x[i]]]=i;
    for(int k=1;k<=nn;k<<=1)
    {
        int p=0;
        for(int i=n-k;i<n;i++) y[p++]=i;
        for(int i=0;i<n;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
        for(int i=0;i<m;i++) c[i]=0;
        for(int i=0;i<n;i++) c[x[y[i]]]++;
        for(int i=1;i<m;i++) c[i]+=c[i-1];
        for(int i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
        swap(x,y);p=1;x[sa[0]]=0;
        for(int i=1;i<n;i++)
            x[sa[i]]=y[sa[i-1]+k]==y[sa[i]+k]&&y[sa[i-1]]==y[sa[i]]?p-1:p++;
        if(p>=n) break;m=p; 
    }
} 

int main()
{
    scanf("%s",s);nn=strlen(s);
    for(int i=0;i<nn-1;i++) s[i+nn]=s[i];
    n=nn*2-1;m=256;
    Build_Sa();
    for(int i=0;i<n;i++) 
        if(sa[i]<nn) printf("%c",s[sa[i]+nn-1]);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值