题意

给你长度为$n$的序列,序列中的每个元素$i$有一个区间限制$[l_i,r_i]$,你从中选出一个子序列,并给它们标号$x_i$,要求满足 $,∀i<j,x_i<x_j$,且$, ∀i,x_i∈[l_i,r_i]$。 问满足条件子序列的长度最长为多少?

其中$1\leq n\leq3\times 10^5\ 1\leq l_i\leq r_i\leq 10^9$

题解

不妨设$f[i][j]$表示已经选到第$i$个,其中最大值为$j$最多能选几个。

显然是开不下的...但是还记得$O(nlogn)$的$LIS$吗?它利用了二分栈!

所以我们不妨考虑设$f[i][j]$表示当前选到第$i$个,已经选了$j$个的末尾最小元素。

如果不选,显然$f[i+1][j]=f[i][j]$

如果$r[i+1]>f[i][j]$,那么,则有$f[i+1][j+1]=max(f[i][j]+1,l[i+1])$

不难发现可以用滚动数组滚掉第一维,所以第一个方程直接作废:

$ f[j+1]=max(f[j]+1,l[i+1]) $

但是时间还是$O(n^2)$的啊,枚举一个$i$,枚举一个$j$。

不急,咱们分开考虑,对于一个$l[i]<f<r[i]$

它的当前$dp$值可以直接$+1$,那么位置也将$+1$因为多选了一个嘛。

剩下的话,就是$f<l[i]$的情况,直接就等于$l[i]$了。

用一颗$FHQ-Treap$维护一下就行了。

$PS$:注意数组开两倍(最开始$n$个以及$dp$中的$n$个新节点)

#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm> template<typename T>
void read(T &x) {
int flag = 1; x = 0; char ch = getchar();
while(ch < '0' || ch > '9') { if(ch == '-') flag = -flag; ch = getchar(); }
while(ch >= '0' && ch <= '9') x = x * 10 + ch - '0', ch = getchar(); x *= flag;
} const int N = 6e5 + 10, Inf = 1e9 + 7;
int n, l, r;
int rt, tot, val[N], add[N], lc[N], rc[N], pri[N]; inline int node(int x) { val[++tot] = x, pri[tot] = rand(); return tot; }
inline void pushdown(int o) {
if(add[o]) {
val[o] += add[o];
if(lc[o]) add[lc[o]] += add[o];
if(rc[o]) add[rc[o]] += add[o];
add[o] = 0;
}
}
int merge(int l, int r) {
if(!l || !r) return l + r;
pushdown(l), pushdown(r);
if(pri[l] < pri[r]) { rc[l] = merge(rc[l], r); return l; }
else { lc[r] = merge(l, lc[r]); return r; }
}
void split(int o, int k, int &l, int &r) {
if(!o) { l = r = 0; return ; } pushdown(o);
if(val[o] <= k) l = o, split(rc[o], k, rc[o], r);
else r = o, split(lc[o], k, l, lc[o]);
}
int min(int o) { pushdown(o); return lc[o] ? min(lc[o]) : val[o]; }
int calc(int o) {
if(!o) return 0; pushdown(o);
return calc(lc[o]) + calc(rc[o]) + (val[o] < Inf);
} int main () {
read(n), srand(19260817); int x, y, k, p;
for(int i = 1; i <= n; ++i) rt = merge(rt, node(Inf));
for(int i = 1; i <= n; ++i) {
read(l), read(r);
split(rt, l - 1, x, y), split(y, r - 1, k, y);
if(k) ++add[k];
split(y, min(y), p, y);
rt = merge(merge(merge(x, node(l)), k), y);
} printf("%d\n", calc(rt));
return 0;
}

最新文章

  1. jQuery EasyUI 使用笔记
  2. 北航 编译实践 PL/0文法
  3. Integer封装与拆箱
  4. 黄聪:走进wordpress do_action函数
  5. css div 垂直居中
  6. Vijos_1792_摆花_(动态规划,多重集组合数)
  7. [译]内存中的Java数组是怎么样的
  8. 通过HttpClient 调用ASP.NET Web API
  9. Java模拟http请求调用远程接口工具类
  10. jenkins grunt 自动构建流程
  11. python 基础 列表
  12. python中class的序列化和反序列化
  13. HTML前期学习总结
  14. php数据类型之自动转换和强制转换
  15. Echarts折线图点击事件
  16. SElinux测试及排错
  17. Spring IOC-ContextLoaderListener
  18. SEO之H1,H2,H3,H4....STRONG使用方法
  19. 为 Tomcat 安装 apr
  20. asp.net mvc部署iis常见问题

热门文章

  1. java 7修改文件权限
  2. java提取SVN提交log
  3. 暑假集训——cf热身赛部分题有感加其题解
  4. fundamentals of the jQuery library
  5. docker 环境
  6. 使用Sysmon分析宏病毒(Macros Downloader)
  7. python中BeautifulSoup模块
  8. python基础===Sublime Text 3 快捷键
  9. python基础===discover函数介绍
  10. js-xlsx操作excel表格