验证中...
10月20日,北京【人工智能】源创会火热报名中,点击报名收获 AI 训练指南
【BZOJ1692】【USACO2007】队列变换
原始数据 复制代码
#include<bits/stdc++.h>
using namespace std;
const int N=400010;
int x[N],c[N],y[N],sa[N],rk[N];
int s[N];
void get_sa(int n,int m){
for(int i=1;i<=n;i++) c[x[i]=s[i]]++;
for(int i=2;i<=m;i++) c[i]+=c[i-1];
for(int i=n;i;i--) sa[c[x[i]]--]=i;
for(int k=1;k<=n;k<<=1){
int p=0;
for(int i=n-k+1;i<=n;i++) y[++p]=i;
for(int i=1;i<=n;i++) if(sa[i]>k) y[++p]=sa[i]-k;
for(int i=1;i<=m;i++) c[i]=0;
for(int i=1;i<=n;i++) c[x[i]]++;
for(int i=2;i<=m;i++) c[i]+=c[i-1];
for(int i=n;i;i--) sa[c[x[y[i]]]--]=y[i];
swap(x,y);
x[sa[1]]=1,p=1;
for(int i=2;i<=n;i++)
x[sa[i]]=(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k])?p:++p;
if(n==p) break;
m=p;
}
for(int i=1;i<=n;i++) rk[sa[i]]=i;
}
bool cmp(int a,int b){return s[a]==s[b]?rk[a]<rk[b]:s[a]<s[b];}
int main(){
int n;
scanf("%d",&n);
for(int i=1;i<=n;i++){
s[i]=getchar();
while(!isalpha(s[i])) s[i]=getchar();
}
s[1+n]='#';
int n1=(n<<1)+2;
for(int i=n+2;i<n1;i++) s[i]=s[n-i+n+2];
s[n1]='#';
get_sa(n1,255);
int l=1,r=n+2;
for(int i=1;i<=n;i++){
if(cmp(l,r)){
putchar(s[l]);
l++;
}else{
putchar(s[r]);
r++;
}
if(i%80==0) putchar('\n');
}
return 0;
}

评论列表( 0 )

你可以在登录后,发表评论

搜索帮助