【BZOJ】【4066】简单题(强制在线)
KD-Tree
KD-Tree的进阶姿势戳这里 http://zyfzyf.is-programmer.com/posts/92431.html
为啥有种线段树&平衡树的即视感……(树形结构的相似性?)
每次插入之后,判断下如果某个子树的size>父亲size*0.7,那么重构一下……(替罪羊树的即视感)
查询的时候,如果当前点表示的坐标范围被查询范围完全包含,则直接返回sum;
否则:当前点若在范围内则更新答案,左子树若不全在范围外则递归进入查询,右子树同理(线段树的即视感)
TLE:rebuild的时候,我没清当前点的min和max就直接build子树&Push_up(o)了……活该跪……
UPD:2015-05-23 17:20:26
今天看论文的时候,发现这一做法在陈老师的论文中已经讲出来了……orz,无限ym陈老师!
/**************************************************************
Problem: 4066
User: Tunix
Language: C++
Result: Accepted
Time:28400 ms
Memory:11444 kb
****************************************************************/ //BZOJ 4066
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
typedef long long LL;
inline int getint(){
int r=,v=; char ch=getchar();
for(;!isdigit(ch);ch=getchar()) if (ch=='-') r=-;
for(; isdigit(ch);ch=getchar()) v=v*-''+ch;
return r*v;
}
const int N=,INF=1e9;
/*******************template********************/
int n,m,ans,D,cnt,p[N],tot,root;
struct node{
int d[],mn[],mx[],l,r,D,size,sum,v;
int& operator [] (int x){return d[x];}
}t[N],now;
inline bool cmp(int x,int y){return t[x][D]<t[y][D];} #define L t[o].l
#define R t[o].r
#define mid (l+r>>1)
inline void Push_up(int o){
F(i,,){
t[o].mn[i]=min(t[o].mn[i],min(t[L].mn[i],t[R].mn[i]));
t[o].mx[i]=max(t[o].mx[i],max(t[L].mx[i],t[R].mx[i]));
}
t[o].sum=t[L].sum+t[R].sum+t[o].v;
t[o].size=t[L].size+t[R].size+;
}
inline int build(int l,int r,int dir){
D=dir;
nth_element(p+l,p+mid,p+r+,cmp);
int o=p[mid];
t[o].D=dir;
F(i,,) t[o].mn[i]=t[o].mx[i]=t[o][i];
t[o].sum=t[o].v;
L=l<mid ? build(l,mid-,dir^) : ;
R=mid<r ? build(mid+,r,dir^) : ;
Push_up(o);
return o;
}
inline void dfs(int o){
if (!o) return;
dfs(L);
p[++cnt]=o;
dfs(R);
}
inline void rebuild(int &o){
cnt=;
dfs(o);
o=build(,cnt,t[o].D);
}
inline void Insert(int &o,int dir){
if (!o){
o=++tot; t[o]=now;
F(i,,) t[o].mn[i]=t[o].mx[i]=t[o][i];
t[o].D=dir;
t[o].size=;
t[o].sum=t[o].v;
return;
}
if (now[dir]<t[o][dir]){
Insert(L,dir^);
Push_up(o);
if ((double)t[L].size>(double)t[o].size*0.7) rebuild(o);
}
else{
Insert(R,dir^);
Push_up(o);
if ((double)t[R].size>(double)t[o].size*0.7) rebuild(o);
}
}
int query(int o,int x1,int y1,int x2,int y2){
if (!o) return ;
if (t[o].mn[]>=x1 && t[o].mn[]>=y1 && t[o].mx[]<=x2 && t[o].mx[]<=y2)
return t[o].sum;
else{
int ans=;
if (t[o][]>=x1 && t[o][]<=x2 && t[o][]>=y1 && t[o][]<=y2) ans+=t[o].v;
if (t[L].mn[]>x2 || t[L].mx[]<x1 || t[L].mn[]>y2 || t[L].mx[]<y1) ;
else ans+=query(L,x1,y1,x2,y2);
if (t[R].mn[]>x2 || t[R].mx[]<x1 || t[R].mn[]>y2 || t[R].mx[]<y1) ;
else ans+=query(R,x1,y1,x2,y2);
return ans;
}
}
void print(int o){
if (!o) return;
printf("%d t[o].mn[0]=%d t[o].mn[1]=%d t[o].mx[0]=%d t[o].mx[1]=%d\n",o,t[o].mn[],t[o].mn[],t[o].mx[],t[o].mx[]);
print(L);
print(R);
}
int main(){
#ifndef ONLINE_JUDGE
freopen("4066.in","r",stdin);
freopen("4066.out","w",stdout);
#endif
F(i,,) t[].mn[i]=INF,t[].mx[i]=-INF;
t[].size=t[].sum=t[].v=; n=getint(); int cmd;
while(scanf("%d",&cmd)!=EOF && cmd<){
if (cmd==){
now[]=getint()^ans,now[]=getint()^ans,now.v=getint()^ans;
Insert(root,);
}else{
int x1=getint()^ans,y1=getint()^ans,x2=getint()^ans,y2=getint()^ans;
printf("%d\n",ans=query(root,x1,y1,x2,y2));
}
}
return ;
}
4066: 简单题
Time Limit: 50 Sec Memory Limit: 20 MB
Submit: 297 Solved: 99
[Submit][Status][Discuss]
Description
命令 |
参数限制 |
内容 |
1 x y A |
1<=x,y<=N,A是正整数 |
将格子x,y里的数字加上A |
2 x1 y1 x2 y2 |
1<=x1<= x2<=N 1<=y1<= y2<=N |
输出x1 y1 x2 y2这个矩形内的数字和 |
3 |
无 |
终止程序 |
Input
Output
Sample Input
1 2 3 3
2 1 1 3 3
1 1 1 1
2 1 1 0 7
3
Sample Output
5
HINT
Source
【BZOJ】【4066】简单题(强制在线)的更多相关文章
- bzoj 4066: 简单题 kd-tree
4066: 简单题 Time Limit: 50 Sec Memory Limit: 20 MBSubmit: 234 Solved: 82[Submit][Status][Discuss] De ...
- BZOJ 4066 简单题(KD树)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=4066 [题目大意] 要求维护矩阵内格子加点和矩阵查询 [题解] 往KD树上加权值点,支 ...
- bzoj 4066 简单题——KDtree(带重构)
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4066 带部分重构的KDtree.就是那个替罪羊树思想的. 写了对拍,调了半天,发现忘了 re ...
- bzoj 4066: 简单题 K-D树
题目大意: http://www.lydsy.com/JudgeOnline/problem.php?id=4066 题解 我们把每次的修改操作都当作二维平面上多了一个权值点 对于每组询问可以看做求一 ...
- BZOJ 4066 简单题 ——KD-Tree套替罪羊树
[题目分析] 直接x,y二维轮番划分,暴力即可. 套上替罪羊,打碎重构,对于时间复杂度有了保证. 写起来好麻烦,重构的技巧很棒! [代码] #include <cstdio> #inclu ...
- bzoj 4066: 简单题
#include<cstdio> #include<iostream> #include<cstdlib> #include<algorithm> #d ...
- BZOJ 2683: 简单题
2683: 简单题 Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 913 Solved: 379[Submit][Status][Discuss] ...
- BZOJ 3687: 简单题 bitset
3687: 简单题 Time Limit: 10 Sec Memory Limit: 512 MB[Submit][Status][Discuss] Description 小呆开始研究集合论了,他 ...
- bzoj 4066 & bzoj 2683 简单题 —— K-D树(含重构)
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4066 https://www.lydsy.com/JudgeOnline/problem.p ...
- BZOJ 2683: 简单题(CDQ分治 + 树状数组)
BZOJ2683: 简单题(CDQ分治 + 树状数组) 题意: 你有一个\(N*N\)的棋盘,每个格子内有一个整数,初始时的时候全部为\(0\),现在需要维护两种操作: 命令 参数限制 内容 \(1\ ...
随机推荐
- 纯CSS仿制Google女生节Doodle
看到google今天的女生节Doodle,自己用纯css仿制一个,送给老妈.老婆.女儿. 大家可以点这里在线观看效果,点这里下载收藏效果. 实现原理 1.利用checkbox侦听处理单击事件. 2.单 ...
- CSUOJ 1560 图书管理员的表白方式
Description 小V是中南大学图书馆的图书管理员,每天要整理很多同学们还回来的书.久而久之,他认识了很多常来图书馆的同学,比如说小L.简而言之吧,就是小V喜欢上了小L,并且想在下一次她来还书的 ...
- Docker应用系列(六)| 如何去掉sudo及避免权限问题
一.如何在使用docker时去掉sudo 1.添加账户 $ sudo groupadd docker 2.授权给docker账户 sudo gpasswd -a yourname docker 3.重 ...
- @NamedEntityGraphs --JPA按实体类对象参数中的字段排序问题得解决方法
JPA按实体类对象参数中的字段排序问题得解决方法@Entity @Table(name="complaints") @NamedEntityGraphs({ @NamedEntit ...
- vue.js-过滤器 filters使用详细示例
什么也不说了,直接上干货: 1.首先,获取后台数据到页面,并调用过滤器 在<script>中添加 onRefreshItems (currentPage, perPage) { if (t ...
- 从Table 表中取出第 m 条到第 n 条的记录
* FROM Table id FROM Table )) --从TABLE表中取出第m到n条记录 (Exists版本) * FROM TABLE AS a WHERE Not Exists ( * ...
- Java 的类加载顺序
Java 的类加载顺序 一.加载顺序:先父类后子类,先静态后普通 1.父类的静态成员变量初始化 2.父类的静态代码块 3.子类的静态成员变量初始化 4.子类的静态代码块 5.父类的普通成员变量初始化 ...
- BZOJ 3572: [Hnoi2014]世界树 虚树 树形dp
https://www.lydsy.com/JudgeOnline/problem.php?id=3572 http://hzwer.com/6804.html 写的时候参考了hzwer的代码,不会写 ...
- nginx_lua vs nginx+php 应用场景
在我的印象中很多人还是选择nginx+php这种组合搭配,你的选择是nginx+lua,那么nginx+lua比和php的组合优势在哪里?清无:首先,Nginx+php之间是要有进程之间通信的,这样以 ...
- 20162307 课堂测试 hash
20162307 课堂测试 hash 作业要求 利用除留余数法为下列关键字集合的存储设计hash函数,并画出分别用开放寻址法和拉链法解决冲突得到的空间存储状态(散列因子取0.75) 关键字集合:85, ...