后缀数组二倍增模板解析

原创 2015年07月09日 22:29:49

后缀数组模板解析

    自己的理解和分析 如有错误或不同的理解请指出。

    int cmp(int *r, int a, int b, int k)

    {

        return r[a] == r[b]&&r[a + k] == r[b + k];

    }

    void getsa(int *r, int *sa, int n, int m)

    {

        int i, j, p, *x = wa, *y = wb, *t;  

        for(i = 0;i < m;i++)

            wsf[i] = 0;

        for(i = 0;i < n;i++)

            wsf[x[i] = r[i]]++;

        for(i = 1;i < m;i++)

            wsf[i] += wsf[i - 1];

        for(i = n - 1;i >= 0;i--)

            sa[--wsf[x[i]]] = i;

      基基排序。(基数排序)

      X[i]数组记录的是一开始的r[i]的值。 

      Wsf[i]数组记录的是值为i和小于i的元素有几个。

      最后一个循环则是记录排第几的数的位置。

 

        j = 1;

        p = 1;

        for(;p < n;j *= 2, m = p)

        {

            for(p = 0, i = n - j;i < n;i++)

                y[p++] = i;

            for(i = 0;i < n;i++)

                 if(sa[i] >= j)

                     y[p++] = sa[i] - j;

      Y[i]记录的是排名为i的第二要素的第一要素位置。

      Sa[i]记录的是上一次排序的元素排名第i的元素的位置。

      因为是二倍增排序,j就是第二要素据第一要素的距离,所以n-j第一元素的第二元素都为零,所以排名直接从零开始。

      第二个for循环就是单纯的记录第二要素的排名。至于为什么sa[i]>=j是因为j前面的元素不肯能作为第二元素。

            for(i = 0;i < n;i++)

                wv[i] = x[y[i]];

            for(i = 0;i < m;i++)

                wsf[i] = 0;

            for(i = 0;i < n;i++)

                 wsf[wv[i]]++;

            for(i = 1;i < m;i++)

                 wsf[i] += wsf[i - 1];

            for(i = n - 1;i >= 0;i--)

                sa[--wsf[wv[i]]] = y[i];

    与一开始一样 是基基排序。至于第一个for循环则是根据第二元素进行了排序然后赋值。因为最后一个for循环是从后往前的所以先扫到的排名低,所以要构成一个新的数列。

    同样的 因为第一个for循环 所以最后一个for循环就是=y[i].因为y[i]表示的是第二元素排名第i位的第一元素的位置.sa[i]数组记录的是排名第i位的原数组的位置。

            t = x;

            x = y;

            y = t;

            x[sa[0]] = 0;

            for(i = 1, p = 1;i < n;i++)

                 x[sa[i]] = cmp(y, sa[i - 1], sa[i], j) ? p - 1: p++;

        }

    }

    数组交换.

    因为一开始人为的加入了x[末尾+1] = 0;0肯定是排在第0位的,所以x[sa[0]] = 0;只是再加一遍而已。这里p 的值表示的是不同关键值得数量。

    void getheight(int *r, int n)

    {

        int i, j, k = 0;

        for(i = 1; i <= n;i++)

             rank[sa[i]] = i;

      求出rank数组,rank[i]的意思是位置为i的后缀的排名,与sa[i]数组正好互逆。

      for(i = 0; i < n;i++)

      {

            if(k)

                k--;

            else

                k = 0;

            j = sa[rank[i] - 1];

            while(r[i + k] == r[j + k])

                k++;

            height[rank[i]] = k;

        }

    }

      J是位置为i的后缀排名前一名的位置。如果这两个字符相等边加一。

                                  

                                                                                             2015.7.9 by Eson

[Tyvj 1860] 后缀数组 倍增O(nlogn) 求sa[],height[]模板,单调栈维护,求特征值

#include #include #include using namespace std; const int maxn=200005; char text[maxn]; int SA[maxn]...
  • ALPS233
  • ALPS233
  • 2016年03月08日 15:06
  • 437

SA(后缀数组)——倍增法模板

倍增法原理:一步一步把n个后缀按照它们的前 2 ^ k 的字符排序。

后缀数组(倍增)【模板+注释】

看了两天的后缀数组,终于算是理解了。为了方便复习,加了些注释。

后缀数组倍增算法实现

  • 2012年09月02日 16:24
  • 2KB
  • 下载

后缀数组的倍增法实现

  • 2013年12月13日 20:24
  • 17KB
  • 下载

POJ 2774 Long Long Message(后缀数组:倍增算法)

字符串的任何一个子串都是该字符串某个后缀的前缀。 两个字符串的最长公共子串可转化为,求两者后缀两两比较的最长公共前缀的最大值。 定义: 后缀数组(sa[]),名次数组(rank[]): 子串st...

2倍倍增算法构造后缀数组

用2
  • kamsau
  • kamsau
  • 2014年09月15日 23:18
  • 834

后缀数组倍增法

那夜风,宫灯昏暗,小楼听雨灯辉摇梦

【倍增后缀数组】 poj1743 &2406& 3261 & 3294 & 2774

曾经有几次想写后缀数组,结果被论文中缩成豆腐块的程序吓到了.......         现在看懂了之后, 被这个程序写的丑陋吓到了。         丑虽丑,但是写的还是很精巧的, 照着打吧。 ...
  • cjoilmd
  • cjoilmd
  • 2012年01月15日 22:13
  • 628

POJ 3693 后缀数组 重复次数最多的连续重复子串 倍增法以及D3法

n以第k个字符开始的后缀称为后缀k。   【后缀数组SA】后缀数组保存的是一个字符串的所有后缀的排序结果。其中SA[i]保存的是字符串所有的后缀中第i小的后缀的开头位置。 sa记录的是排序...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:后缀数组二倍增模板解析
举报原因:
原因补充:

(最多只允许输入30个字)