乐曲主题Musical Themes
SA例题
对于串 \(S\) 的两个子串 \(A\) 和 \(B\) ,满足 \(k = |A| = |B|\),\(\exists c \forall i\, a_i + c=b_i\),且 \(A\) 和 \(B\) 在原串中没有重叠,求最大的满足条件的 \(k\)
先考虑 \(\forall i\, a_i=b_i\) 的情况,二分 \(k\) ,如果在\(height\)上存在一段满足 \(\forall L \le i \le R \,height_i\ge k\),且存在 \(SA[p_1]-SA[p_2] \ge k\),说明 \(k\) 满足条件
接下来将上述做法用在原数组的差分数组上即可,注意判定时条件要改为\(SA[p_1]-SA[p_2] > k\),且要在第一个字符前添加特殊字符,否则判定时可能出问题
#include<cstdio>
using namespace std;
const int MAXN = 5001;
int str[MAXN], sa[MAXN], rank[MAXN], ht[MAXN], buc[MAXN];
void Swap(int*& a, int*& b){
int *tmp = a; a = b; b = tmp;
}
void getSA(int s[], int len){
int *rk = rank, *tp = ht;
int crd = 180, p;
for (int i = 1; i <= len; ++i)
rk[i] = s[i], tp[i] = i;
for (int i = 0; i <= crd; ++i)
buc[i] = 0;
for (int i = 1; i <= len; ++i)
++buc[rk[i]];
for (int i = 1; i <= crd; ++i)
buc[i] += buc[i - 1];
for (int i = len; i >= 1; --i)
sa[buc[rk[tp[i]]]--] = tp[i];
for (int w = 1; p != len; w <<= 1, crd = p){
p = 0;
for (int i = len - w + 1; i <= len; ++i)
tp[++p] = i;
for (int i = 1; i <= len; ++i)
if (sa[i] > w)
tp[++p] = sa[i] - w;
for (int i = 0; i <= crd; ++i)
buc[i] = 0;
for (int i = 1; i <= len; ++i)
++buc[rk[i]];
for (int i = 1; i <= crd; ++i)
buc[i] += buc[i - 1];
for (int i = len; i >= 1; --i)
sa[buc[rk[tp[i]]]--] = tp[i];
Swap(rk, tp);
rk[sa[1]] = p = 1;
for (int i = 2; i <= len; ++i)
rk[sa[i]] = tp[sa[i]] == tp[sa[i - 1]] && tp[sa[i] + w] == tp[sa[i - 1] + w] ? p : ++p;
}
for (int i = 1; i <= len; ++i)
rank[sa[i]] = i;
}
void getHeight(int s[], int len){
ht[1] = 0;
for (int i = 1, j = 0, pos; i <= len; ++i, --j){
if (rank[i] == 1)
continue;
if (j < 0)
j = 0;
pos = sa[rank[i] - 1];
while (s[i + j] == s[pos + j] && i + j <= len && pos + j <= len)
++j;
ht[rank[i]] = j;
}
}
bool judge(int k, int len){
int lp, rp;
for (int i = 1, j; i < len;){
if (ht[i + 1] >= k){
lp = rp = sa[i];
j = i + 1;
while (ht[j] >= k && j <= len){
if (sa[j] < lp)
lp = sa[j];
if (sa[j] > rp)
rp = sa[j];
++j;
}
if (rp - lp > k)
return 1;
i = j;
}
else
++i;
}
return 0;
}
int main(){
int slen;
scanf("%d", &slen);
for (int i = 1, pre = 89, tmp; i <= slen; ++i){
scanf("%d", &tmp);
str[i] = tmp - pre + 89; pre = tmp;
}
getSA(str, slen);
getHeight(str, slen);
int L = 1, R = slen >> 1, mid;
while (L + 1 < R){
mid = L + R >> 1;
if (judge(mid, slen))
L = mid;
else
R = mid;
}
int ans;
if (judge(R, slen))
ans = R;
else
ans = L;
if (ans < 4)
printf("0");
else
printf("%d", ans + 1);
return 0;
}
乐曲主题Musical Themes的更多相关文章
- 洛谷P2743 乐曲主题Musical Themes [USACO5.1] SA
正解:SA 解题报告: 传送门 这题三个条件嘛,那就一个个考虑下都解决了就把这题解决了嘛QwQ 那就直接分别针对三个条件写下各个击破就欧克辣? 1)长度大于等于5:求出答案之后和5比大小 2)不能有公 ...
- Luogu P2743 [USACO5.1]乐曲主题Musical Themes
链接 \(Click\) \(Here\) 人生第一道后缀数组的题目.首先要对输入的串进行差分处理,差分后长度为(\(n - 1\))的相同子段就是原串中长度为\(n\)的相同(可变调)子段.求出来\ ...
- [USACO5.1] 乐曲主题Musical Themes
题目链接:戳我 Emmm......hash怎么做啊不会啊 这里是SA后缀数组版本的 就是先两两做差分,作为要处理后缀的数组.普通地求出来h数组之后,我们二分这个答案,然后判定是否合法就行了.是否合法 ...
- P2743(poj1743) Musical Themes[差分+后缀数组]
P2743 乐曲主题Musical Themes(poj1743) 然后呢这题思路其实还是蛮简单的,只是细节特别多比较恶心,忘记了差分带来的若干疏漏.因为转调的话要保证找到相同主题,只要保证一段内相对 ...
- [转]jQuery EasyUI 扩展-- 主题(Themes)
主题(Themes)允许您改变站点的外观和感观.使用主题可以节省设计的时间,让您腾出更多的时间进行开发.您也可以创建一个已有主题的子主题. 主题生成器(Theme Builder) jQuery UI ...
- USACO Section 5.1 Musical Themes(枚举)
直接枚举O(n^3)会TLE,只要稍微加点优化,在不可能得到更优解时及时退出.其实就是道水题,虽说我提交了6次才过= =..我还太弱了 -------------------------------- ...
- USACO 5.1 Musical Themes(哈希+二分)
Musical ThemesBrian Dean A musical melody is represented as a sequence of N (1 <= N <= 5000) n ...
- Android应用开发中的风格和主题(style,themes)
http://www.cnblogs.com/playing/archive/2011/04/01/2002469.html 越来越多互联网企业都在Android平台上部署其客户端,为了提升用户体验, ...
- [USACO 5.1.3]乐曲主题
Description 我们用N(1 <= N <=5000)个音符的序列来表示一首乐曲,每个音符都是1..88范围内的整数,每个数表示钢琴上的一个键.很不幸这种表示旋律的方法忽略了音符的 ...
随机推荐
- Don't Panic! KRACK 没你想象的那么糟
上海交通大学密码与计算机安全实验室(LoCCS)软件安全小组(GoSSIP)版权所有,转载请与作者取得联系! 著名的计算机学术安全会议CCS在2017年录用了一篇名为Key Reinstallatio ...
- Python3.2官方文档翻译--作用域和命名空间实例
6.2.1 作用域和命名空间实例 以下的实例主要用来示范怎样引用不同的作用域和命名空间,keywordglobal和nonlocalru怎样影响变量绑定. 实例执行结果是: After local a ...
- nyoj 1077 小博弈 【另类巴什博奕】
分析:分析当整除(a+b)的时候肯定是后者胜利,假设余数不等于0的时候.假设余数大于b肯定是前者胜利,否则后者胜利. 代码: import java.math.*; import java.util. ...
- react webapp 开发小结
1.监听props的方法 componentWillReceiveProps(nextProps) { // } 2.监听state的方法 3.props 传递的方法 <AlarmList {. ...
- Odoo event
使用流程 建立活动 发布到网站 在线销售 订单确认,付款确认 注册.出席 建立活动 设置门票 确认并发布到网站 进入编辑模式,即可在线编辑活动 ...
- 开关电路_MOS和三极管
https://blog.csdn.net/acelit/article/details/70171312 绍过一般的电源开关电路,控制电源的目的是省电,控制静态电流.不过以下的电路存在着几个缺点: ...
- Linux kernel Wikipedia
http://en.wikipedia.org/wiki/Linux_kernel Development model The current development model of the Lin ...
- 实习日记)select option 选择不同的option时, 页面发生不同的变化
怎么在下拉框的选择不同的option时, 页面发生响应的变化 因为option是没有点击事件什么的, 只有select才有, 所以不能通过option的点击事件来完成, 所以开始的尝试都失败了(之前 ...
- Sass编译css/Grunt压缩文件
Sass安装(mac) $ sudo gem install sass scss编译成css文件 $ sass ui.scss ui.css CLI安装 $ npm install -g grunt- ...
- Mvc之Ajax上传图片
MyAjaxForm下载地址,点击此处下载 视图部分: @{ ViewBag.Title = "Index"; Layout = "~/Views/Shared/_Lay ...