【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 第一行一 ...
随机推荐
- Oracle 监听/数据库 启动/关闭
LSNRCTL命令启动.关闭和查看监听器的状态的方法 从lsnrctl status命令的输出中得到监听器状态,包括如下的信息: 监听器的启动时间 监听器的运行时间 监听器参数文件listener.o ...
- 爬虫:Scrapy14 - Telnet 终端(Telnet Console)
Scrapy 提供了内置的 Telnet 终端,以供检查,控制 Scrapy 运行的进程.Telnet 仅仅是一个运行在 Scrapy 进程中的普通 Python 终端.因此你可以在其中做任何是. T ...
- TZOJ3043: 取个标题好难 最长的出现次数>=k的不重复子串长度
3043: 取个标题好难 Time Limit(Common/Java):6000MS/18000MS Memory Limit:65536KByteTotal Submit: 17 ...
- 【python基础】--常用数据结构
list tuple dict set四种常用数据结构 list list 有序的集合,可以随时添加.删除其中元素值; 支持list嵌套模式, >>> p = ['a','b']&g ...
- J2EE的十三种技术——JNDI
背景: 上一篇博客中介绍了J2EE的十三种技术之一--JDBC,主要用于提供了统一访问多种数据库的方式.这篇文章我们继续介绍J2EE的技术--JNDI. JNDI: Java Naming and D ...
- Java实现身份证号码校验
二话不说,直接上代码. package hope.identitycodecheck.demo; import java.text.DateFormat; import java.text.Simpl ...
- nodejs & docker
nodejs & docker https://github.com/xgqfrms-GitHub/Node.js/blob/master/Docker-Nodejs/translation. ...
- 【bzoj5015】[Snoi2017]礼物 矩阵乘法
题目描述 热情好客的请森林中的朋友们吃饭,他的朋友被编号为 1-N,每个到来的朋友都会带给他一些礼物:.其中,第一个朋友会带给他 1 个,之后,每一个朋友到来以后,都会带给他之前所有人带来的礼物个数再 ...
- 基于MPLAB X IDE配置位设置讲解
http://blog.csdn.net/superanters/article/details/8541171 在讲基于MPLAB X IDE 配置位配置前我先讲讲如何配置配置位. 比如PICLF1 ...
- webpack watch模式产生*.hot-update.json文件
webpack --watch会产生*.hot-update.json文件,解决方法如下: output: { path: path.join(root, "dist"), fil ...