后缀数组模板 SuffixArray

 1 #include<bits/stdc++.h>
 2 #define rep(i,a,b) for(i=a;i<b;++i)
 3 #define drep(i,a,b) for(i=b;i>=a;--i)
 4 using namespace std;
 5 int p[N],wc[N],wv[N],wa[N],wb[N],rank[N],sa[N];
 6 char s[N];
 7 inline bool cmp(int *r,int a,int b,int l) {
 8     return r[a]==r[b]&&r[a+l]==r[b+l];
 9 }
10 inline void tate(int *x,int *y,char *r,int n,int m,bool opt) {
11     int i;
12     rep(i,0,m) wc[i]=0;
13     rep(i,0,n) wc[opt?(x[i]=r[i]):wv[i]]++;
14     rep(i,1,m) wc[i]+=wc[i-1];
15     drep(i,n-1,0) sa[--wc[opt?x[i]:wv[i]]]=opt?i:y[i];
16 }
17 void f(char *r,int *sa,int n,int m) {
18     int i,j,p,*x=wa,*y=wb,*t;
19     tate(x,y,r,n,m,1);
20     for (j=1,p=1;p<n;j<<=1,m=p) {
21         p=0;
22         rep(i,n-j,n) y[p++]=i;
23         rep(i,0,n) if(sa[i]>=j) y[p++]=sa[i]-j;
24         rep(i,0,n) wv[i]=x[y[i]];
25         tate(x,y,r,n,m,0);
26         for (t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;++i)
27             x[sa[i]]=cmp(r,sa[i-1],sa[i],j)?p-1:p++;
28     }
29 }

压了下,把基数排序的压进函数,然后看起来还挺好写的,这只是后缀数组的一个阶段啦, 继续努力!

时间: 2024-10-11 06:49:25

后缀数组模板 SuffixArray的相关文章

后缀数组模板第一版

/*---------------倍增算法+RMQ后缀数组模板-------------- 输入:从0开始的字符串g,长度len最大为10^6 输出: sa[]表示:n 个后缀从小到大进行排序之后把排好序的后缀的开头位置顺 次放入 sa 中,sa[i]表示排第i位的字符串开头是sa[i],因为添加了一个结尾0,所以sa[0]=len height 数组(h[]):定义 h[i]=suffix(sa[i-1])和 suffix(sa[i])的最长公 共前缀,也就是排名相邻的两个后缀的最长公共前缀.

后缀数组模板一份

1 /****************** 2     by zhuyuqi      * 3     QQ:1113865149 * 4     name:2-sat    * 5                   * 6 ******************/ 7  8 using namespace std; 9 const int MAX = 1000;10 int r[MAX],*rank;11 int wa[MAX],wb[MAX],ws[MAX],wv[MAX];12 int h

后缀数组模板/LCP模板

1 //后缀数组模板,MANX为数组的大小 2 //支持的操作有计算后缀数组(sa数组), 计算相邻两元素的最长公共前缀(height数组),使用get_height(); 3 //计算两个后缀a, 和b的最长公共前缀,请先使用lcp_init(),再调用get_lcp(a, b)得到 4 //下面的n是输入字符串的长度+1(n = strlen(s) + 1), m是模板的范围 m=128表示在字母,数字范围内,可以扩大也可缩小 5 //s[len] 是插入的一个比输入字符都要小的字符 6 s

后缀数组模板及解释

以前做过后缀数组,直接用模板,最近打算重新认真的学一遍.感觉学一个东西一定要弄懂了,不然到最后还是要重学. int wa[MAXN],wb[MAXN],wv[MAXN],Ws[MAXN]; void da(int *r,int *sa,int n,int m){//n表示字符串长度 + 1,包括添加的那个0,m表示取值的范围 //把单个字符进行基数排序 int *x = wa,*y = wb; for(int i = 0; i < m; i++)Ws[i] = 0; for(int i = 0;

后缀数组模板及一些数组的含义

最近学习了一下后缀数组,模板原理以后再看,先记一下一些数组的含义.用以下这张图做例子: rank(i)代表第i个后缀的字典序排名 sa(i)代表排名为i的字典序对应的位置 lcp(i, j)表示suffix(i)和suffix(j)的公共最长前缀 height(i) = lcp(sa(i-1), sa(i)) 当rank(i)<rank(j),有lcp(i,j) = min(height(k)),   rank(i)<k<=rank(j) 原文地址:https://www.cnblogs

hdu1403(后缀数组模板)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1403 题意: 给出两个字符串, 求他们的最长公共子串 思路: 两个字符串的最长公共子串长度显然就是两个字符串的所有后缀中的最长公共前缀长度. 可以先用一个没有出现的字符(便于后面区分后缀是否属于相同字符串)将两个字符串连成一个字符串,再用后缀数组求其height, SA数组, 对于当前 i, 通过 SA 数组区分后缀 i 和 i - 1 是否在同一个初始字符串中, 若不是则用 height[i] 维

后缀数组 模板题 hdu1403(最长公共(连续)子串)

好气啊,今天没有看懂后缀树和后缀自动机 只能写个后缀数组发泄一下了orz #include <cstdio> #include <cstring> const int N = 100017*2; int wa[N], wb[N], wv[N], ws[N]; int rank[N]; int height[N]; char str[N]; int s[N], sa[N]; int cmp(int *r, int a, int b, int l) { return r[a]==r[b

后缀数组模板(理解)

字符串的处理真可谓是博大精深,后缀数组这种数据结构我花了两天时间才明白了其构造的过程.主要是代码不好理解. 数据结构: 1.sa数组,就是后缀数组,按照字典序排列,其意义为:sa[i]=k,排第i名的子串是从k位开始的. 2.rank名次数组,其意义为:rank[i]=k,以i为起点的子串排名为k. 很容易看出来两者可以相互转化. 求这两个数组的过程是基于基数排序,计数排序的方法. 下面是一个大牛的注释版.其中我在补充点: 1.对于求y[]这个数组,因为已经知道上次长为j的子串比较结果,那么这次

[UOJ#35] [UOJ后缀数组模板题] 后缀排序 [后缀数组模板]

后缀数组,解决字符串问题的有利工具,本题代码为倍增SA算法 具体解释详见2009年国家集训队论文 1 #include <iostream> 2 #include <algorithm> 3 #include <cstdio> 4 #include <cstdlib> 5 #include <cstring> 6 #include <cmath> 7 #include <ctime> 8 9 using namespace