bzoj4066: 简单题

链接

bzoj

思路

强制在线。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的更多相关文章

  1. [BZOJ2683][BZOJ4066]简单题

    [BZOJ2683][BZOJ4066]简单题 试题描述 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x ...

  2. Bzoj4066 简单题

    Time Limit: 50 Sec  Memory Limit: 20 MBSubmit: 2185  Solved: 581 Description 你有一个N*N的棋盘,每个格子内有一个整数,初 ...

  3. BZOJ4066 简单题(KD-Tree)

    板子题. #include<iostream> #include<cstdio> #include<cmath> #include<cstdlib> # ...

  4. BZOJ4066:简单题(K-D Tree)

    Description 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作:   命令 参数限制 内容 1 x y A 1<=x,y<=N,A是正整数 ...

  5. 【kd-tree】bzoj4066 简单题

    同p1176. #include<cstdio> #include<cmath> #include<algorithm> using namespace std; ...

  6. [bzoj4066/2683]简单题_KD-Tree

    简单题 bzoj-4066 题目大意:n*n的棋盘,开始为均为0,支持:单点加权值,查询矩阵权值和,强制在线. 注释:$1\le n\le 5\cdot 10^5$,$1\le m \le 2\cdo ...

  7. 【BZOJ4066】简单题(KD-Tree)

    [BZOJ4066]简单题(KD-Tree) 题面 BZOJ 题解 如果这题不卡空间,并且不强制在线的话 显然可以用\(CDQ\)分治做 但是它又卡空间又强制在线,于是我们欢快的来用\(KD-Tree ...

  8. 【BZOJ4066】简单题 KDtree

    [BZOJ4066]简单题 Description 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x,y& ...

  9. Bzoj2683 简单题

    Time Limit: 50 Sec  Memory Limit: 128 MBSubmit: 1071  Solved: 428 Description 你有一个N*N的棋盘,每个格子内有一个整数, ...

随机推荐

  1. JQuery的使用案例(二级联动,隔行换色,轮播图,广告插入)

    JQuery的使用案例 (一)利用JQuery完成省市二级联动 第一步:确定事件(change事件),在绑定的函数里面获取用户选择的省份 第二步:创建二维数组来存储省份和城市 第三步:遍历二维数组中的 ...

  2. 原生js数值开根算法

    不借助Math函数求开根值 1.二分迭代法求n开根后的值 思路: left=0 right=n mid=(left+right)/2 比较mid^2与n大小 =输出: >改变范围,right=m ...

  3. Unity AsyncGPUReadback 接口测试

    Unity2018新加入了该接口,可以做到异步RenderTexture->像素数据和异步的ComputeBuffer.GetData 那么写了几个例子来测试下. 1.RenderTexture ...

  4. xcode模拟器使用常用的命令。

    1.查看模拟器的udid用的 xcrun instruments -s xcrun simctl list 2.启动这个模拟器: xcrun instruments -w 'B39EC2FF-8A8B ...

  5. SAP HANA学习资料大全 Simple Finane + Simple Logisitic [非常完善的学习资料汇总]

    Check out this SDN blog if you plan to write HANA Certification exam http://scn.sap.com/community/ha ...

  6. MybatisGenerator生成SSM的dao层

    官网下载 mybatis generator 下载generator的release版本mybatis-generator-core-1.4.0-bundle.zip https://github.c ...

  7. js中 !==和 !=的区别是什么

    1.比较结果上的区别 !=返回同类型值比较结果. !== 不同类型不比较,且无结果,同类型才比较. 2.比较过程上的区别 != 比较时,若类型不同,会偿试转换类型. !== 只有相同类型才会比较. 3 ...

  8. Android图像滤镜框架GPUImage从配置到应用

    GPUImage简介 GPUImage 是iOS下一个开源的基于GPU的图像处理库,提供各种各样的图像处理滤镜,并且支持照相机和摄像机的实时滤镜.GPUImage for Android是它在Andr ...

  9. idea在maven中引入了jar包依赖,但是编译过程中报出XXX程序包不存在,已解决

    idea在maven中引入了jar包依赖,但是编译过程中报出XXX程序包不存在 1. 报错具体情况 2. Project Structure中的Libraries没有任何红色波浪线 3. 发现自己要引 ...

  10. 利用 AWS DMS 在线迁移 MongoDB 到 Amazon Aurora

    将数据从一种数据库迁移到另一种数据库通常都非常具有挑战性,特别是考虑到数据一致性.应用停机时间.以及源和目标数据库在设计上的差异性等因素.这个过程中,运维人员通常都希望借助于专门的数据迁移(复制)工具 ...