博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
洛谷 P4093: bzoj 4553: [HEOI2016/TJOI2016]序列
阅读量:5174 次
发布时间:2019-06-13

本文共 2149 字,大约阅读时间需要 7 分钟。

题目传送门:。

题意简述:

给定一个长度为 \(n\) 的序列 \(a\)。

同时这个序列还可能发生变化,每一种变化 \((x_i,y_i)\) 对应着 \(a_{x_i}\) 可能变成 \(y_i\)。

不会同时发生两种变化。

需要找出一个最长的子序列,使得这个子序列在任意一种变化下都是不降的。

只需要求出这个子序列的长度即可。

注意:可以不发生任何变化。

题解:

记 \(f[i]\) 为以第 \(i\) 项结尾的子序列最长长度。

则有转移:\(f[i]=\max_{j<i}(f[j])+1\),同时还要满足 \(maxval_j\le a_i\) 和 \(a_j\le minval_i\)。

按照项从小到大转移,形成了天然的时间顺序,同时还要满足两个偏序限制。

其中 \(maxval_i\) 表示第 \(i\) 项最大能变成的值,\(minval_i\) 表示第 \(i\) 项最小能变成的值。

算上时间顺序,这是一个三维偏序问题,用 CDQ 分治 + 数据结构(我用了树状数组)就能解决。

1 #include 
2 #include
3 using namespace std; 4 5 const int MN = 100005; 6 const int MC = 100000; 7 8 int N, M; 9 int A[MN], Mx[MN], Mn[MN];10 int f[MN], Ans;11 int p[MN];12 inline bool cmp1(int i, int j) { return Mx[i] < Mx[j]; }13 inline bool cmp2(int i, int j) { return A[i] < A[j]; }14 15 int B[MN];16 inline void Ins(int i, int x) { for (; i <= MC; i += i & -i) B[i] = max(B[i], x); }17 inline void Clr(int i) { for (; i <= MC; i += i & -i) B[i] = 0; }18 inline int Qur(int i) { int A = 0; for (; i; i -= i & -i) A = max(A, B[i]); return A;}19 20 void CDQ(int lb, int rb) {21 if (lb == rb) {22 f[lb] = max(f[lb], 1);23 return;24 }25 int mid = lb + rb >> 1;26 CDQ(lb, mid);27 for (int i = lb; i <= rb; ++i)28 p[i] = i;29 sort(p + lb, p + mid + 1, cmp1);30 sort(p + mid + 1, p + rb + 1, cmp2);31 int j = lb;32 for (int i = mid + 1; i <= rb; ++i) {33 while (j <= mid && Mx[p[j]] <= A[p[i]]) {34 Ins(A[p[j]], f[p[j]]);35 ++j;36 }37 f[p[i]] = max(f[p[i]], Qur(Mn[p[i]]) + 1);38 }39 for (int i = lb; i <= mid; ++i)40 Clr(A[i]);41 CDQ(mid + 1, rb);42 }43 44 int main() {45 int x, y;46 scanf("%d%d", &N, &M);47 for (int i = 1; i <= N; ++i)48 scanf("%d", &A[i]),49 Mx[i] = Mn[i] = A[i];50 for (int i = 1; i <= M; ++i)51 scanf("%d%d", &x, &y),52 Mx[x] = max(Mx[x], y),53 Mn[x] = min(Mn[x], y);54 CDQ(1, N);55 for (int i = 1; i <= N; ++i)56 Ans = max(Ans, f[i]);57 printf("%d\n", Ans);58 return 0;59 }

 

转载于:https://www.cnblogs.com/PinkRabbit/p/10096763.html

你可能感兴趣的文章
ArcGIS REST 缓存清除(地图空白不显示的问题 )
查看>>
第0次作业
查看>>
"类" 库添加继承
查看>>
ucos在s3c2410上运行过程整体剖析之基础知识-与UCOS运行有关的ARM9芯片知识--续 ...
查看>>
存储器的寻址问题 分类: 计算机组成原理 2011-...
查看>>
DDRmenu(翻译)
查看>>
atitit.文件上传带进度条的实现原理and组件选型and最佳实践总结O7
查看>>
Atitit 架构的原则attilax总结
查看>>
ASP.Net之数据绑定
查看>>
Android自动化测试第三季第二讲Toast控件文字获取
查看>>
Google Analytics的能与不能
查看>>
Ubuntu 基本操作
查看>>
JAVA数组的定义及用法
查看>>
18寒假第七测
查看>>
帧中继
查看>>
105:MyBatis常见实用面试题整理
查看>>
Base on QC Automation Framework v1.0
查看>>
bzoj 3261: 最大异或和 (可持久化trie树)
查看>>
UVA 11440 Help Tomisu
查看>>
bzoj千题计划258:bzoj3123: [Sdoi2013]森林
查看>>