#1407 : 后缀数组二·重复旋律2 (不可重叠最长重复子串问题)
题目链接
思路
求不可重叠最长重复子串,也可以利用height[i]height[i]height[i]。
- 二分枚举答案KKK
- 当height[i]≥Kheight[i] \ge Kheight[i]≥K可能存在解,然后判断min(sa[i],sa[i−1]),max(sa[i],sa[i−1])min(sa[i], sa[i-1]),max(sa[i], sa[i-1])min(sa[i],sa[i−1]),max(sa[i],sa[i−1])的关系,看能否构造成不重叠的公共前缀
- height[i]height[i]height[i]数值相同的分成一组,取小组中的最值计算
#include <bits/stdc++.h>
#define LL long long
#define P pair<int, int>
#define lowbit(x) (x & -x)
#define mem(a, b) memset(a, b, sizeof(a))
#define mid ((l + r) >> 1)
#define lc rt<<1
#define rc rt<<1|1
#define endl '\n'
const int maxn = 1e5 + 5;
const int inf = 0x3f3f3f3f;
const int mod = 1e9 + 7;
using namespace std;
int cntA[maxn], cntB[maxn], A[maxn], B[maxn];
int Sa[maxn], tsa[maxn], height[maxn], Rank[maxn];
int s[maxn];
int n;
void suffixArray () {for (int i = 0; i < 1000; ++i) cntA[i] = 0;for (int i = 1; i <= n; ++i) cntA[s[i]]++;for (int i = 1; i < 1000; ++i) cntA[i] += cntA[i-1];for (int i = n; i >= 1; --i) Sa[ cntA[s[i]]-- ] = i;Rank[ Sa[1] ] = 1;for (int i = 2; i <= n; ++i) {Rank[Sa[i]] = Rank[Sa[i-1]];if (s[Sa[i]] != s[Sa[i-1]]) Rank[Sa[i]]++;}for (int l = 1; Rank[Sa[n]] < n; l <<= 1) {for (int i = 0; i <= n; ++i) cntA[i] = 0;for (int i = 0; i <= n; ++i) cntB[i] = 0;for (int i = 1; i <= n; ++i) {cntA[ A[i] = Rank[i] ]++;cntB[ B[i] = (i + l <= n) ? Rank[i+l] : 0]++;}for (int i = 1; i <= n; ++i) cntB[i] += cntB[i-1];for (int i = n; i >= 1; --i) tsa[ cntB[B[i]]-- ] = i;for (int i = 1; i <= n; ++i) cntA[i] += cntA[i-1];for (int i = n; i >= 1; --i) Sa[ cntA[A[tsa[i]]]-- ] = tsa[i];Rank[ Sa[1] ] = 1;for (int i = 2; i <= n; ++i) {Rank[Sa[i]] = Rank[Sa[i-1]];if (A[Sa[i]] != A[Sa[i-1]] || B[Sa[i]] != B[Sa[i-1]]) Rank[Sa[i]]++;}}for (int i = 1, j = 0; i <= n; ++i) {if (j) --j;int tmp = Sa[Rank[i] - 1];while (i + j <= n && tmp + j <= n && s[i+j] == s[tmp+j]) ++j;height[Rank[i]] = j;}
}
int check(int x) {int minsa = 1e9, maxsa = -1e9;for (int i = 1; i <= n; ++i) {if (height[i] >= x) {minsa = min(minsa, Sa[i]);maxsa = max(maxsa, Sa[i]);if (maxsa - minsa >= x) return 1;}else {minsa = Sa[i];maxsa = Sa[i];}}return 0;
}int main () {ios::sync_with_stdio(0);cin.tie(0), cout.tie(0);cin >> n;for (int i = 1; i <= n; ++i) {cin >> s[i];}suffixArray();int l = 1, r = n;while (l <= r) {if (check(mid)) l = mid + 1;else r = mid - 1;} cout << r << endl;return 0;
}
#1407 : 后缀数组二·重复旋律2 (不可重叠最长重复子串问题)相关推荐
- 后缀数组--(可重叠最长重复子串问题)
问题描述:给定一个字符串,求最长重复子串,这两个子串可以重叠. 其实问题可以转化为height数组的最大值.至于为什么是这样,我可以这样解释: 求可重叠最长重复子串等价于求两个后缀的最长公共前缀的最大 ...
- POJ 1743 (后缀数组+不重叠最长重复子串)
题目链接: http://poj.org/problem?id=1743 题目大意:楼教主の男人八题orz.一篇钢琴谱,每个旋律的值都在1~88以内.琴谱的某段会变调,也就是说某段的数可以加减一个旋律 ...
- I - 后缀数组二·重复旋律2 HihoCoder - 1407
描述 小Hi平时的一大兴趣爱好就是演奏钢琴.我们知道一个音乐旋律被表示为长度为 N 的数构成的数列.小Hi在练习过很多曲子以后发现很多作品自身包含一样的旋律. 旋律可以表示为一段连续的数列,相似的旋律 ...
- #1403 : 后缀数组一·重复旋律 (可重叠最长重复K次子串问题)
题目链接 思路 K=2K = 2K=2时就是求height[i]height[i]height[i]的最大值. K≠2K =\not 2K≠2时就是连续K−1K-1K−1个height[i]he ...
- poj1743(后缀数组+二分--不可重叠最长重复子串)
题目:Musical Theme #include <stdio.h> #include <string.h> #define N 1000010 int wa[N],wb[N ...
- 字符串-后缀树和后缀数组详解
文章目录 后缀树 后缀数组 概念 sa[] rk[] height[] 例题 HDU-1403最长公共子串 洛谷P2408 不同子串个数 HDU-5769Substring 后缀树 建议先了解一下字典 ...
- 后缀数组求最长重复子串
问题描述 给定一个字符串,求出其最长重复子串 例如:abcdabcd 最长重复子串是 abcd,最长重复子串可以重叠 例如:abcdabcda,这时最长重复子串是 abcda,中间的 a 是被重叠的. ...
- 字符串相关处理kmp,前缀数,后缀树,后缀数组,最长回文串,最长重复字串,最长非重复字串
1. 最长回文串 一般用后缀数组或者后缀树可以解决, 用此方法:http://blog.csdn.net/v_july_v/article/details/6897097 预处理后缀树,使得查询LCA ...
- 试题 算法训练 后缀数组——最长重复子串
资源限制 时间限制:100ms 内存限制:256.0MB 问题描述 给定一个长度为n的数串,求至少出现k 次的最长重复子串的长度,这k 个子串可以重叠.保证有子串出现至少k次. 输入格式 第一行:两个 ...
最新文章
- Java集合的概念及其相关类的简单使用
- [原][osg][gdal]两种方式修改tiff高程
- Microsoft Office SharePoint Server 2007 Starter Page
- 聚类Introducion
- java swing removeall_Java Swing Combobox removeAllItems调用ItemStateChange...
- pandas读取大文件csv容易出现的错误
- Cocos2dx游戏开发系列笔记8:开搞一个射击游戏《战神传说》//就个打飞机的
- 建立域用户时遇到的问题。
- 误删数据库怎么办?mysql 回滚,撤销操作,恢复数据
- 温度冲击试验箱军标GJB 150与GB 2423区别
- LeCun、Bengio、Hinton三巨头曾合体,Nature发文综述深度学习(论文精华)
- 电信dns地址是多少?
- WORD邮件合并打印EXCEL数据制作大量奖证、奖状、准考证、成绩单、明信片、信封等个人报表
- 【原创】超级详细的iPhone铃声制作、导入、设置全过程教程
- 上海交通大学电子与通信工程819考研上岸经验分享
- 重庆航天职业技术学院计算机系在哪个校区,重庆航天职业技术学院地址在哪里,哪个城市...
- MOSFET开通特性(1)——输入特性
- pycharm启动图片修改
- Python 编程训练2:二维向量
- 拉美外贸: 哥伦比亚网银支付PSE
热门文章
- idea maven 创建webapp项目没有src目录
- VS.net中快捷键收缩和展开代码段
- hdu 3047 Zjnu Stadium(并查集)
- winform TreeView 节点选择
- Linux系统管理技术手册——第6章 添加新用户
- 继承练习 :开发一个系统时 	需要对员工进行建模 员工包含3个属性 姓名 工号 工资 功能 work
- [android] No resource found that matches the given name 'Theme.AppCompat.Light'
- Windows 10每天定时执行脚本
- 信息学奥赛一本通(C++)在线评测系统——基础(二)基础算法 —— 1312:【例3.4】昆虫繁殖
- 信息学奥赛一本通(C++)在线评测系统——基础(一)C++语言——1010:计算分数的浮点数值