bzoj4066: 简单题 K-Dtree
bzoj4066: 简单题
链接
思路
强制在线。k-dtree。
卡常啊。空间开1e6就T了。
代码
#include <bits/stdc++.h>
#define my_min(a,b) (a<b?a:a=b)
#define my_max(a,b) (a>b?a:a=b)
using namespace std;
const double alpha=0.75;
const int N=2e5+7;
int read() {
int x=0,f=1;
char s=getchar();
for(; s>'9'||s<'0'; s=getchar()) if(s=='-') f=-1;
for(; s>='0'&&s<='9'; s=getchar()) x=x*10+s-'0';
return x*f;
}
int WD,lastans;
struct point {
int x[2],v;
bool operator < (const point &b) const {
return x[WD]<b.x[WD];
}
} p[N];
struct node {
int ls,rs,mi[2],ma[2],sum,siz;
point a;
} e[N];
int rub[N],top,cnt;
int newnode() {
if(top) return rub[top--];
return ++cnt;
}
void update(int x,int y) {
my_min(e[x].mi[0],e[y].mi[0]);
my_min(e[x].mi[1],e[y].mi[1]);
my_max(e[x].ma[0],e[y].ma[0]);
my_max(e[x].ma[1],e[y].ma[1]);
}
void up(int u) {
int ls=e[u].ls,rs=e[u].rs;
e[u].mi[0]=e[u].ma[0]=e[u].a.x[0];
e[u].mi[1]=e[u].ma[1]=e[u].a.x[1];
e[u].sum=e[ls].sum+e[rs].sum+e[u].a.v;
e[u].siz=e[ls].siz+e[rs].siz+1;
if(ls) update(u,ls);
if(rs) update(u,rs);
}
int build(int l,int r,int wd) {
if(l>r) return 0;
int mid=(l+r)>>1;
WD=wd;
nth_element(p+l,p+mid,p+r+1);
int u=newnode();
e[u].a=p[mid];
e[u].ls=build(l,mid-1,wd^1);
e[u].rs=build(mid+1,r,wd^1);
up(u);
return u;
}
void pia(int u,int num) {
if(e[u].ls) pia(e[u].ls,num);
p[num+e[e[u].ls].siz+1]=e[u].a;
rub[++top]=u;
if(e[u].rs) pia(e[u].rs,num+e[e[u].ls].siz+1);
}
int check(int &u,int wd) {
if(alpha*e[u].siz<e[e[u].ls].siz||alpha*e[u].siz<e[e[u].rs].siz)
pia(u,0),u=build(1,e[u].siz,wd);
}
void insert(point a,int &rt,int wd) {
if(!rt) {
rt=newnode();
e[rt].a=a;
e[rt].ls=e[rt].rs=0;
up(rt);
return;
}
if(e[rt].a.x[wd]<a.x[wd]) insert(a,e[rt].rs,wd^1);
else insert(a,e[rt].ls,wd^1);
up(rt);
check(rt,wd);
}
int pd(node a,int x,int y,int X,int Y) {
if(x<=a.mi[0]&&a.ma[0]<=X&&y<=a.mi[1]&&a.ma[1]<=Y) return 1;
if(a.ma[0]<x||X<a.mi[0]) return 0;
if(a.ma[1]<y||Y<a.mi[1]) return 0;
return 2;
}
void query(int rt,int x,int y,int X,int Y) {
if(x<=e[rt].a.x[0]&&e[rt].a.x[0]<=X
&&y<=e[rt].a.x[1]&&e[rt].a.x[1]<=Y) lastans+=e[rt].a.v;
int ls=e[rt].ls,rs=e[rt].rs;
if(ls) {
int PD=pd(e[ls],x,y,X,Y);
if(PD==1) lastans+=e[ls].sum;
if(PD==2) query(ls,x,y,X,Y);
}
if(rs) {
int PD=pd(e[rs],x,y,X,Y);
if(PD==1) lastans+=e[rs].sum;
if(PD==2) query(rs,x,y,X,Y);
}
}
int main() {
int rt=0,n=read();
while(233) {
int opt=read();
if(opt==1) {
point a;
a.x[0]=read()^lastans,a.x[1]=read()^lastans,a.v=read()^lastans;
insert(a,rt,0);
} else if(opt==2) {
int x=read()^lastans,y=read()^lastans,X=read()^lastans,Y=read()^lastans;
lastans=0;
query(rt,x,y,X,Y);
printf("%d\n",lastans);
} else break;
}
return 0;
}
bzoj4066: 简单题 K-Dtree的更多相关文章
- [BZOJ2683][BZOJ4066]简单题
[BZOJ2683][BZOJ4066]简单题 试题描述 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x ...
- Bzoj4066 简单题
Time Limit: 50 Sec Memory Limit: 20 MBSubmit: 2185 Solved: 581 Description 你有一个N*N的棋盘,每个格子内有一个整数,初 ...
- BZOJ4066 简单题(KD-Tree)
板子题. #include<iostream> #include<cstdio> #include<cmath> #include<cstdlib> # ...
- BZOJ4066:简单题(K-D Tree)
Description 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x,y<=N,A是正整数 ...
- 【kd-tree】bzoj4066 简单题
同p1176. #include<cstdio> #include<cmath> #include<algorithm> using namespace std; ...
- [bzoj4066/2683]简单题_KD-Tree
简单题 bzoj-4066 题目大意:n*n的棋盘,开始为均为0,支持:单点加权值,查询矩阵权值和,强制在线. 注释:$1\le n\le 5\cdot 10^5$,$1\le m \le 2\cdo ...
- 【BZOJ4066】简单题(KD-Tree)
[BZOJ4066]简单题(KD-Tree) 题面 BZOJ 题解 如果这题不卡空间,并且不强制在线的话 显然可以用\(CDQ\)分治做 但是它又卡空间又强制在线,于是我们欢快的来用\(KD-Tree ...
- 【BZOJ4066】简单题 KDtree
[BZOJ4066]简单题 Description 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x,y& ...
- Bzoj2683 简单题
Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 1071 Solved: 428 Description 你有一个N*N的棋盘,每个格子内有一个整数, ...
随机推荐
- 【spring boot】【redis】spring boot基于redis的LUA脚本 实现分布式锁
spring boot基于redis的LUA脚本 实现分布式锁[都是基于redis单点下] 一.spring boot 1.5.X 基于redis 的 lua脚本实现分布式锁 1.pom.xml &l ...
- K8S学习笔记之K8S日志搜集实战
详细参考这篇文章,几乎覆盖了了K8S的各种日志搜集方案 https://juejin.im/post/5b6eaef96fb9a04fa25a0d37#heading-8
- K8S学习笔记之Grafana App for Kubernetes的配置
Grafana有一套针对Kubernetes监控的APP,和Grafana-Zabbix App类似,但是配置咋一看比较麻烦,主要参数都是来自K8S. 这款APP的详细介绍请参考Grafana App ...
- listener中@Autowired无法注入bean的一种解决方法
背景:使用监听器处理业务,需要使用自己的service方法: 错误:使用@Autowired注入service对象,最终得到的为null: 原因:listener.fitter都不是Spring容器管 ...
- js浏览器对象模型【BOM】(十三)
一.时间定时器1.超时调用setTimeout(fun,time) [返回一个唯一标识该超时调用的ID数值]参数:fun:要执行的函数time:设置第多少毫秒后执行fun函数 clearTime ...
- vs2017 添加 mysql EF实体数据模型闪退
1.查看vs2017安装路径找到Mysql.Data.dll版本号与MySQL Connector Net版本是否一致 历史版本下载地址 http://mysql.inspire.net.nz/Dow ...
- JavaScript 流程控制(一)顺序结构与分支结构
语句:语句可以理解为一个行为,循环语句和判断语句就是典型的语句.一个程序有很多个语句组成,一般情况下分号;分割一个一个的语句:如果省略分号,则由解析器确定语句的结尾(不推荐使用) 一.流程控制 流程控 ...
- DataPipeline CTO陈肃:构建批流一体数据融合平台的一致性语义保证
文 | 陈肃 DataPipelineCTO 交流微信 | datapipeline2018 本文完整PPT获取 | 关注公众号后,后台回复“陈肃” 首先,本文将从数据融合角度,谈一下DataPipe ...
- MYSQL使用mysqldump导出表的部分数据
MySQLdump是MySQL自带的导出数据工具,通常我们用它来导出MySQL中,但是有时候我们需要导出MySQL数据库中某个表的部分数据,这时该怎么办呢? mysqldump命令中带有一个 --wh ...
- spring boot 的request.getServletContext().getRealPath路径获取问题
默认情况下springboot中request.getServletContext().getRealPath 返回的是一个临时文件夹的地址 通过查看源代码 位置在 org.springframewo ...