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

来源:互联网 发布:java hashmap原理 编辑:程序博客网 时间:2024/05/17 00:01

先把串复制一遍

然后对这个2n的串求后缀数组

然后输出sa值小于n的


开始我只是对n做后缀数组 wa不停 耽误了很长时间

如果只是对n求后缀数组的话 不能求出后面n的rank值 就不能倍增!!!


#include <cstring>#include <cstdio>#include <iostream>using namespace std;const int MAXN=200000+100;char s[MAXN];int n,sa[MAXN],rank[MAXN],tmp[MAXN],cnt[MAXN],sec[MAXN];void get_sa(){for(int i=1;i<=n;i++) cnt[s[i]]++;for(int i=1;i<=200;i++) cnt[i]+=cnt[i-1];for(int i=1;i<=n;i++) sa[cnt[s[i]]--]=i;int nn=rank[sa[1]]=1;for(int i=2;i<=n;i++){if(s[sa[i]]!=s[sa[i-1]]) nn++;rank[sa[i]]=nn;}for(int i=0;i<=200;i++) cnt[i]=0;for(int k=1;k<=n&&nn<n;k<<=1){int top=0;for(int i=n-k+1;i<=n;i++) sec[++top]=i;for(int i=1;i<=n;i++) if(sa[i]>k) sec[++top]=sa[i]-k;for(int i=0;i<=nn;i++) cnt[i]=0;for(int i=1;i<=n;i++) cnt[rank[i]]++;for(int i=1;i<=nn;i++) cnt[i]+=cnt[i-1];for(int i=n;i>=1;i--) sa[cnt[rank[sec[i]]]--]=sec[i];top=tmp[sa[1]]=1;for(int i=2;i<=n;i++){if(rank[sa[i]]!=rank[sa[i-1]]||rank[sa[i]+k]!=rank[sa[i-1]+k]) top++;tmp[sa[i]]=top;}nn=top;for(int i=1;i<=n;i++) rank[i]=tmp[i];}}int main(){scanf("%s",s+1);n=strlen(s+1);int t=n;for(int i=1;i<=n;i++) s[i+n]=s[i];n*=2;get_sa();n=t;for(int i=1;i<=n*2;i++) if(sa[i]<=n)printf("%c",s[sa[i]+n-1]);printf("\n");return 0;}


0 0
原创粉丝点击