题目描述

虽不能至,心向往之。

  $Treap=Tree+Heap$
  椎$=$树$+$堆
  小$\pi$学习了计算机科学中的数据结构$Treap$。
  小$\pi$知道$Treap$指的是一种树。
  小$\pi$还知道$Treap$节点上有两个权值$k$和$w$,其中$k$满足二叉搜索树性质。$w$满足堆性质。
  小$\pi$还知道在$k$和$w$都各不相同的时候,$Treap$的形态是固定的。
  但是小$\pi$不知道这道题目的做法。
  这道题目要求你维护一个大根堆$Treap$,要求支持$n$个操作:
  $0\ k\ w$:插入一个关键字为$k$,权值为$w$的点。
  $1\ k$:删除一个关键字为$k$的点。
  $2\ ku\ kv$:返回关键字分别为$k_u$和$k_v$两个节点的距离
  保证任意时刻树中结点$key$和$weight$都是两两不同的。不会删除当前$Treap$中不存在的点
  小$\pi$找到了学信息学竞赛的你


输入格式

  输入的第一行包含一个正整数$n$,表示操作数。
  下面$n$行,为上述操作之一。


输出格式

  对于每个$2$操作,返回对应答案。


样例

样例输入:

12
0 4 17
0 10 38
0 2 21
0 1 61
2 1 10
2 10 4
1 10
1 1
0 8 42
2 8 4
2 8 2
2 4 2

样例输出:

1
2
2
1
1


数据范围与提示

  对于$20\%$的数据,$n\leqslant 100$。
  对于$40\%$的数据,$n\leqslant 1,000$。
  对于$60\%$的数据,$n\leqslant 10,000$。
  对于$100\%$的数据,$n\leqslant 200,000,0\leqslant k,w,k_u,k_v\leqslant 2^{31}$。


题解

一句话题意:求$Treap$上两点的距离。

$Treap$毕竟还是一棵树,考虑对一棵树如何求两个点$(u,v)$之间的距离,无非就是$depth[u]+depth[v]-2\times depth[lca]$。

那么对于这道题,就是要想办法求出$lca$和$depth$就好了。

因为$key$满足二叉搜索树性质,也就是中序遍历从小到大;那么不妨将其拍扁放在序列上,也就是说维护一个$key$从小到大的序列。

注意空间问题,因为这道题$key$值只与相对大小有关,而与具体值无关,所以可以离散化,并在序列中提前留好坑即可。

那么来考虑$lca$是什么。

画个图便可以知道,$k_u$和$k_v$的$lca$就是在这个序列上$k_u$和$k_v$之间的$weight$值的最大值,用线段树维护即可。

接着考虑$depth$。

对于点$u$,其在序列上左(右)边第一个$weight$比它大的一定是它的祖先;而一个点的$depth=$其祖先的个数,对于这道题也就是$u$左边和右边的上升序列长度之和(注意这里的上升序列是指碰见比它大的就选上,而不是求最长上升子序列)。

考虑平常维护上升序列可以用单调栈,这道题可以用线段树维护单调栈(然而考场并不会,花了一个小时$zzyy$了一下)。

知道了$lca$和$depth$,这道题也就解决了。

时间复杂度:$\Theta(n\log^2 n)$。

期望得分:$100$分。

实际得分:$100$分。


代码时刻

