当前位置: 首页 > 编程日记 > 正文

【bzoj2770】YY的Treap 权值线段树

题目描述

志向远大的YY小朋友在学完快速排序之后决定学习平衡树,左思右想再加上SY的教唆,YY决定学习Treap。友爱教教父SY如砍瓜切菜般教会了YY小朋友Treap(一种平衡树,通过对每个节点随机分配一个priority,同时保证这棵平衡树关于priority是一个小根堆以保证效率)。这时候不怎么友爱的510跑了出来,他问了YY小朋友一个极不和谐的问题:怎么求Treap中两个点之间的路径长度。YY秒了之后决定把这个问题交给你来做,但只要求出树中两点的LCA。

输入

第一行两个整数n,m

第二行n个整数表示每个元素的key

第三行n个整数表示每个元素的priority

接下m行,每行一条命令

I A B,插入一个元素,key为A, priority为B

D A,删除一个元素,key为A

Q A B,询问key分别为A和B的LCA的key

输出

对于每个Q输出一个整数。

样例输入

2 1
1 2
4 5
Q 1 2

样例输出

1


题解

权值线段树

一个小结论:不妨设$a\le b$,则Treap中权值为$a$、$b$两点的LCA为权值在$[a,b]$之间,优先级最小的点。

证明:

1.LCA的权值在$[a,b]$之间:考虑从$a$、$b$找到LCA的最后一步:一定是$a$在LCA的非右子树,$b$在LCA的非左子树。否则$a$、$b$在同子树内则LCA可以为更优的该儿子节点。

2.权值在$[a,b]$之间的节点一定都在LCA的子树内:如果不在LCA的子树内,那么节点如果在LCA右侧则一定大于$b$,或在LCA左侧则一定小于$a$。

3.LCA的子树中所有节点的优先级都小于等于LCA:证明显然。

4.LCA一定在LCA的子树内:证明显然。

因此由1、2、3、4得证。

于是只需要对每个数的key维护权值线段树,维护权值在某区间内的数的优先级最小值及其位置。查询时直接区间查询即可。

为了避免一些细节问题(比如两个int加起来爆int之类的),代码中使用了离线离散化。

时间复杂度$O(n\log n)$。

#include <cstdio>
#include <utility>
#include <algorithm>
#define N 100010
#define lson l , mid , x << 1
#define rson mid + 1 , r , x << 1 | 1
#define id(x) lower_bound(v + 1 , v + tot + 1 , x) - v
#define inf 0x7fffffff
using namespace std;
typedef pair<int , int> pr;
int key[N] , pri[N] , opt[N * 3] , qx[N * 3] , qy[N * 3] , v[N << 2] , tot;
pr mn[N << 4];
char str[5];
void build(int l , int r , int x)
{mn[x] = pr(inf , inf);if(l == r) return;int mid = (l + r) >> 1;build(lson) , build(rson);
}
void insert(int p , int v , int l , int r , int x)
{if(l == r){mn[x] = pr(v , p);return;}int mid = (l + r) >> 1;if(p <= mid) insert(p , v , lson);else insert(p , v , rson);mn[x] = min(mn[x << 1] , mn[x << 1 | 1]);
}
void erase(int p , int l , int r , int x)
{if(l == r){mn[x] = pr(inf , inf);return;}int mid = (l + r) >> 1;if(p <= mid) erase(p , lson);else erase(p , rson);mn[x] = min(mn[x << 1] , mn[x << 1 | 1]);
}
pr query(int b , int e , int l , int r , int x)
{if(b <= l && r <= e) return mn[x];int mid = (l + r) >> 1;pr ans(inf , inf);if(b <= mid) ans = min(ans , query(b , e , lson));if(e > mid) ans = min(ans , query(b , e , rson));return ans;
}
int main()
{int n , m , i;scanf("%d%d" , &n , &m);for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &key[i]) , v[++tot] = key[i];for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &pri[i]);for(i = 1 ; i <= m ; i ++ ){scanf("%s%d" , str , &qx[i]);if(str[0] == 'I') opt[i] = 1 , scanf("%d" , &qy[i]) , v[++tot] = qx[i];else if(str[0] == 'D') opt[i] = 2;else opt[i] = 3 , scanf("%d" , &qy[i]);}sort(v + 1 , v + tot + 1);build(1 , tot , 1);for(i = 1 ; i <= n ; i ++ ) insert(id(key[i]) , pri[i] , 1 , tot , 1);for(i = 1 ; i <= m ; i ++ ){if(opt[i] == 1) insert(id(qx[i]) , qy[i] , 1 , tot , 1);else if(opt[i] == 2) erase(id(qx[i]) , 1 , tot , 1);else if(qx[i] < qy[i]) printf("%d\n" , v[query(id(qx[i]) , id(qy[i]) , 1 , tot , 1).second]);else printf("%d\n" , v[query(id(qy[i]) , id(qx[i]) , 1 , tot , 1).second]);}return 0;
}

