题目描述 Description
天凯是MIT的新生。Prof. HandsomeG给了他一个长度为n的由小写字母构成的字符串,要求他把该字符串的n个后缀(suffix)从小到大排序。
何谓后缀?假设字符串是S=S1S2……Sn,定义Ti=SiSi+1……Sn。T1, T2, …, Tn就叫做S的n个后缀。
关于字符串大小的比较定义如下(比较规则和PASCAL中的定义完全相同,熟悉PASCAL的同学可以跳过此段):
若A是B的前缀,则A<B;否则令p满足:A1A2…Ap-1=B1B2…Bp-1,Ap<>Bp。如果Ap<Bp,则A<B;否则A>B。
输入描述 Input Description
第一行一个整数n(n<=15000)
第二行是一个长度为n字串。
输出描述 Output Description
输出文件包含n行,第i行是一个整数pi。表示所有的后缀从小到大排序后是Tp1, Tp2, …, Tpn。
样例输入 Sample Input
4
abab
样例输出 Sample Output
3
1
4
2
数据范围及提示 Data Size & Hint
说明:后缀排序后的顺序是T3=”ab”, T1=”abab”, T4=”b”, T2=”bab”。所以输出是3, 1, 4, 2。
1 /* 2 后缀数组板子题,注意第29行写法 3 */ 4 include<iostream> 5 #include<cstdio> 6 #include<cstring> 7 #define maxn 15010 8 using namespace std; 9 int m=26,n,sa[maxn],rank[maxn],c[maxn],s[maxn],x[maxn],y[maxn]; 10 char ch[maxn]; 11 int main(){ 12 freopen("Cola.txt","r",stdin); 13 scanf("%d%s",&n,ch); 14 for(int i=0;i<n;i++)s[i]=ch[i]-'a'; 15 16 for(int i=0;i<n;i++)c[x[i]=s[i]]++; 17 for(int i=1;i<m;i++)c[i]+=c[i-1]; 18 for(int i=n-1;i!=-1;i--)sa[--c[x[i]]]=i; 19 int k=1,p=0; 20 while(k<=n){ 21 for(int i=n-k;i<n;i++)y[p++]=i; 22 for(int i=0;i<n;i++)if(sa[i]>=k)y[p++]=sa[i]-k; 23 for(int i=0;i<m;i++)c[i]=0; 24 for(int i=0;i<n;i++)c[x[y[i]]]++; 25 for(int i=1;i<m;i++)c[i]+=c[i-1]; 26 for(int i=n-1;i!=-1;i--)sa[--c[x[y[i]]]]=y[i]; 27 swap(x,y);p=1;x[sa[0]]=0; 28 for(int i=1;i<n;i++) 29 if(y[sa[i-1]]==y[sa[i]]&&((y[sa[i-1]+k]==y[sa[i]+k]&&sa[i-1]+k<n&&sa[i]+k<n)||(sa[i-1]+k>=n&&sa[i]+k>=n))) x[sa[i]]=p-1; 30 else x[sa[i]]=p++; 31 if(p>=n)break; 32 m=p;p=0;k*=2; 33 } 34 for(int i=0;i<n;i++)printf("%d ",sa[i]+1); 35 }