#include<bits/stdc++.h>
#define L(x) x<<1
#define R(x) x<<1|1
using namespace std;
struct rec{int opt,k,w;}e[300000];
unordered_map<int,int>mp;
int n;
int a[300000];
int trmx[1000000],trl[1000000],trr[1000000],trns[1000000],pos[1000000],cnt;
int pushupl(int x,int l,int r,int w)
{
if(l==r)return trmx[x]>w;int mid=(l+r)>>1;
if(trmx[L(x)]<=w)return pushupl(R(x),mid+1,r,w);
return trl[x]-trl[L(x)]+pushupl(L(x),l,mid,w);
}
int pushupr(int x,int l,int r,int w)
{
if(l==r)return trmx[x]>w;int mid=(l+r)>>1;
if(trmx[R(x)]<=w)return pushupr(L(x),l,mid,w);
return trr[x]-trr[R(x)]+pushupr(R(x),mid+1,r,w);
}
void pushup(int x,int l,int r)
{
if(trmx[L(x)]<=trmx[R(x)]){trns[x]=trns[R(x)];trmx[x]=trmx[R(x)];}
else{trns[x]=trns[L(x)];trmx[x]=trmx[L(x)];}int mid=(l+r)>>1;
trl[x]=trl[L(x)]+pushupl(R(x),mid+1,r,trmx[L(x)]);
trr[x]=trr[R(x)]+pushupr(L(x),l,mid,trmx[R(x)]);
}
void add(int x,int l,int r,int k,int w)
{
if(l==r){trns[x]=k;trmx[x]=w;trl[x]=trr[x]=1;return;}
int mid=(l+r)>>1;
if(k<=mid)add(L(x),l,mid,k,w);
else add(R(x),mid+1,r,k,w);
pushup(x,l,r);
}
void del(int x,int l,int r,int k)
{
if(l==r){trns[x]=k;trmx[x]=-0x3f3f3f3f;trl[x]=trr[x]=0;return;}
int mid=(l+r)>>1;
if(k<=mid)del(L(x),l,mid,k);
else del(R(x),mid+1,r,k);
pushup(x,l,r);
}
pair<int,int> lca(int x,int l,int r,int res1,int res2)
{
if(res1<=l&&r<=res2)return make_pair(trmx[x],trns[x]);
int mid=(l+r)>>1;
if(res2<=mid)return lca(L(x),l,mid,res1,res2);
if(mid<res1)return lca(R(x),mid+1,r,res1,res2);
pair<int,int> flag1=lca(L(x),l,mid,res1,res2);
pair<int,int> flag2=lca(R(x),mid+1,r,res1,res2);
if(flag1<flag2)return flag2;return flag1;
}
pair<int,int> askl(int x,int l,int r,int k,int w)
{
if(!k)return make_pair(0,0);
if(r<=k)return make_pair(max(w,trmx[x]),pushupr(x,l,r,w));
int mid=(l+r)>>1;
if(k<=mid)return askl(L(x),l,mid,k,w);
pair<int,int> flag=askl(R(x),mid+1,r,k,w);
pair<int,int> res=askl(L(x),l,mid,k,max(w,flag.first));
res.second+=flag.second;
return res;
}
pair<int,int> askr(int x,int l,int r,int k,int w)
{
if(l>cnt)return make_pair(0,0);
if(k<=l)return make_pair(max(w,trmx[x]),pushupl(x,l,r,w));
int mid=(l+r)>>1;
if(mid<k)return askr(R(x),mid+1,r,k,w);
pair<int,int> flag=askr(L(x),l,mid,k,w);
pair<int,int> res=askr(R(x),mid+1,r,k,max(w,flag.first));
res.second+=flag.second;
return res;
}
int getdep(int x){return askl(1,1,cnt,x-1,pos[x]).second+askr(1,1,cnt,x+1,pos[x]).second;}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%d",&e[i].opt);
switch(e[i].opt)
{
case 0:scanf("%d%d",&e[i].k,&e[i].w);a[++a[0]]=e[i].k;break;
case 1:scanf("%d",&e[i].k);break;
case 2:scanf("%d%d",&e[i].k,&e[i].w);if(e[i].k>e[i].w)swap(e[i].k,e[i].w);break;
}
}
sort(a+1,a+a[0]+1);
for(int i=1;i<=a[0];i++){if(a[i]!=a[i-1])cnt++;mp[a[i]]=cnt;}
memset(trmx,-0x3f,sizeof(trmx));
for(int i=1;i<=n;i++)
switch(e[i].opt)
{
case 0:add(1,1,cnt,mp[e[i].k],e[i].w);pos[mp[e[i].k]]=e[i].w;break;
case 1:del(1,1,cnt,mp[e[i].k]);pos[mp[e[i].k]]=-0x3f3f3f3f;break;
case 2:printf("%d\n",getdep(mp[e[i].k])+getdep(mp[e[i].w])-2*getdep(lca(1,1,cnt,mp[e[i].k],mp[e[i].w]).second));break;
}
return 0;
}

rp++

