uoj35: 后缀排序

原创 2017年01月03日 10:55:32

题目链接
后缀数组模板题。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

const int N = 100010;

int n, c[N], sa[N];
int rank[N], height[N];
char str[N];

inline void buildSa(int m) {
    int *x = rank, *y = height;
    for (int i = 0; i < m; i++) c[i] = 0;
    for (int i = 0; i < n; i++) c[x[i] = str[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;
    }
}
inline void calcHeight() {
    int k = 0;
    for (int i = 0; i < n; i++) rank[sa[i]] = i;
    for (int i = 0; i < n - 1; i++) {
        if (k) k--;
        int j = sa[rank[i] - 1];
        while (str[i + k] == str[j + k]) k++;
        height[rank[i]] = k;
    }
}

int main() {
    scanf("%s", str);
    n = strlen(str);
    for (int i = 0; i < n; i++)
        str[i] -= 'a' - 1;
    str[n++] = 0;//补$不能省略
    buildSa(30);
    calcHeight();
    for (int i = 1; i < n; i++)
        printf("%d ", sa[i] + 1); puts("");
    for (int i = 2; i < n; i++)
        printf("%d ", height[i]); puts("");
    return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

[UOJ35]后缀排序

后缀自动机 后缀树 后缀数组
  • hbhcy98
  • hbhcy98
  • 2016年04月07日 19:16
  • 238

uoj#35: 后缀排序

题意: 给出一个字符串; 求这个字符串的后缀数组和height数组; 字符串长度 题解: 后缀数组纯裸题; 不得不说众人的力量是伟大的,VFK的毒瘤功底是有目共睹的; A掉这道题...

uoj P35 后缀排序

传送门 妥妥的后缀数组裸题(主要是楼主菜只会打这种题目),还有拓展lucas还没调对qwq菜飞了。。。 求sa和lcp,就是后缀数组求sa和height。 代码: #include #i...

Uoj#35. 后缀排序

#include #include #include #include #include #include #include #include #include #include #define T ...
  • ws_yzy
  • ws_yzy
  • 2016年01月10日 16:48
  • 342

uoj#35. 后缀排序 后缀数组

uoj#35. 后缀排序 裸的后缀数组…… 敲了2遍…… 希望……不要忘记QAQ 参考了lrj代码和某个AC神犇参考了lrj代码后写的代码……#include #include #includ...
  • yxr0105
  • yxr0105
  • 2016年02月04日 00:10
  • 188

UOJ #35 后缀排序(后缀数组)

题目链接 题解:后缀数组板子 之前学后缀数组只是学了个大概,而且写的是二维的。现在从头理一遍,改成了算法导论上的写法。。 #include #include #include #include #...

UOJ #35. 后缀排序(后缀数组模板题)

测试一下模板 #include #include using namespace std; const int maxn =2e5+10; struct Suffix { int r[ma...

Uoj #35. 后缀排序(后缀数组)

35. 后缀排序 统计 描述 提交 自定义测试 这是一道模板题。 读入一个长度为 nn 的由小写英文字母组成的字符串,请把这个字符串的所有非空后缀按字典序...

BZOJ 3676 UOJ 103 APIO 2014 后缀自动机 Manacher

对于所有的本质不同的回文子串,在SAM中查找对应节点,由于我们已知回文子串的右端点,因此可以在Parent树中快速定位我们要的状态。3676: [Apio2014]回文串Time Limit: 20 ...

UOJ #9 vfk的数据 排序 字典序比较

UOJ #9 vfk的数据 排序题目描述UOJ上,由于一群人开心的在hack着,有的题的测试数据都有1000010000组了!为了处理这些数据,vfk决定把这道题所有输入数据都弄出来,每个数据给个编号...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:uoj35: 后缀排序
举报原因:
原因补充:

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