BZOJ4066:简单题
浅谈\(K-D\) \(Tree\):https://www.cnblogs.com/AKMer/p/10387266.html
题目传送门:https://lydsy.com/JudgeOnline/problem.php?id=4066
裸的强制在线插入点,像替罪羊树那样重建不平衡的子树即可。在回收要重建的子树的节点的时候记得把结点信息初始化掉。
时间复杂度:\(O(n\sqrt{n})\)
空间复杂度:\(O(n)\)
代码如下:
#include <cstdio>
#include <algorithm>
using namespace std;
const double alpha=0.75;
const int maxn=2e5+5,inf=2e9;
bool need_rebuild;
int n,lstans,pps,x,y,A,x1,x2,y1,y2;
int read() {
int x=0,f=1;char ch=getchar();
for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';
return x*f;
}
struct kd_tree {
int id[maxn];
int tot,root,cnt;
struct point {
int c[2],mn[2],mx[2];
int val,sum,siz,ls,rs;
bool operator<(const point &a)const {
return c[pps]<a.c[pps];
}
}p[maxn],node[maxn];
void prepare() {
p[0].mn[0]=p[0].mn[1]=inf;
p[0].mx[0]=p[0].mx[1]=-inf;
}
void update(int u) {
int ls=p[u].ls,rs=p[u].rs;
p[u].siz=p[ls].siz+1+p[rs].siz;
p[u].sum=p[ls].sum+p[u].val+p[rs].sum;
for(int i=0;i<2;i++) {
int mn=min(p[ls].mn[i],p[rs].mn[i]);
p[u].mn[i]=min(p[u].c[i],mn);
int mx=max(p[ls].mx[i],p[rs].mx[i]);
p[u].mx[i]=max(p[u].c[i],mx);
}
}
void ins(int &u,int d) {
if(!u) {
u=++tot;need_rebuild=1;
p[u].mn[0]=p[u].mx[0]=p[u].c[0]=x;
p[u].mn[1]=p[u].mx[1]=p[u].c[1]=y;
p[u].val=p[u].sum=A,p[u].siz=1;return;
}
if(p[u].c[0]==x&&p[u].c[1]==y) {
p[u].val+=A,p[u].sum+=A;return;
}
int tmp=d?y:x;
if(tmp<p[u].c[d])ins(p[u].ls,d^1);
else ins(p[u].rs,d^1);
update(u);
}
void recycle(int u) {
if(!u)return;
int ls=p[u].ls,rs=p[u].rs;
p[u].siz=1,p[u].ls=p[u].rs=p[u].sum=0;
id[++cnt]=u,node[cnt]=p[u];
recycle(ls),recycle(rs);
}
int rebuild(int l,int r,int d) {
int mid=(l+r)>>1,u=id[mid];pps=d;
nth_element(node+l,node+mid,node+r+1);
p[u]=node[mid];
if(l<mid)p[u].ls=rebuild(l,mid-1,d^1);
if(r>mid)p[u].rs=rebuild(mid+1,r,d^1);
update(u);return u;
}
void check(int &u,int d) {
int siz1=p[p[u].ls].siz,siz2=p[p[u].rs].siz;
if(max(siz1,siz2)>1.0*p[u].siz*alpha) {
recycle(u),u=rebuild(1,cnt,d),cnt=0;return;
}
if(p[u].c[0]==x&&p[u].c[1]==y)return;
int tmp=d?y:x;
if(tmp<p[u].c[d])check(p[u].ls,d^1);
else check(p[u].rs,d^1);
}
void query(int u) {
if(!u)return;
if(x2<p[u].mn[0]||x1>p[u].mx[0])return;
if(y2<p[u].mn[1]||y1>p[u].mx[1])return;
bool bo1=(x1<=p[u].mn[0]&&p[u].mx[0]<=x2);
bool bo2=(y1<=p[u].mn[1]&&p[u].mx[1]<=y2);
if(bo1&&bo2) {lstans+=p[u].sum;return;}
bo1=(x1<=p[u].c[0]&&p[u].c[0]<=x2);
bo2=(y1<=p[u].c[1]&&p[u].c[1]<=y2);
if(bo1&&bo2)lstans+=p[u].val;
if(p[u].ls)query(p[u].ls);
if(p[u].rs)query(p[u].rs);
}
}T;
int main() {
n=read();T.prepare();
while(1) {
int opt=read();
if(opt==1) {
x=read()^lstans,y=read()^lstans,A=read()^lstans;
T.ins(T.root,0);if(need_rebuild)T.check(T.root,0);
}
if(opt==2) {
x1=read()^lstans,y1=read()^lstans;
x2=read()^lstans,y2=read()^lstans;
lstans=0,T.query(T.root),printf("%d\n",lstans);
}
if(opt==3)break;
}
return 0;
}
BZOJ4066:简单题的更多相关文章
- [BZOJ2683][BZOJ4066]简单题
[BZOJ2683][BZOJ4066]简单题 试题描述 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x ...
- bzoj4066: 简单题 K-Dtree
bzoj4066: 简单题 链接 bzoj 思路 强制在线.k-dtree. 卡常啊.空间开1e6就T了. 代码 #include <bits/stdc++.h> #define my_m ...
- 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> # ...
- 【kd-tree】bzoj4066 简单题
同p1176. #include<cstdio> #include<cmath> #include<algorithm> using namespace std; ...
- BZOJ4066:简单题(K-D Tree)
Description 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x,y<=N,A是正整数 ...
- 【BZOJ4066】简单题(KD-Tree)
[BZOJ4066]简单题(KD-Tree) 题面 BZOJ 题解 如果这题不卡空间,并且不强制在线的话 显然可以用\(CDQ\)分治做 但是它又卡空间又强制在线,于是我们欢快的来用\(KD-Tree ...
- 【BZOJ4066】简单题 KDtree
[BZOJ4066]简单题 Description 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x,y& ...
- [bzoj4066/2683]简单题_KD-Tree
简单题 bzoj-4066 题目大意:n*n的棋盘,开始为均为0,支持:单点加权值,查询矩阵权值和,强制在线. 注释:$1\le n\le 5\cdot 10^5$,$1\le m \le 2\cdo ...
- BZOJ 2683: 简单题
2683: 简单题 Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 913 Solved: 379[Submit][Status][Discuss] ...
随机推荐
- @MarkFan 口语练习录音 20140423 [风雨哈佛路.Homeless To Harvard口语录音]
世界在转动,你只是一粒尘埃 没有你,世界照样在转 现实不会按照你的意识去改变的 一些人的需求 一些人的意志要比你更强 严酷的生活会让人不知所措 所以他们久久地困在挫败中 我们生气地抱怨,而对整体的形势 ...
- node做验证码
使用了ccap插件 1.安装: 通用方法:npm install ccap 2. cnst ccap= require('ccap')({ width: 128, height: 40, offset ...
- windows7 安装Apache2时出现failed to open the winNT service manager 提示
因为电脑实在太慢了,C盘的空间所剩无几,要想再安装大一点的软件的话,可能性很小.加之系统已经好久没有重装过了,于是重新安装windows7旗舰版,系统装好后,免不了一堆软件的重装和开发环境配置,首要的 ...
- poj 1265 Area【计算几何:叉积计算多边形面积+pick定理计算多边形内点数+计算多边形边上点数】
题目:http://poj.org/problem?id=1265 Sample Input 2 4 1 0 0 1 -1 0 0 -1 7 5 0 1 3 -2 2 -1 0 0 -3 -3 1 0 ...
- 轻量级的Java 开发框架 Spring
Spring是一个开源框架,Spring是于2003 年兴起的一个轻量级的Java 开发框架,由Rod Johnson 在其著作Expert One-On-One J2EE Development a ...
- 更新gcc_Install gcc 4.7.x/4.8.x on CentOS
参考一:How to Install gcc 4.7.x/4.8.x on CentOS 参考二:CentOS 升级 gcc 和 g++ 的方法 参考三:linux yum命令详解 参考四:每天一个l ...
- Kafka详解一:Kafka简介
问题导读 1.Kafka有何特性?2.Kafka有哪些组件? 背景: 当今社会各种应用系统诸如商业.社交.搜索.浏览等像信息工厂一样不断的生产出各种信息,在大数据时代,我们面临如下几个挑战: ...
- 修改SpringBoot 默认的小叶子图标
Springboot 项目,在浏览器中访问时,浏览器上导航栏的图标是一片绿色的叶子,我们可以修改它. 将格式为.ico的图片放入以下任一项目文件夹即可.但,图片命名必须为favicon.ico 1.类 ...
- Spring初学之bean之间的关系和bean的作用域
一.bean之间的关系 <?xml version="1.0" encoding="UTF-8"?> <beans xmlns="h ...
- SPOJ375 QTREE - Query on a tree
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ ...