#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 100005
using namespace std;
char s[maxn];
int sa[maxn],t[maxn],t2[maxn],c[maxn],Rank[maxn],height[maxn],n;
void build_sa(int m){
int i,*x=t,*y=t2,*T,p ;
n++;
for(i=0;i<m;++i)c[i]=0;
for(i=0;i<n;++i)++c[x[i]=s[i]];
for(i=1;i<m;++i)c[i]+=c[i-1];
for(i=n-1;i>=0;--i)sa[--c[x[i]]]=i;
for(int k=1;k<=n;k<<=1)
{
p=0;
for(i=n-1;i>=n-k;--i)y[p++]=i;
for(i=0;i<n;++i)if(sa[i]>=k)y[p++]=sa[i]-k;
for(i=0;i<m;++i)c[i]=0;
for(i=0;i<n;++i)++c[x[y[i]]];
for(i=1;i<m;++i)c[i]+=c[i-1];
for(i=n-1;i>=0;--i)sa[--c[x[y[i]]]]=y[i];
swap(x,y);
x[sa[0]]=0;p=1;
for(i=1;i<n;++i)
x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
if(p>=n)break;
m=p;
}
n--;
for(int i=1;i<=n;++i)printf("%d ",sa[i]+1);
printf("\n");
}
void LCP(){
int i,j,k=0;
for(i=1;i<=n;i++)Rank[sa[i]]=i;
for(int i=0;i<n;++i)
{
j=sa[Rank[i]-1];//h[i-1]
if(k)k--;
while(s[i+k]==s[j+k])k++;
height[Rank[i]]=k;//h[i]
}
for(int i=2;i<=n;++i)printf("%d ",height[i]);
}
int main(){
scanf("%s",s);
n=strlen(s);
build_sa(256);
LCP();
return 0;
}
后缀数组
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- Openjudge原题链接 题意输入一个长度N(<=20000)的数组，求出其重复K次的最长可重叠子串 思路由于N...
- 先%罗DADA建议按照论文手推，更易明白再%kuangbin大神 1、什么是后缀数组 后缀数组是后缀树的替代品，十...
