<pre name="code" class="cpp">#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=200090;
int n,height[N],sa[N],s[N],rank[N],c[N],t1[N],t2[N];
char ch[N];
void build_sa()
{
int m=27,*x=t1,*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<=n;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);
x[sa[0]]=0;p=1;
for (int i=1;i<n;i++)
x[sa[i]]= y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?p-1:p++;
if (p>n) break;
m=p;
}
}
void build_height()
{
int j,k=0;
for (int i=0;i<n;i++) rank[sa[i]]=i;
for (int i=0;i<n-1;i++)
{
if (k) k--;
j=sa[rank[i]-1];
while (s[j+k]==s[i+k]) k++;
height[rank[i]]=k;
}
}
int main()
{
scanf("%s",ch);
n=strlen(ch);
for (int i=0;i<n;i++) s[i]=ch[i]-'a'+1;
s[n++]=0;//最小字符
build_sa();
build_height();
for (int i=1;i<n;i++) printf("%d ",sa[i]+1);//+1,因为数组是从0开始的
//从一开始,因为,第0名肯定是最后一个那个最小字符
printf("\n");
for (int i=1;i<n;i++) printf("%d ",height[i]);
//排在第i名的和第i-1名的最长公共前缀
return 0;
}
最后给出几题习题(改天要练)
•Poj 2774 –最长公共连续子串,入门题目
•Poj1743—最长不重叠重复子串
•Hint:二分的判定要小心点,这题有点特别。
•Poj3294—出现次数超过一半的最长子串
•Hint:判断组中不同串出现次数的技巧很关键
•Poj3261—重复k次可重叠子串。
•Hint:会了上面两题,这题应该很简单,可以试试用单调栈。
•Poj2758—后缀数组+rmq
•Hint:这题难度不在rmq,而在于写代码的能力和查询的算法实现。