[CSP-S模拟测试]:椎(线段树维护区间最值和单调栈)的更多相关文章

  1. 滑动窗口(poj,线段树维护区间最值)

    题目描述 现在有一堆数字共N个数字(N<=10^6),以及一个大小为k的窗口.现在这个从左边开始向右滑动,每次滑动一个单位,求出每次滑动后窗口中的最大值和最小值. 例如: The array i ...

  2. 【8.26校内测试】【重构树求直径】【BFS模拟】【线段树维护DP】

    题目性质比较显然,相同颜色联通块可以合并成一个点,重新建树后,发现相邻两个点的颜色一定是不一样的. 然后发现,对于一条链来说,每次把一个点反色,实际上使点数少了2个.如下图 而如果一条链上面有分支,也 ...

  3. POJ.2763 Housewife Wind ( 边权树链剖分 线段树维护区间和 )

    POJ.2763 Housewife Wind ( 边权树链剖分 线段树维护区间和 ) 题意分析 给出n个点,m个询问,和当前位置pos. 先给出n-1条边,u->v以及边权w. 然后有m个询问 ...

  4. 【uoj#164】[清华集训2015]V 线段树维护历史最值

    题目描述 给你一个长度为 $n$ 的序列,支持五种操作: $1\ l\ r\ x$ :将 $[l,r]$ 内的数加上 $x$ :$2\ l\ r\ x$ :将 $[l,r]$ 内的数减去 $x$ ,并 ...

  5. 【bzoj3064】Tyvj 1518 CPU监控 线段树维护历史最值

    题目描述 给你一个序列,支持4种操作:1.查询区间最大值:2.查询区间历史最大值:3.区间加:4.区间赋值. 输入 第一行一个正整数T,表示Bob需要监视CPU的总时间. 然后第二行给出T个数表示在你 ...

  6. CF213E Two Permutations 线段树维护哈希值

    当初竟然看成子串了$qwq$,不过老师的$ppt$也错了$qwq$ 由于子序列一定是的排列,所以考虑插入$1$到$m$到$n-m+1$到$n$; 如何判断呢?可以用哈希$qwq$: 我们用线段树维护哈 ...

  7. Can you answer these queries V SPOJ - GSS5 (分类讨论+线段树维护区间最大子段和)

    recursion有一个整数序列a[n].现在recursion有m次询问,每次她想知道Max { A[i]+A[i+1]+...+A[j] ; x1 <= i <= y1 , x2 &l ...

  8. 线段树维护区间前k小

    线段树维护区间前k小 $ solution: $ 觉得超级钢琴太麻烦?在这里线段树提供一条龙服务 . 咳咳,开始讲正题!这道题我们有一个和超级钢琴复杂度一样 $ ~O(~\sum x\times lo ...

  9. CodeForces - 587E[线段树+线性基+差分] ->(线段树维护区间合并线性基)

    题意:给你一个数组,有两种操作,一种区间xor一个值,一个是查询区间xor的结果的种类数 做法一:对于一个给定的区间,我们可以通过求解线性基的方式求出结果的种类数,而现在只不过将其放在线树上维护区间线 ...

随机推荐

  1. TableLayoutPanel

    1.简单属性 百度经验:https://jingyan.baidu.com/article/9113f81b7966df2b3314c775.html Name属性:tableLayoutPanel经 ...

  2. c++primer chapter one

    一个函数的定义包含四个部分:返回类型(return type),函数名(function name),一个括号包含的形参列表(parameter,允许为空)以及函数体(function body). ...

  3. C++中的单例类模板

    1,本节课讲述单例类模式,实现并抽取相关代码实现单例类模板,在以后开发工作 中,如果想要使用单例模式,那么直接使用今天开发的单例类模板就可以: 2,需求的提出: 1,在架构设计时,某些类在整个系统生命 ...

  4. c++实现直接插入排序

    基本概念 直接插入排序是一种最简单的排序方法,排序过程为:先将第一个元素看作是只有一个元素的有序子表,然后从第二个元素开始,将待排序元素依次插入到前面有序的子表中,直到全部排序完毕.在整个过程中,前面 ...

  5. Ubuntu下软件打开时状态图标与原始图标不重合问题解决

    问题描述 如下图书所示,pycharm打开时,运行的pycharm图标与原始的在收藏夹中的图标不重合.而其他软件不会这样,运行软件的图标与原始图标会重合, 解决办法 把鼠标悬浮在打开的pycharm上 ...

  6. 【BZOJ-4289】Tax 最短路 + 技巧建图(化边为点)

    题意 给出一个N个点M条边的无向图,经过一个点的代价是进入和离开这个点的两条边的边权的较大值,求从起点1到点N的最小代价.起点的代价是离开起点的边的边权,终点的代价是进入终点的边的边权N<=10 ...

  7. IDEA一些有用的功能

    使用 Type Info 如果你想要更多的关于符号的信息,例如从哪里或它的类型是什么, 快速文档可以很好的帮到您,您可以按下 Ctrl+Q 来调用它,然后你会看到一个包含这些细节的弹出窗口.如果您不需 ...

  8. 查看ftp创建的用户

    先su然后查看cat /etc/passwd,500以后的是添加的使用者账户.

  9. MySql 5.7关键字和保留字-附表

    现在使用navicat图形界面或者Hibernate做映射生成表的时候,渐渐的会忽视掉关键字这个问题,而后续也会不断的产生错误提示,一遍遍的查询代码无果,甚至开始怀疑人生,但是其实很多情况下只是使用了 ...

  10. 关于mysql的查询优化

    由于工作原因,最近甲方客户那边多次反应了他们那边的系统查询速度慢,经过排除之后,发现他们那边的数据库完全没有用到索引,简直坑得一笔,通过慢查询日志分析,为数据表建立了适当的索引之后,查询速度明显的提高 ...