后缀数组模版

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<string>

using namespace std;
const int maxn=100009;
int s[maxn],sa[maxn],len,c[maxn],t[maxn],t2[maxn];

void build_sa(int *s,int *sa,int n,int m)//s 原数组，sa 后缀数组，n，原数组长度，m基数排序，数字的上线
{
	int i,*x=t,*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]+c[i];
	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]+c[i];
		for (int i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
		swap(x,y);
		p=1;
		x[sa[0]]=0;
		for (int 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;
	}
}
int main()
{
	cin>>len;
	for (int i=0;i<len;i++)
	{
		char ch;
		cin>>ch;
		s[i]=ch-'a'+1;
	}
	s[len]=0;//末尾加个最小的字符
	build_sa(s,sa,len+1,26);//末尾字符，所以len+1
	//for (int i=1;i<=len;i++) cout<<sa[i]<<' ';
	return 0;
}

最后给出几题习题（改天要练）

•Poj 2774 –最长公共连续子串，入门题目

•Poj1743—最长不重叠重复子串

•Hint：二分的判定要小心点，这题有点特别。

•Poj3294—出现次数超过一半的最长子串

•Hint：判断组中不同串出现次数的技巧很关键

•Poj3261—重复k次可重叠子串。

•Hint：会了上面两题，这题应该很简单，可以试试用单调栈。

•Poj2758—后缀数组+rmq

•Hint：这题难度不在rmq，而在于写代码的能力和查询的算法实现。