转载于:https://www.cnblogs.com/GXZlegend/p/7570637.html

相关文章:

第二章 数据类型、运算符与表达式

#include <stdio.h> #include <stdlib.h>int Add(int a,int b) {return ab; }int main() {int x,y,sum 0;printf("Input two integers:\n");scanf("%d%d",&x,&y);sum Add(x,y);printf("sum %d\n",sum);return 0; }计算圆…

从设计原则谈软件开发(二)

最近一直在一个培训公司做着极为无聊的培训&#xff0c;所以一直都没有时间上网。今天突然发现这里可以上无线&#xff0c;嘿嘿&#xff0c;就上来继续把这个文章完成。 上次说到了设计原则中的单一职责原则&#xff0c;今天时间比较紧&#xff0c;我就继续往下写&#xff0c;也…

Npm环境依赖重置

需要彻底删除node_modules 然后npm install 如果无法删除node_modules文件&#xff0c;可以试试这个&#xff1a; 安装RimRaf&#xff1a; npm install rimraf -g 并在项目文件夹中删除node_modules文件夹&#xff1a; rimraf node_modules 然后你可以去npm install转载于:http…

MySQL的information_schema

在一次清空一张比较大的表时&#xff08;在清空前占用400多兆&#xff09;&#xff0c;发现该表中记录为0条但是空间并没有被释放&#xff0c;采用下面方式可查看占用情况 -- 查询各个数据库占用磁盘的情况 select TABLE_SCHEMA, concat(truncate(sum(data_length)/1024/1024,2…

编写spring应用

测试类自动注入失败&#xff1a;RunWith(SpringRunner.class)详解 CtrlAltDelete键&#xff0c;打开任务管理器&#xff0c;结束占据8080端口的Tomcat进程。 HomeController.java <!DOCTYPE html> <html xmlns"http://www.w3.org/1999/xhtml"xmlns:th&quo…

[导入]Learning.ASP.NET 2.0.with.AJAX.pdf(14.14 MB)

ASP.NET 2.0的AJAX无疑是最快&#xff0c;最有效&#xff0c;最可靠和最佳的方式支持创建交互式Web应用程序上市。结合开发工具&#xff0c;可以从Microsoft &#xff0c;免费和商业&#xff0c;这是难以置信轻松地创建网站&#xff0c;看看伟大的表现良好。最重要的是&#xf…

c# IO线程 打造 定时打开指定程序

用IO以及线程轻松实现 定时器 &#xff0c;在指定的时间打开指定的程序&#xff1a;&#xff09; 首先是如何实现定时&#xff1f;这可以单独的用个线程&#xff0c;在时间到的时候打开程序 然后是如何打开程序 &#xff0c;用Process.Start就可以了 最后就是如何把程序列表保存…

Python操作 RabbitMQ、Redis、Memcache、SQLAlchemy

Memcached Memcached 是一个高性能的分布式内存对象缓存系统&#xff0c;用于动态Web应用以减轻数据库负载。它通过在内存中缓存数据和对象来减少读取数据库的次数&#xff0c;从而提高动态、数据库驱动网站的速度。Memcached基于一个存储键/值对的hashmap。其守护进程&#xf…

jQuery中的常用内容总结(一)

