【后缀数组之应用】【待续】
2015-08-31 11:10
387 查看
【最长重复子串问题】
可重叠最长重复子串 -- POJ 3261
Disc: 给出包含n个元素的数组a,问其中最长可重叠重复子串的长度,要求该子串至少重复k次;输入保证至少存在一个重复k次的最长子串;
Tips: 二分查找子串的长度,注意对该长度的子串是否存在K个重复子串的判定方法;
代码:
/* Problem: POJ 2774 Tips: 最长公共连续子串 Date: 2015.8.30 */ #include<iostream> #include<cstdio> #include<cstring> using namespace std; const int maxn = 200002; int sa[maxn], rank[maxn], height[maxn]; int wa[maxn], wb[maxn], wv[maxn], wd[maxn]; int cmp(int *r, int a, int b, int l){ return r[a] == r[b] && r[a+l] == r[b+l]; } void da(int *r, int n, int m){ // 倍增算法 r为待匹配数组 n为总长度 m为字符范围 int i, j, p, *x = wa, *y = wb, *t; for(i = 0; i < m; i ++) wd[i] = 0; for(i = 0; i < n; i ++) wd[x[i]=r[i]] ++; for(i = 1; i < m; i ++) wd[i] += wd[i-1]; for(i = n-1; i >= 0; i --) sa[-- wd[x[i]]] = i; for(j = 1, p = 1; 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; for(i = 0; i < n; i ++) wv[i] = x[y[i]]; for(i = 0; i < m; i ++) wd[i] = 0; for(i = 0; i < n; i ++) wd[wv[i]] ++; for(i = 1; i < m; i ++) wd[i] += wd[i-1]; for(i = n-1; i >= 0; i --) sa[-- wd[wv[i]]] = y[i]; for(t = x, x = y, y = t, p = 1, x[sa[0]] = 0, i = 1; i < n; i ++){ x[sa[i]] = cmp(y, sa[i-1], sa[i], j) ? p - 1: p ++; } } } void calheight(int *r, int n){ // 求height数组。 int i, j, k = 0; for(i = 1; i <= n; i ++) rank[sa[i]] = i; for(i = 0; i < n; height[rank[i ++]] = k){ for(k ? k -- : 0, j = sa[rank[i]-1]; r[i+k] == r[j+k]; k ++); } } void solve(int n, int l1, int l2) { int maxl = 0; for(int i = 2; i < n; i++) { if(height[i] > maxl) { if(sa[i-1] >= 0 && sa[i-1] < l1 && sa[i] > l1) maxl = height[i]; if(sa[i] >= 0 && sa[i] < l1 && sa[i-1] > l1) maxl = height[i]; } } printf("%d\n", maxl); } int main() { char s1[maxn], s2[maxn]; int s[maxn], r[maxn]; scanf("%s%s", s1, s2); int l1 = strlen(s1), l2 = strlen(s2); int n = 0; for(int i = 0; i < l1; i++) s[n++] = (s1[i]-'a'+1); s[n++] = 28; for(int i = 0; i < l2; i++) s[n++] = (s2[i]-'a'+1); s = 0; da(s,n+1,30); calheight(s,n); solve(n, l1, l2); return 0; }View Code
相关文章推荐
- 解决firefox弱临时 Diffie-Hellman 密钥的错误
- HDU 3743 Frosh Week(逆序对-BIT)
- hdoj 4322 Candy 【最大费用最大流】【经典题目】【最大流时 维护费用的最大效益】
- liunx-fastboot命令行的使用方法
- Hadoop常见错误及解决办法汇总
- 构造url测试action与dao实例
- 【笔试】37、顺时针打印矩阵
- 数据库常用函数大全(数学函数、字符串函数、其他函数)
- 二叉搜索树的后序遍历序列
- Linux:chmod -R 777 * 是什么意思?
- SQL ---指令实例语句
- 初到Hp实训基地感受
- WEB中配置错误页面
- 详解C语言中的wait()函数和waitpid()函数
- UCenter应用同步登录分析
- Java数组学习练习3冒泡排序数组代码
- LearningSpark6.1
- ARM 的指令集
- 使用Python的requests库进行接口测试——session对象的妙用
- 关于改变Bitmap透明度