后缀数组模版

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