LG4848 崂山白花蛇草水
崂山白花蛇草水
神犇 Aleph 在 SDOI Round2 前立了一个 flag:如果进了省队,就现场直播喝崂山白花蛇草水。凭借着神犇 Aleph 的实力,他轻松地进了山东省省队,现在便是他履行诺言的时候了。蒟蒻 Bob 特地为他准备了 999,999,999,999,999,999 瓶崂山白花蛇草水,想要灌神犇 Aleph。神犇 Aleph 求(跪着的)蒟蒻 Bob 不要灌他,由于神犇 Aleph 是神犇,蒟蒻 Bob 最终答应了他的请求,但蒟蒻 Bob 决定将计就计,也让神犇 Aleph 回答一些问题。
具体说来,蒟蒻 Bob 会在一个宽敞的广场上放置一些崂山白花蛇草水(可视为二维平面上的一些整点),然后询问神犇 Aleph 在矩形区域\(x_1\le x\le x_2,y_1\le y\le y_2\)中,崂山白花蛇草水瓶数第 \(k\) 多的是多少。为了避免麻烦,蒟蒻 Bob 不会在同一个位置放置两次或两次以上的崂山白花蛇草水,但蒟蒻 Bob 想为难一下神犇Aleph,希望他能在每次询问时立刻回 答出答案。
神犇 Aleph 不屑于做这种问题,所以把这个问题交给了你。
对于所有数据,\(n\le500000,q\le100000,1\le x, y\le n,1\le v\le 10^9,1\le x_1\le x_2\le n1,1\le y_1\le y_2\le n,1\le k\le q\)。
Claris的题解
外层维护权值线段树,内层维护kd-tree。
修改的时候只往右儿子里插入(因为只会查询右儿子的矩形内点的个数),不平衡的时候替罪羊式重构。
查询的时候在外层线段树上走,在内层kd-tree上查询矩形内点数即可。
\(Q(n)=\sqrt{n}+Q(n/2)=O(\sqrt{n})\)
\(C(n)=\log^2 n + C(n/2)=O(\log^3 n)\)
时间复杂度\(O(q\log v (\log^2 q+\sqrt{q}))\)。
代码
终于找到我梦寐以求的好的替罪羊树的写法了,记录插入进去的那条链,先判深度,可能不平衡再重构一次。
插入和重构分开进行便可去掉冗余操作。
#include<bits/stdc++.h>
#define rg register
#define il inline
#define co const
template<class T>il T read(){
rg T data=0,w=1;
rg char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') w=-1;
ch=getchar();
}
while(isdigit(ch))
data=data*10+ch-'0',ch=getchar();
return data*w;
}
template<class T>il T read(rg T&x){
return x=read<T>();
}
typedef long long ll;
using namespace std;
co double A=0.8;
co int N=1e5+1,M=16e5;
int n,ans,cmp_d,op,X1,Y1,X2,Y2,k;
int tmp[N],deep,need[N],cnt,cur;
int T[M],l[M],r[M],tot=1;
struct node{int d[2],l,r,Max[2],Min[2],size;}t[M];
bool cmp(int a,int b){return t[a].d[cmp_d]<t[b].d[cmp_d];}
void umax(int&a,int b){if(a<b)a=b;}
void umin(int&a,int b){if(a>b)a=b;}
void up(int x){
t[x].size=1+t[t[x].l].size+t[t[x].r].size;
if(t[x].l){
umax(t[x].Max[0],t[t[x].l].Max[0]),umax(t[x].Max[1],t[t[x].l].Max[1]);
umin(t[x].Min[0],t[t[x].l].Min[0]),umin(t[x].Min[1],t[t[x].l].Min[1]);
}
if(t[x].r){
umax(t[x].Max[0],t[t[x].r].Max[0]),umax(t[x].Max[1],t[t[x].r].Max[1]);
umin(t[x].Min[0],t[t[x].r].Min[0]),umin(t[x].Min[1],t[t[x].r].Min[1]);
}
}
int build(int l,int r,int D){
int mid=(l+r)>>1;
cmp_d=D;
nth_element(need+l,need+mid,need+r+1,cmp);
int x=need[mid];
t[x].Max[0]=t[x].Min[0]=t[x].d[0];
t[x].Max[1]=t[x].Min[1]=t[x].d[1];
t[x].l=l<mid?build(l,mid-1,D^1):0;
t[x].r=mid<r?build(mid+1,r,D^1):0;
up(x);
return x;
}
void dfs(int x){if(x) need[++cnt]=x,dfs(t[x].l),dfs(t[x].r);}
void ins(int&root,int now){
if(!root) {root=now;return;}
for(int D=deep=0,x=root;;D^=1){
tmp[++deep]=x;
umax(t[x].Max[0],t[now].Max[0]),umax(t[x].Max[1],t[now].Max[1]);
umin(t[x].Min[0],t[now].Min[0]),umin(t[x].Min[1],t[now].Min[1]);
++t[x].size;
if(t[now].d[D]>=t[x].d[D]){
if(!t[x].r) {t[x].r=now;break;}
else x=t[x].r;
}
else{
if(!t[x].l) {t[x].l=now;break;}
else x=t[x].l;
}
}
tmp[++deep]=now;
if(deep<log(t[root].size)/log(1/A)) return;
while(t[t[now].l].size<A*t[now].size&&t[t[now].r].size<A*t[now].size) now=tmp[--deep];
if(!now) return;
if(now==root){
cnt=0,dfs(root);
root=build(1,cnt,0);
return;
}
int y=tmp[--deep];
cnt=0,dfs(now);
int k=build(1,cnt,deep&1);
if(t[y].l==now) t[y].l=k;
else t[y].r=k;
}
void ask(int x){
if(!x||t[x].Max[0]<X1||t[x].Min[0]>X2||t[x].Max[1]<Y1||t[x].Min[1]>Y2||ans>=k) return;
if(t[x].Min[0]>=X1&&t[x].Max[0]<=X2&&t[x].Min[1]>=Y1&&t[x].Max[1]<=Y2) {ans+=t[x].size;return;}
if(t[x].d[0]>=X1&&t[x].d[0]<=X2&&t[x].d[1]>=Y1&&t[x].d[1]<=Y2) ++ans;
ask(t[x].l),ask(t[x].r);
}
void add(){
int x=1,a=1,b=1e9,mid,flag=1;
while(1){
if(flag){
++cur;
t[cur].Max[0]=t[cur].Min[0]=t[cur].d[0]=X1;
t[cur].Max[1]=t[cur].Min[1]=t[cur].d[1]=Y1;
t[cur].size=1;
ins(T[x],cur);
}
if(a==b) return;
mid=(a+b)>>1;
if(k<=mid){
if(!l[x]) l[x]=++tot;
x=l[x],b=mid,flag=0;
}
else{
if(!r[x]) r[x]=++tot;
x=r[x],a=mid+1,flag=1;
}
}
}
void query(){
ans=0,ask(T[1]);
if(ans<k){
puts("NAIVE!ORZzyz.");
ans=0;
return;
}
int x=1,a=1,b=1e9,mid;
while(a<b){
mid=(a+b)>>1;
ans=0,ask(T[r[x]]);
if(ans>=k) x=r[x],a=mid+1;
else k-=ans,x=l[x],b=mid;
}
printf("%d\n",ans=a);
}
int main(){
// freopen(".in","r",stdin);
// freopen(".out","w",stdout);
read(n),read(n);
while(n--){
read(op),read(X1),read(Y1);X1^=ans,Y1^=ans;
if(op==1){
read(k);k^=ans;
add();
}
else{
read(X2),read(Y2),read(k);X2^=ans,Y2^=ans,k^=ans;
query();
}
}
return 0;
}
LG4848 崂山白花蛇草水的更多相关文章
- bzoj 4605: 崂山白花蛇草水
Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实 力,他轻松地进了山东省省队,现在便是他履行诺言的时 ...
- BZOJ4605:崂山白花蛇草水
浅谈\(K-D\) \(Tree\):https://www.cnblogs.com/AKMer/p/10387266.html 题目传送门:https://lydsy.com/JudgeOnline ...
- 【BZOJ4605】崂山白花蛇草水 权值线段树+kd-tree
[BZOJ4605]崂山白花蛇草水 Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实力,他轻松地进了 ...
- 【bzoj4605】崂山白花蛇草水 权值线段树套KD-tree
题目描述 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实力,他轻松地进了山东省省队,现在便是他履行诺言的时候了.蒟蒻Bob ...
- 崂山白花蛇草水 权值线段树套KDtree
Description 神犇Aleph在SDOI Round2前立了一个flag:如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇Aleph的实 力,他轻松地进了山东省省队,现在便是他履行诺言的时 ...
- bzoj4605: 崂山白花蛇草水 权值线段树套KDtree
bzoj4605: 崂山白花蛇草水 链接 bzoj loj 思路 强制在线,那就权值线段树套KDtree好了,没啥好讲的. KDtree要加平衡因子来重构.另外,那水真难喝. 错误 树套树一边写过了, ...
- 洛谷P4848 崂山白花蛇草水 权值线段树+KDtree
题目描述 神犇 \(Aleph\) 在 \(SDOI\ Round2\) 前立了一个 \(flag\):如果进了省队,就现场直播喝崂山白花蛇草水.凭借着神犇 \(Aleph\) 的实力,他轻松地进了山 ...
- BZOJ4605 : 崂山白花蛇草水
外层维护权值线段树,内层维护kd-tree. 修改的时候只往右儿子里插入,不平衡的时候替罪羊式重构. 查询的时候在外层线段树上走,在内层kd-tree上查询矩形内点数即可. 时间复杂度$O(q\log ...
- BZOJ 4605 崂山白花蛇草水(权值线段树+KD树)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=4605 [题目大意] 操作 1 x y k 表示在点(x,y)上放置k个物品, 操作 2 ...
随机推荐
- maven 引入jar包
问题描述:自己的项目需要引入jar包,已知jar包名字,怎么在maven中添加依赖,使其能自动导入? 第一次使用:本文作为记录! 首先,找到maven仓库的网址!如下: http://mvnrepos ...
- UVA-1623 Enter The Dragon (贪心)
题目大意:有n个装满水的湖,m天.每天可能下雨也可能晴天,只要下雨就会把湖填满,若已满,则发洪水.有一台只能在晴天使用的抽水机,每次抽水只能抽一个湖,并且全部抽光.问是否存在一种使得不发洪水的抽水方案 ...
- 在Linux下安装JDK环境
解压java安装包到指定目录下,这里为/data/software/java/, 然后在/etc/profile添加以下内容: JAVA_HOME =/data/software/java/jdk1. ...
- IOS-Quartz2D
一.画基本图形 // // BWView.m // IOS_0221_Quartz2D画矩形 // // Created by ma c on 16/2/21. // Copyright © 2016 ...
- tkinter模块常用参数
tkinter模块常用参数 1.使用tkinter.Tk() 生成主窗口(root=tkinter.Tk()):root.title('标题名') 修改框体的名字,也可在创建时使用class ...
- RabbitMQ消息队列(九)RPC开始应用吧
一 简单应用 RPC——远程过程调用,通过网络调用运行在另一台计算机上的程序的函数\方法,是构建分布式程序的一种方式.RabbitMQ是一个消息队列系统,可以在程序之间收发消息.利用RabbitMQ可 ...
- linux中的redis缓存服务器
Linux中的Redis缓存服务器 一.Redis基础部分: 1.redis介绍与安装比mysql快10倍以上 *****************redis适用场合**************** 1 ...
- [C#]ref,out关键字的作用
ref是传递参数的地址,out是返回值,两者有一定的相同之处,不过也有不同点. 使用ref前必须对变量赋值,out不用 out的函数会清空变量,即使变量已经赋值也不行,退出函数时所有out引用的变量都 ...
- gitblit-禁用用户
使用管理员账号(此处演示使用admin)登录到系统 右上角,找到用户菜单 再用户列表,找到要进行权限管理的用户,(演示alex) 禁用用户,取消权限 验证已经禁用 pull 权限被禁止
- Java中用反码表示整数,及int到byte的转换
Java中的整数都是通过补码表示的,当整数为正数时,补码和原码相同,当整数为负数时,补码为原码按位取反(符号位除外),然后加1. 通过一些例子可以比较明显的看出: int a = -1; System ...