【bzoj3173】【Tjoi2013】【最长上升子序列】treap+dp二分优化
[pixiv] https://www.pixiv.net/member_illust.php?mode=medium&illust_id=61560361
向大(hei)佬(e)实力学(di)习(tou)
Description
给定一个序列,初始为空。现在我们将1到N的数字插入到序列中,每次将一个数字插入到一个特定的位置。每插入一个数字,我们都想知道此时最长上升子序列长度是多少?
Input第一行一个整数N,表示我们要将1到N插入序列中,接下是N个数字,第k个数字Xk,表示我们将k插入到位置Xk(0<=Xk<=k-1,1<=k<=N)
OutputN行,第i行表示i插入Xi位置后序列的最长上升子序列的长度是多少。
Sample Input3
0 0 2
Sample Output1
1
2
HINT100%的数据 n<=100000
其实这道题算是一道裸题了,拿来练习一下模板。看到这种活动的插入,就能立即想到平衡树。至于最长上升子序列,如果每插入一次就求一次,再优化都要超时。于是我们就找找性质,发现:由于是有序插入,后插入的一定大,那么每插入一个点,以此点结尾的最长上升子序列就不受后面点的干扰。那么在插入完毕后求一次最长上升子序列就可以了。之后的输出需注意要求的是整个数列的答案,还要在小于等于插入数的dp数组中求max。因为这个我又wa了一遍
至于二分优化,其实也算是练手。简单梳理一下,也便自己理解更透彻。要开一个单调栈(?),存两个元素,一个为值,一个为下标。
以 6 2 1 4 3 为例,从左往右扫
6进栈 >6(1)
dp[1]=1;
2比栈顶元素小,要去替换,二分到6,替换>2(2)
dp[2]=1;
1同理>1(3)
dp[3]=1;
4>1,4直接进栈>1(3) 4(4)
dp[4]=dp[3]+1=2;
3二分换4,>1(3) 3(5)
dp[5]=dp[3]+1=2;
最后再for一遍找max,总o(log n)
这道题因为比较特殊,数字不重复,且按序插入,所以直接以数字本身作为下标。
求最长上升子序列还可以用树状数组优化,哪天再遇到这样的题再复习树状数组优化吧!
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std ;
const int N=100000+5;
int n,x;
struct Node{
Node *ch[2];
int v,r,siz;
int cmp(int x){
return x <= ch[0]->siz ? 0 : 1;
}
}*root,*null,pool[N],*tail=pool;
int c[N],top=0,dp[N];
Node *newnode(){
Node *rt=++tail;
rt->ch[0]=rt->ch[1]=null;
rt->v=rt->r=rt->siz=0;
return rt;
}
void rotate(Node *&nd,int d){
Node *tp=nd->ch[d];
nd->ch[d]=tp->ch[d^1];tp->ch[d^1]=nd;
nd->siz=nd->ch[0]->siz+nd->ch[1]->siz+1;
tp->siz=tp->ch[0]->siz+tp->ch[1]->siz+1;
nd=tp;
}
void insert(Node *&nd,int val,int pos){
if(nd==null){
nd=newnode();
nd->v=val;
nd->siz=1;
nd->r=rand();
return ;
}
int d=nd->cmp(pos);
insert(nd->ch[d],val,d ? pos - nd->ch[0]->siz - 1: pos );
if(nd->ch[d]->r > nd->r ){
rotate(nd,d);
}
nd->siz=nd->ch[0]->siz+nd->ch[1]->siz+1;
}
int get_dp(int x){
int rt;
if(x>c[top]){
c[++top]=x;
return dp[c[top-1]]+1;
}
int le=0,ri=top;
while(le<ri){
int mid=(le+ri)>>1;
if(c[mid]>x) ri=mid;
if(c[mid]<=x) le=mid+1;
}
c[le]=x;
return dp[c[le-1]]+1;
}
void query(Node *nd){
if(nd==null) return ;
query(nd->ch[0]);
dp[nd->v]=get_dp(nd->v);
query(nd->ch[1]);
}
int main(){
srand(23425546);
null=++tail;
null->ch[0]=null->ch[1]=null;
null->v=null->r=null->siz=0;
root=null;
scanf("%d",&n);
for(int i=1;i<=n;i++){
scanf("%d",&x);
insert(root,i,x);
}
query(root);
for(int i=1;i<=n;i++){
dp[i]=max(dp[i-1],dp[i]);
printf("%d\n",dp[i]);
}
return 0;
}
总结:
1、对题的辨析能力还需加强,要多想多观察,找到不被影响的东西
【bzoj3173】【Tjoi2013】【最长上升子序列】treap+dp二分优化的更多相关文章
- BZOJ3173 TJOI2013最长上升子序列(Treap+ZKW线段树)
传送门 Description 给定一个序列,初始为空.现在我们将1到N的数字插入到序列中,每次将一个数字插入到一个特定的位置.每插入一个数字,我们都想知道此时最长上升子序列长度是多少? Input ...
- [BZOJ3173][Tjoi2013]最长上升子序列
[BZOJ3173][Tjoi2013]最长上升子序列 试题描述 给定一个序列,初始为空.现在我们将1到N的数字插入到序列中,每次将一个数字插入到一个特定的位置.每插入一个数字,我们都想知道此时最长上 ...
- BZOJ 3173: [Tjoi2013]最长上升子序列 [splay DP]
3173: [Tjoi2013]最长上升子序列 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1613 Solved: 839[Submit][St ...
- bzoj3173[Tjoi2013]最长上升子序列 平衡树+lis
3173: [Tjoi2013]最长上升子序列 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2253 Solved: 1136[Submit][S ...
- bzoj3173: [Tjoi2013]最长上升子序列(树状数组+二分倒推)
3173: [Tjoi2013]最长上升子序列 题目:传送门 题解: 好题! 怎么说吧...是应该扇死自己...看错了两次题: 每次加一个数的时候,如果当前位置有数了,是要加到那个数的前面,而不是直 ...
- 【bzoj3173】[Tjoi2013]最长上升子序列 Treap
题目描述 给定一个序列,初始为空.现在我们将1到N的数字插入到序列中,每次将一个数字插入到一个特定的位置.每插入一个数字,我们都想知道此时最长上升子序列长度是多少? 输入 第一行一个整数N,表示我们要 ...
- bzoj千题计划316:bzoj3173: [Tjoi2013]最长上升子序列(二分+树状数组)
https://www.lydsy.com/JudgeOnline/problem.php?id=3173 插入的数是以递增的顺序插入的 这说明如果倒过来考虑,那么从最后一个插入的开始删除,不会对以某 ...
- bzoj3173: [Tjoi2013]最长上升子序列(fhqtreap)
这题用fhqtreap可以在线. fhqtreap上维护以i结尾的最长上升子序列,数字按从小到大加入, 因为前面的数与新加入的数无关, 后面的数比新加入的数小, 所以新加入的数对原序列其他数的值没有影 ...
- BZOJ3173:[TJOI2013]最长上升子序列(Splay)
Description 给定一个序列,初始为空.现在我们将1到N的数字插入到序列中,每次将一个数字插入到一个特定的位置.每插入一个数字,我们都想知道此时最长上升子序列长度是多少? Input 第一行一 ...
随机推荐
- 珍藏版 Python 开发工程师面试试题
珍藏版 Python 开发工程师面试试题 说明:不拿到几家公司的offer,那就是卑鄙的浪费 一.Python_基础语法 1.可变与不可变类型: 2.浅拷贝与深拷贝的实现方式.区别:deepcopy如 ...
- CSU-2110 Keeping Cool
题目链接 http://acm.csu.edu.cn:20080/csuoj/problemset/problem?pid=2110 题目 Description Kevin has just got ...
- vmware设置静态ip(复制)
一.安装好虚拟后在菜单栏选择编辑→ 虚拟网络编辑器,打开虚拟网络编辑器对话框,选择Vmnet8 Net网络连接方式,随意设置子网IP,点击NAT设置页面,查看子网掩码和网关,后面修改静态IP会用到. ...
- JavaScript中的parseInt和Number函数
函数作用: parseInt将字符串(String)类型转为整数类型. Number() 函数把对象(Object)的值转换为数字. 语法不同: parseInt(string, [radix]) s ...
- CSS 的Hack 问题
1.什么是CSS hack? CSS hack是通过在CSS样式中加入一些特殊的符号,让不同的浏览器识别不同的符号(什么样的浏览器识别什么样的符号是有标准的,CSS hack就是让你记住这个标准),以 ...
- html5 download blob
html5 download blob https://stackoverflow.com/questions/19327749/javascript-blob-filename-without-li ...
- 实时流处理Storm、Spark Streaming、Samza、Flink孰优孰劣
对于一个成熟的消息中间件而言,消息格式不仅关系到功能维度的扩展,还牵涉到性能维度的优化.随着Kafka的迅猛发展,其消息格式也在不断的升级改进,从0.8.x版本开始到现在的1.1.x版本,Kafka的 ...
- [poj] 2749 building roads
原题 2-SAT+二分答案! 最小的最大值,这肯定是二分答案.而我们要2-SATcheck是否在该情况下有可行解. 对于目前的答案limit,首先把爱和恨连边,然后我们n^2枚举每两个点通过判断距离来 ...
- BZOJ1296 [SCOI2009]粉刷匠 【dp】
题目 windy有 N 条木板需要被粉刷. 每条木板被分为 M 个格子. 每个格子要被刷成红色或蓝色. windy每次粉刷,只能选择一条木板上一段连续的格子,然后涂上一种颜色. 每个格子最多只能被粉刷 ...
- bzoj1396&&2865 识别子串 后缀自动机+线段树
Input 一行,一个由小写字母组成的字符串S,长度不超过10^5 Output L行,每行一个整数,第i行的数据表示关于S的第i个元素的最短识别子串有多长. Sample Input agoodco ...