jQuery中的常用内容总结(一) 前言 不好意思(✿◠‿◠)&#xff0c;由于回家看病以及处理一些其它事情耽搁了&#xff0c;不然这篇博客本该上上周或者上周写的&#xff1b;同时闲谈几句&#xff1a;在这里建议各位开发的童鞋&#xff0c;如果有疾病尽快治疗&#xff0c;不要拖&a…

你需要眼光和资格

机上遇到一男人&#xff0c;操北京口音&#xff0c;三十二三&#xff0c;婚否不详&#xff0c;容貌体面。 优势&#xff1a;技术好&#xff0c;聪明&#xff0c;没坏心&#xff0c;乐观 劣势&#xff1a;有点懒&#xff0c;自傲&#xff0c;责任心与意志力指数一般 其所谓“恰当…

嵌入式系统基础了解

1.一个启动 .s文件&#xff08;start.s&#xff09;&#xff0c;至少需要包含三个段: ;//堆栈段 ; //中断向量表 ; //代码段 参考代码&#xff1a; stack_size EQU 0x200 ;…

BFD与IGP快速收敛应用测试

cqmmx&#xff0c;2008-9-10 1、背景介绍 目前对网络稳定性影响较大的一般是链路中断、节点失效等故障&#xff0c;而常规的慢Hello机制检测耗时较长&#xff0c;且常用IGP&#xff08;ISIS和OSPF&#xff09;在默认配置情况下&#xff0c;收敛速度很慢&#xff0c;一般需要几十…

rsa证书ssh登陆服务器

好久不用&#xff0c;又生疏了。 今晚实操了一下&#xff0c;作一个记录。 使用rsa的密钥对登陆linux服务器&#xff0c;主要是为了安全。 这种证书级别的登陆&#xff0c;比最复杂的root用户名和帐号的安全性都要高一个等级。 至少服务器不会被暴破(暴力破解)。 ~~~~~~~~~~~~~…

简单数据结构(队列 栈 树 堆 )

基础知识 基本概念 程序 算法 数据结构数据结构是计算机存储、组织数据的方式。数据结构是指相互之间存在一种或多种特定关系的数据元素的集合。通常情况下&#xff0c;精心选择的数据结构可以带来更高的运行或者存储效率。数据结构往往同高效的检索算法和索引技术有关。 常见…

cannot access a closed file

用上传控件上传文件时&#xff0c;当执行SaveAs()时出错&#xff0c;异常为cannot access a closed file&#xff1b; 当传小文件没有异常&#xff0c;当超过80k时就出上述异常&#xff0c;后来发现要将Web.config里增加 <!--设置上传附近的大小--> <httpRuntime …

实验四-常用图像增强方法

1、采用二维中值滤波函数medfilt2 对受椒盐噪声干扰的图像滤波&#xff0c;窗口分别采用3*3,5*5,7*7 i imread(D:\study\third_down\ImageProcessing\work\work_one\flower.jpg); I rgb2gray(i); J imnoise(I,salt & pepper,0.04); K1 medfilt2(J,[3 3]);%对矩阵I进行…

商贸通服装鞋帽版客户端无法连接服务器的问题(自己遇到的,已解决)

今天给一客户装“商贸通服装鞋帽版”&#xff0c;客户机是xp&#xff0c;服务器是win2003,装好后在连接服务器是总是报“无法连接服务器”之类的错&#xff0c;但是客户机和服务器双向都可以ping的通&#xff0c;而且双方的防火墙都已关闭&#xff0c;没有第三方防火墙&#xf…

【转】Visual C#创建和使用ActiveX组件

开发基于.net平台上的程序员是很难从本质上把Visual C#和ActiveX组件联起来&#xff0c;虽然在使用Visual C#开发应用程序时&#xff0c;有时为了快速开发或者由于.Net Framework SDK的不完整&#xff0c;还需要借助ActiveX。但即使如此&#xff0c;也很难把二者联系起来。其中…

Why Sleeping May Be More Important Than Studying

Why Sleeping May Be More Important Than Studying转载于:https://www.cnblogs.com/Lamfai/p/10441451.html

测试开发板与主机之间通过串口收发数据(uart.c/uart.h )

