题目链接

题意分析

我们假设每一个数都有一个变动范围\([L_i,R_i]\)

那么我们令\(dp[i]\)表示以\(i\)结尾的最长不下降子序列的长度

那么就是\(dp[i]=max\{dp[j]+1\}\)

转移的条件是

\(1.j<i\)

\(2.r_j≤a_i\)

\(3.a_j≤l_i\)

所以这就是一个三位偏序问题

同时需要注意的是

我们当前对于\([l,r]\)区间

我们由于必须要遵循从左向右更新的性质

所以先处理好\([l,mid]\)有\([l,mid]\)向\((mid,r]\)转移

然后再处理\((mid,r]\)

至于具体的处理

我们对于\([l,mid]\)需要按照\(r_i\)从小到大排序

对于\((mid,r]\)需要按照\(a_i\)从小到大排序

然后最后一维使用树状数组处理即可

CODE:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<string>
#include<queue>
#include<map>
#include<stack>
#include<list>
#include<set>
#include<deque>
#include<vector>
#include<ctime>
#define ll long long
#define inf 0x7fffffff
#define N 500008
#define IL inline
#define M 100011
#define D double
#define maxn 100002
#define R register
using namespace std;
template<typename T>IL void read(T &_)
{
T __=0,___=1;char ____=getchar();
while(!isdigit(____)) {if(____=='-') ___=0;____=getchar();}
while(isdigit(____)) {__=(__<<1)+(__<<3)+____-'0';____=getchar();}
_=___ ? __:-__;
}
/*-------------OI使我快乐-------------*/
int n,m,ans;
int tre[M],dp[M];
struct Node{
int id,lx,rx,num;
}e[M];
IL bool cmp1(const Node &A,const Node &B)
{return A.rx<B.rx;}
IL bool cmp2(const Node &A,const Node &B)
{return A.num<B.num;}
IL bool cmp3(const Node &A,const Node &B)
{return A.id<B.id;}
IL void add(int x,int y)
{for(;x<=maxn;x+=x&-x) tre[x]=max(tre[x],y);}
IL int qury(int x)
{int res=0;for(;x;x-=x&-x) res=max(res,tre[x]);return res;}
IL void clean(int x)
{for(;x<=maxn;x+=x&-x) tre[x]=0;}
IL void solve(int le,int ri)
{
if(le==ri) return;
int mid=(le+ri)>>1;
solve(le,mid);
sort(e+le,e+mid+1,cmp1);sort(e+mid+1,e+ri+1,cmp2);
int cdy=le,wzy=mid+1;
for(;wzy<=ri;++wzy)
{
for(;cdy<=mid&&e[cdy].rx<=e[wzy].num;++cdy)
add(e[cdy].num,dp[e[cdy].id]);
dp[e[wzy].id]=max(dp[e[wzy].id],qury(e[wzy].lx)+1);
}
for(R int i=le;i<cdy;++i) clean(e[i].num);
sort(e+mid+1,e+ri+1,cmp3);
solve(mid+1,ri); }
int main()
{
// freopen(".in","r",stdin);
// freopen(".out","w",stdout);
read(n);read(m);
for(R int i=1;i<=n;++i)
{
read(e[i].num);e[i].id=i;;
e[i].lx=e[i].rx=e[i].num;
}
for(R int i=1,x,y;i<=m;++i)
{
read(x);read(y);
e[x].lx=min(e[x].lx,y);
e[x].rx=max(e[x].rx,y);
}
// for(R int i=1;i<=n;++i)
// printf("%d %d %d %d\n",e[i].num,e[i].lx,e[i].rx,e[i].id);
for(R int i=1;i<=n;++i) dp[i]=1;
solve(1,n);
for(R int i=1;i<=n;++i) ans=max(ans,dp[i]);
printf("%d\n",ans);
// fclose(stdin);
// fclose(stdout);
return 0;
}

HEOI 2019 RP++

