后缀数组模版

<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，而在于写代码的能力和查询的算法实现。