bzoj 3600 没有人的算术——二叉查找树动态标号
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3600
已知 l 和 r 的排名,想快速知道 k 的排名。那么建一个 BIT ,用已知的排名做比较,走到一个位置,就知道自己的排名了。
但这样会让很多点的排名改变。如果用实数表示标号就可以方便地生成一个 “两个排名之间的排名” 而不影响其他点。
BIT 的每个节点代表一个曾经出现过的序列里的点,有自己的实数标号和 ( l , r ) 的组成。
如果把实数标号不是记在 BIT 节点上而是记在序列上,即 v[ i ] 表示序列第 i 个位置的标号,那么修改 k 位置的话,其他曾经调用 k 而组成的点也会收到影响。
所以实数标号记在 BIT 节点上,序列记一个 dy[ ] 表示它对应 BIT 哪个点。序列一个位置的 ( l , r ) 记录成 BIT 上对应的点即可。
用替罪羊树实现平衡。那么在重构的时候要注意,因为 dy[ ] 不会改变,所以 BIT 上点标之间的大小关系也不应改变。只要在回收的时候不是回收 Node 而是回收 int 即可(int 是那个 Node 的角标)。
最开始的时候只有一个 0 。那么 BIT 里唯一的节点表示 “只有一个 0 ” ,其 l 就是自己, r 是 0 (BIT 节点从1标号)。v[ 1 ] 会是 0.5 (如果标号总范围是 ( 0 , 1 ) 的话), v[ 0 ] 会是 0 ,这样就表示了 “没有” 是最小的。
#include<cstdio>
#include<cstring>
#include<algorithm>
#define db double
#define ls Ls[cr]
#define rs Rs[cr]
using namespace std;
int rdn()
{
int ret=;bool fx=;char ch=getchar();
while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
return fx?ret:-ret;
}
const int N=1e5+,M=5e5+;
const db sl=0.75;
int n,rt,tot,c0[M],c1[M],siz[M],dy[N]; db v[M];
int Ls[N<<],Rs[N<<],mx[N<<];
int rcr,rfa,sta[M],top; bool rfx; db rL,rR;
struct Node{
int l,r;
Node(int l=,int r=):l(l),r(r) {}
bool operator< (const Node &b)const
{ return v[l]==v[b.l]?v[r]<v[b.r]:v[l]<v[b.l];}
bool operator== (const Node &b)const
{ return v[l]==v[b.l]&&v[r]==v[b.r];}
}a[M];
int Mx(int u,int v){return u>v?u:v;}
Node Mx(Node u,Node v){return u<v?v:u;}
void build(int l,int r,int cr)
{
mx[cr]=l; if(l==r)return;
int mid=l+r>>;
ls=++tot; build(l,mid,ls);
rs=++tot; build(mid+,r,rs);
}
void mdfy(int l,int r,int cr,int p)
{
if(l==r)return; int mid=l+r>>;
if(p<=mid)mdfy(l,mid,ls,p);
else mdfy(mid+,r,rs,p);
if(a[dy[mx[ls]]]<a[dy[mx[rs]]])mx[cr]=mx[rs];
else mx[cr]=mx[ls];
}
int qry(int l,int r,int cr,int L,int R)
{
if(l>=L&&r<=R)return mx[cr];
int mid=l+r>>;
if(R<=mid)return qry(l,mid,ls,L,R);
if(mid<L)return qry(mid+,r,rs,L,R);
int u=qry(l,mid,ls,L,R), v=qry(mid+,r,rs,L,R);
if(a[dy[u]]<a[dy[v]])return v;
else return u;
}
void chk(int cr,db L,db Md,db R)
{
if(Mx(siz[c0[cr]],siz[c1[cr]])>=siz[cr]*sl)rcr=cr;
if(rcr==c0[cr])
{ rfa=cr; rfx=; rL=L; rR=Md;}
if(rcr==c1[cr])
{ rfa=cr; rfx=; rL=Md; rR=R;}
}
int ins(int &cr,Node k,db L,db R)
{
db Md=(L+R)/;
if(!cr)
{
cr=++tot; a[cr]=k;
v[cr]=Md;siz[cr]=;return cr;
}
if(k==a[cr])return cr;
int ret=; siz[cr]++;
if(k<a[cr])ret=ins(c0[cr],k,L,Md);
else ret=ins(c1[cr],k,Md,R);
chk(cr,L,Md,R); return ret;
}
void del_dfs(int cr)//sorted!!
{
if(c0[cr])del_dfs(c0[cr]);
sta[++top]=cr;
if(c1[cr])del_dfs(c1[cr]);
}
void bld(int l,int r,int &cr,db L,db R)
{
db Md=(L+R)/; int mid=l+r>>;
cr=sta[mid];
c0[cr]=c1[cr]=; siz[cr]=; v[cr]=Md;
if(l<mid)bld(l,mid-,c0[cr],L,Md);
if(mid<r)bld(mid+,r,c1[cr],Md,R);
siz[cr]=siz[c0[cr]]+siz[c1[cr]]+;
}
void rebld(int &cr,db L,db R)
{
top=; del_dfs(cr); bld(,top,cr,L,R);
}
int main()
{
n=rdn();int m=rdn(); char ch[];
tot=;build(,n,);
tot=; ins(rt,Node(,),,);//v[1]=0.5//v[0]=0 is min
for(int i=;i<=n;i++)dy[i]=;
for(int i=,l,r,k;i<=m;i++)
{
scanf("%s",ch); l=rdn();r=rdn();
if(ch[]=='C')
{
k=rdn();
dy[k]=ins(rt,Node(dy[l],dy[r]),,);
Node tp=a[dy[k]];
mdfy(,n,,k);
if(rcr)
{
if(rcr==rt)rebld(rt,,);
else if(!rfx)rebld(c0[rfa],rL,rR);
else rebld(c1[rfa],rL,rR);
rcr=;
}
}
else printf("%d\n",qry(,n,,l,r));
}
return ;
}
bzoj 3600 没有人的算术——二叉查找树动态标号的更多相关文章
- bzoj 3600 没有人的算术 - 替罪羊树 - 线段树
题目都是图片,就不给了,就给链接好了 由于bzoj比较慢,就先给[vjudge传送门] 有兴趣的可以去逛bzoj[bzoj传送门] 题目大意 有n个数a[1],a[2],...,a[n],它们开始都是 ...
- bzoj 3600: 没有人的算术
Description Solution 我们可以给每一个数钦定一个权值 , 这样就可以 \(O(1)\) 比较大小了. 考虑怎么确定权值: 用平衡树来维护 , 我们假设根节点管辖 \([1,2^{6 ...
- 【题解】BZOJ 3600: 没有人的算术——替罪羊树、线段树
题目传送门 题意 具体的自己去上面看吧...反正不是权限题. 简单来说,就是定义了一类新的数,每个数是0或者为 \((x_L, x_R)\) ,同时定义比较大小的方式为:非零数大于零,否则按字典序比较 ...
- 「BZOJ3600」没有人的算术 替罪羊树+线段树
题目描述 过长--不想发图也不想发文字,所以就发链接吧-- 没有人的算术 题解 \(orz\)神题一枚 我们考虑如果插入的数不是数对,而是普通的数,这就是一道傻题了--直接线段树一顿乱上就可以了. 于 ...
- bzoj 2594: [Wc2006]水管局长数据加强版 动态树
2594: [Wc2006]水管局长数据加强版 Time Limit: 25 Sec Memory Limit: 128 MBSubmit: 934 Solved: 291[Submit][Sta ...
- BZOJ 3531 [Sdoi2014]旅行 树链剖分+动态开点线段树
题意 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天面条神教.隐形独角兽教.绝地教都是常见的信仰. 为了方便,我们用 ...
- bzoj3600: 没有人的算术
题意:太难说了..手动去看吧反正不是权限题. 膜拜VFK大爷的神题! 其实一开始思路挺清楚的,如果我们能做到用一个实数去代表"数",这就是裸的动态区间最值查询. 关键是怎么用实数去 ...
- BZOJ 3924: [Zjoi2015]幻想乡战略游戏(动态点分治)
这种动态点分治嘛,GDKOI时听打到了,也有同学讲到了,所以印象比较深刻也就想出来了,然后就在实现方面卡了好久= = 不得不说CLJ说得真的太简单了,实现方面根本没提. 首先我们可以先用树分治构建出这 ...
- BZOJ 2055: 80人环游世界 [上下界费用流]
2055: 80人环游世界 题意:n个点带权图,选出m条路径,每个点经过val[i]次,求最小花费 建图比较简单 s拆点限制流量m 一个点拆成两个,限制流量val[i],需要用上下界 图中有边的连边, ...
随机推荐
- CSS学习笔记之样式规划
大家都知道规范灵活的代码布局对提升程序员开发和后期维护效率至关重要,因为css同一元素可能被不同偏重度的选择器命中,相同元素不同的选择器表达式的样式冲突导致的显示异常,再加上不规范的代码,经常让前端代 ...
- FCC JS基础算法题(6):Truncate a string(截断字符串)
先看一下题目描述: 如果字符串的长度比指定的参数num长,则把多余的部分用...来表示.切记,插入到字符串尾部的三个点号也会计入字符串的长度.但是,如果指定的参数num小于或等于3,则添加的三个点号不 ...
- python scrapy爬虫存储数据库方法带去重步骤
import pymongo import requests import random import time import pymysql db = pymongo.MongoClient()[' ...
- ES6 Module的语法
下面比较一下默认输出和正常输出. // 第一组 export default function crc32() { // 输出 // ... } import crc32 from 'crc32'; ...
- this -- apply或call调用
apply 和 call 这两个方法可以切换函数执行的上下文环境,也就是改变this绑定的对象.apply和call比较类似,区别在于传入参数时一个要求是数组,一个要求是分开传入.所以我们以apply ...
- SQL注入之Sqli-labs系列第九关和第十关(基于时间盲注的注入)
开始挑战第九关(Blind- Time based- Single Quotes- String)和第十关( Blind- Time based- Double Quotes- String) gog ...
- Codeforce 9C - Hexadecimal's Numbers
One beautiful July morning a terrible thing happened in Mainframe: a mean virus Megabyte somehow got ...
- BT原理分析
BT全名為BitTorrent,是一個p2p軟件,你在下載download的同時,也在為其他用戶提供上傳upload,因為大家是「互相幫助」,所以不會隨著用戶數的增加而降低下載速度. 其實跟ED也十分 ...
- Gtk-WARNING **: cannot open display: :0.0
Gtk-WARNING **: cannot open display: :0.0 https://blog.csdn.net/Rong_Toa/article/details/80365932
- linux下history命令显示执行时间
想在输入history命令之后,显示自己历史的命令执行的时间,需要在用户目录下~/.bashrc的文件末尾追加添加如下几行 之前一直想看一下自己历史命令执行的时间,找了很多教程都没有卵用,最终参考了如 ...