P4093 [HEOI2016/TJOI2016]序列的更多相关文章

  1. 洛谷 P4093 [HEOI2016/TJOI2016]序列 解题报告

    P4093 [HEOI2016/TJOI2016]序列 题目描述 佳媛姐姐过生日的时候,她的小伙伴从某宝上买了一个有趣的玩具送给他.玩具上有一个数列,数列中某些项的值可能会变化,但同一个时刻最多只有一 ...

  2. 洛谷 P4093 [HEOI2016/TJOI2016]序列 CDQ分治优化DP

    洛谷 P4093 [HEOI2016/TJOI2016]序列 CDQ分治优化DP 题目描述 佳媛姐姐过生日的时候,她的小伙伴从某宝上买了一个有趣的玩具送给他. 玩具上有一个数列,数列中某些项的值可能会 ...

  3. Luogu P4093 [HEOI2016/TJOI2016]序列 dp套CDQ

    题面 好久没写博客了..最近新学了CDQ...于是就来发一发一道CDQ的练习题 看上去就是可以dp的样子. 设\(dp_{i}\)为以i结尾的最长不下降序列. 易得:\(dp_{i}\)=\(max( ...

  4. BZOJ4553/洛谷P4093 [HEOI2016/TJOI2016]序列 动态规划 分治

    原文链接http://www.cnblogs.com/zhouzhendong/p/8672434.html 题目传送门 - BZOJ4553 题目传送门 - 洛谷P4093 题解 设$Li$表示第$ ...

  5. 洛谷P4093 [HEOI2016/TJOI2016]序列

    题目描述 佳媛姐姐过生日的时候,她的小伙伴从某宝上买了一个有趣的玩具送给他.玩具上有一个数列,数列中某些项的值可能会变化,但同一个时刻最多只有一个值发生变化.现在佳媛姐姐已经研究出了所有变化的可能性, ...

  6. 洛谷 P4093 [HEOI2016/TJOI2016]序列(Cdq+dp)

    题面 luogu 题解 \(Cdq分治+dp\) \(mx[i],mn[i]\)分别表示第\(i\)位最大,最小能取到多少 那么有 \(j < i\) \(mx[j] \le a[i]\) \( ...

  7. cdq分治(hdu 5618 Jam's problem again[陌上花开]、CQOI 2011 动态逆序对、hdu 4742 Pinball Game、hdu 4456 Crowd、[HEOI2016/TJOI2016]序列、[NOI2007]货币兑换 )

    hdu 5618 Jam's problem again #include <bits/stdc++.h> #define MAXN 100010 using namespace std; ...

  8. 洛谷 P4093: bzoj 4553: [HEOI2016/TJOI2016]序列

    题目传送门:洛谷P4093. 题意简述: 给定一个长度为 \(n\) 的序列 \(a\). 同时这个序列还可能发生变化,每一种变化 \((x_i,y_i)\) 对应着 \(a_{x_i}\) 可能变成 ...

  9. 题解 P4093 【[HEOI2016/TJOI2016]序列】

    这道题原来很水的? noteskey 一开始以为是顺序的 m 个修改,然后选出一段最长子序列使得每次修改后都满足不降 这 TM 根本不可做啊! 于是就去看题解了,然后看到转移要满足的条件的我发出了黑人 ...

随机推荐

  1. Android 重写EditText回车事件

    之前遇到的问题没来得及记录下来,趁今晚有空就重新回忆并写下了. 我们在用到EditText这个空间时经常需要重写软键盘中的回车事件以配合我们接下来的响应,比如点击回车变成搜索.发送.完成等. Edit ...

  2. 智能合约调用另一合约中的payable方法

    参考链接: https://ethereum.stackexchange.com/questions/9705/how-can-you-call-a-payable-function-in-anoth ...

  3. Nmap扫描参数解析

    Nmap是一款开源免费的网络发现(Network Discovery)和安全审计(Security Auditing)工具.软件名字Nmap是Network Mapper的简称.Nmap最初是由Fyo ...

  4. 浅谈 js中parseInt函数的解析[转]

    首先还是从很热门的实例parseInt("09")==0说起. parseInt(number,type)这个函数后面如果不跟第2个参数来表示进制的话,默认是10进制. 比如说pa ...

  5. [Training Video - 4] [Groovy] Initializing log inside class with constructor

    TestService s = new TestService(log,context,testRunner) s.xyz() class TestService{ def log def conte ...

  6. docker镜像存出与载入

    尝试从官网上下载ubuntu镜像,太慢下载不下来. 使用daocloud加速器进行加速之后,由于公司网络不好,仍然下载不下来. 没办法,只能从别的环境上搞一个已经存在的ubuntu镜像,折腾到自己的虚 ...

  7. ajax中 XmlHttp的open( )方法

    博客分类:  Ajax XML  open 创建一个新的http请求,并指定此请求的方法.URL以及验证信息 语法 oXMLHttpRequest.open(bstrMethod, bstrUrl, ...

  8. 3.1.7 线程阻塞工具类:LockSupport

    package 第三章.线程阻塞工具LockSupport; import java.util.concurrent.locks.LockSupport; /** * Created by zzq o ...

  9. SIP简介

    说明:以下内容来着之前下载的一份文档,现将概念部分摘录在BLog,如需要完整文档将放在文件中或留言. SIP简介,第1部分:SIP初探 时间:2006-04-07作者:Emmanuel Proulx浏 ...

  10. 从零开始学习前端JAVASCRIPT — 11、JavaScript运动模型及轮播图效果、放大镜效果、自适应瀑布流

    未完待续...... 一.运动原理 通过连续不断的改变物体的位置,而发生移动变化. 使用setInterval实现. 匀速运动:速度值一直保持不变. 多物体同时运动:将定时器绑设置为对象的一个属性. ...