【后缀数组学习中】

来源:互联网 发布:淘宝男装货到付款 编辑:程序博客网 时间:2024/05/16 12:10

小白书 219 页入门


大概懂得后缀数组原理 即后缀排序后的下标所形成的数组


裸的算法是 n*n*logn:直接利用定义 对所有后缀快排........


稍微好点. nlogn的倍增算法 原理懂了


更好点 线性时间的...只是简单涉及了没有


逛逛博客 学习前人优秀代码风格


http://dongxicheng.org/structure/suffix-array/

很好的将几篇国家队论文整合了,题目,介绍都有。好好研究


http://blog.csdn.net/xymscau/article/details/8798046

提供了代码和几个练习题

http://www.cnblogs.com/kuangbin/archive/2013/05/20/3089189.html

kuangbing写的

const int MAXN=300000;char s[MAXN];int sa[MAXN],t[MAXN],t2[MAXN],c[MAXN],n;void build_sa(int m){    int *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];    for(int i=n-1 ; i >= 0;i--) sa[--c[x[i]]] = i;    for(int k = 1 ; k <= n; k <<=1) {        int p=0;        //直接利用sa数组排序第二关键字        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 = 0; i < m; i++)   c[i] += c[i-1];        for(int i = n-1;i >= 0; i--) sa[--c[x[y[i]]]]=y[i];        //根据sa和y数组重新计算新的y数组        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;    }}

以上代码还需好好研究.

---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

罗穗骞

后缀数组——处理字符串的有力工具

仔细研读后缀数组基本能全懂了


一.后缀数组代码

1.为什么要在最后加一个很小的在最后呢?(在第九页解释了)

这里可以看到规定s[n]=0 的好处,如果s[a]=s[b] 说明以s[a]或s[b]开头的长度为l的字符串肯定不包括字符s[n]

所以调用s[a+l] 和s[b+l] 不会导致数组下标越界,这样就不需要做特殊判断。


二.后缀数组的应用

2.1 最长公共前缀

height 数组:定义 height[i]= suffix(sa[i-1]) 和 suffix(sa[i]) 的最长公共前缀,也就是排名相邻的两个后缀的最长公共前缀。

0 0
原创粉丝点击