usart.c: #include "usart.h"// U1_TX: PA9 // U1_RX: PA10 void usart_init(void) {//1. GPIO口的配置RCC_AHB1PeriphClockCmd(RCC_AHB1Periph_GPIOA,ENABLE); GPIO_InitTypeDef p;p.GPIO_Mode GPIO_Mode_AF;p.GPIO_Pin GPIO_Pin_9 | GPIO_Pin_10;p.GPIO_PuPd G…

卷积神经网络--CNN

1.人工神经网络 神经网络由大量的节点&#xff08;或称“神经元”、“单元”&#xff09;和相互连接而成。每个神经元接受输入的线性组合&#xff0c;进行非线性变换&#xff08;亦称激活函数activation function&#xff09;后输出。每两个节点之间的连接代表加权值&#xff0…

基于WinCE的I2C驱动程序设计

http://www.mcu123.com/news/Article/rtos/WinCE/200607/88.html 引言 随着以计算机技术、通信技术和软件技术为核心的信息技术的迅速发展&#xff0c;嵌入式系统在各行业得到了广泛的应用&#xff0c;极大地推动了行业的渗透性应用。嵌入式系统是“以应用为中心、以计算机技…

poj2965-poj2965-The Pilots Brothers' refrigerator

方法同poj1753&#xff0c;但用在这题就TLE了&#xff0c;以下是TLE版本&#xff1a; Code1#include <stdio.h> 2#include <stdlib.h> 3#include <string.h> 4#define MAXSTATE 65536 5#define MAXSIZE 16 6#define ALLOPEN 0 7 8//队列结构体 9type…

sysctl -p详解

个人一般sysctl -p 或sysctl -a比较多使用 sysctl配置与显示在/proc/sys目录中的内核参数&#xff0e;可以用sysctl来设置或重新设置联网功能&#xff0c;如IP转发、IP碎片去除以及源路由检查等。用户只需要编辑/etc/sysctl.conf文件&#xff0c;即可手工或自动执行由sysctl控…

定制简单的Linux系统

定制简单的Linux系统 制作思路&#xff1a; 新加一块硬盘&#xff0c;设置两个分区&#xff0c;一个存/boot&#xff0c;一个存/&#xff0c;创建文件系统并格式化。要注意&#xff0c;现在我们家的硬盘是要可以拔下来安装到其他机器上使用的&#xff0c;否则就没有意义了。试…

UCOS同步与互斥

代码为老师教授。 /* ********************************************************************************************************* * EXAMPLE CODE * * (c) Copyright 2013; Micrium, Inc.; We…

Spring学习八

1&#xff1a; Tomcat容器四个等级&#xff1f; Container&#xff0c; Engine&#xff0c; Servlet容器&#xff0c; Context 真正管理Servlet的容器是Context容器&#xff1a;一个context对应一个web工程。 <Context path"/projectOne " docBase"D:\proje…

作业六:图像编码相关概念

7.1&#xff0e;信息量&#xff1a;信息源发出的所有消息中该消息出现概率的倒数的对数。信息熵&#xff1a;对信息源X的各符号的自信息量取统计平均。 7.6如图所示&#xff1a;哈夫曼编码最终结果为&#xff1a;X11,X201,X3000,X4001。编码效率为95%。 7.8根据公式&#xff…

linux命令find命令详解

find 查找文件 find 哪里 什么类型 什么名字 -maxdepth 最大的深度 查找目录的最大深度 find -maxdepth 1 -type d -type 找什么类型的 f file 文件 d directory 目录 -name 什么名字 -mtime 根据修改时间找出对应的文件 7 7天前 -7 7天后 find命令一般与 |xargs 一起…

一次次小进步,从毕业开始,你到现在飞跃了几次了,程序人生也不容易?

01. 会写最简单的程序&#xff0c;能编译通过了&#xff0c;是一次飞跃。02. 会写C/S程序了&#xff0c;能用那些常用的控件&#xff0c;对属性事件有了解了&#xff0c;会用了&#xff0c;是一次飞跃。03. 会写B/S程序了&#xff0c;也是一次飞跃。04. 你彻底理解了分层的理念…