[BZOJ2683][BZOJ4066]简单题

试题描述

你有一个N*N的棋盘，每个格子内有一个整数，初始时的时候全部为0，现在需要维护两种操作：

命令

参数限制

内容

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

无

终止程序

输入

输入文件第一行一个正整数N。
接下来每行一个操作。每条命令除第一个数字之外，
均要异或上一次输出的答案last_ans，初始时last_ans=0。（BZOJ2683不需要强制在线）

输出

对于每个2操作，输出一个对应的答案。

输入示例（BZOJ2683）





输出示例（BZOJ2683）








输入示例（BZOJ4066）








输出示例（BZOJ4066）








数据规模及约定


1<=N<=500000,操作数不超过200000个，内存限制20M，保证答案在int范围内并且解码之后数据仍合法。


题解


只有权限号可以享受的双倍经验题（雾）。


2683这道数据水，最裸的kd树就卡时限过了。。。


4066这题数据强，需要加一个定期重构，就是加入的元素达到了某一些固定值就暴力把整棵kd树重新构造一遍，这样就不会被卡成“n2 的优秀算法了”。


2683（不强制在线无定期重构慢的要死）：




#include <iostream>

#include <cstdio>

#include <algorithm>

#include <cmath>

#include <stack>

#include <vector>

#include <queue>

#include <cstring>

#include <string>

#include <map>

#include <set>

using namespace std;

const int BufferSize = 1 << 16;

char buffer[BufferSize], *Head, *Tail;

inline char Getchar() {

    if(Head == Tail) {

        int l = fread(buffer, 1, BufferSize, stdin);

        Tail = (Head = buffer) + l;

    }

    return *Head++;

}

int read() {

    int x = 0, f = 1; char c = Getchar();

    while(!isdigit(c)){ if(c == '-') f = -1; c = Getchar(); }

    while(isdigit(c)){ x = x * 10 + c - '0'; c = Getchar(); }

    return x * f;

}

#define maxn 200010

#define oo 2147483647

int root, ToT, lc[maxn], rc[maxn];

struct Node {

	int x[2], mx[2], mn[2], val, sum;

	bool operator == (const Node& t) const { return x[0] == t.x[0] && x[1] == t.x[1]; }

} nodes[maxn];

Node x, y;

void maintain(int o) {

	int l = lc[o], r = rc[o];

	for(int i = 0; i < 2; i++) {

		nodes[o].mx[i] = max(max(nodes[l].mx[i], nodes[r].mx[i]), nodes[o].x[i]);

		nodes[o].mn[i] = min(min(nodes[l].mn[i], nodes[r].mn[i]), nodes[o].x[i]);

	}

	nodes[o].sum = nodes[l].sum + nodes[r].sum + nodes[o].val;

//	printf("maintain(%d): %d %d %d %d %d\n", o, nodes[o].sum, nodes[o].mx[0], nodes[o].mn[0], nodes[o].mx[1], nodes[o].mn[1]);

	return ;

}

void add(int& o, bool cur) {

	if(!o){ nodes[o = ++ToT] = x; return maintain(o); }

	if(nodes[o] == x){ nodes[o].val += x.val; nodes[o].sum += x.val; return maintain(o); }

	add(x.x[cur] < nodes[o].x[cur] ? lc[o] : rc[o], cur ^ 1);

	return maintain(o);

}

bool all(int o) { return x.x[0] <= nodes[o].mn[0] && nodes[o].mx[0] <= y.x[0] && x.x[1] <= nodes[o].mn[1] && nodes[o].mx[1] <= y.x[1]; }

bool has(int o) { return !(nodes[o].mx[0] < x.x[0] || nodes[o].mn[0] > y.x[0] || nodes[o].mx[1] < x.x[1] || nodes[o].mn[1] > y.x[1]); }

int query(int o) {

	if(!o) return 0;

	int ans = 0;

	if(all(lc[o])) ans += nodes[lc[o]].sum;

	else if(has(lc[o])) ans += query(lc[o]);

	if(all(rc[o])) ans += nodes[rc[o]].sum;

	else if(has(rc[o])) ans += query(rc[o]);

	int nx = nodes[o].x[0], ny = nodes[o].x[1];

	if(x.x[0] <= nx && nx <= y.x[0] && x.x[1] <= ny && ny <= y.x[1]) ans += nodes[o].val;

	return ans;

}

int main() {

	nodes[0].mx[0] = nodes[0].mx[1] = -oo;

	nodes[0].mn[0] = nodes[0].mn[1] = oo;

	nodes[0].val = nodes[0].sum = 0;

	int n = read();

	n = read();

	while(n < 3) {

		if(n == 1) {

			x.x[0] = read(); x.x[1] = read(); x.val = read();

			add(root, 1);

		}

		if(n == 2) {

			x.x[0] = read(); x.x[1] = read(); y.x[0] = read(); y.x[1] = read();

			printf("%d\n", query(root));

		}

		n = read();

	}

	return 0;

}





4066（强制在线加定期重构但还是很慢= =）：




#include <iostream>

#include <cstdio>

#include <algorithm>

#include <cmath>

#include <stack>

#include <vector>

#include <queue>

#include <cstring>

#include <string>

#include <map>

#include <set>

using namespace std;

const int BufferSize = 1 << 16;

char buffer[BufferSize], *Head, *Tail;

inline char Getchar() {

    if(Head == Tail) {

        int l = fread(buffer, 1, BufferSize, stdin);

        Tail = (Head = buffer) + l;

    }

    return *Head++;

}

int read() {

    int x = 0, f = 1; char c = Getchar();

    while(!isdigit(c)){ if(c == '-') f = -1; c = Getchar(); }

    while(isdigit(c)){ x = x * 10 + c - '0'; c = Getchar(); }

    return x * f;

}

#define maxn 200010

#define oo 2147483647

int root, ToT, lc[maxn], rc[maxn];

bool Cur;

struct Node {

	int x[2], mx[2], mn[2], val, sum;

	bool operator == (const Node& t) const { return x[0] == t.x[0] && x[1] == t.x[1]; }

	bool operator < (const Node& t) const { return x[Cur] < t.x[Cur]; }

} nodes[maxn];

Node x, y;

void maintain(int o) {

	int l = lc[o], r = rc[o];

	for(int i = 0; i < 2; i++) {

		nodes[o].mx[i] = max(max(nodes[l].mx[i], nodes[r].mx[i]), nodes[o].x[i]);

		nodes[o].mn[i] = min(min(nodes[l].mn[i], nodes[r].mn[i]), nodes[o].x[i]);

	}

	nodes[o].sum = nodes[l].sum + nodes[r].sum + nodes[o].val;

//	printf("maintain(%d): %d %d %d %d %d\n", o, nodes[o].sum, nodes[o].mx[0], nodes[o].mn[0], nodes[o].mx[1], nodes[o].mn[1]);

	return ;

}

void add(int& o, bool cur) {

	if(!o){ nodes[o = ++ToT] = x; return maintain(o); }

	if(nodes[o] == x){ nodes[o].val += x.val; nodes[o].sum += x.val; return maintain(o); }

	add(x.x[cur] < nodes[o].x[cur] ? lc[o] : rc[o], cur ^ 1);

	return maintain(o);

}

void build(int& o, int L, int R, int cur) {

	if(L > R){ o = 0; return ; }

	int M = L + R >> 1; o = M;

	Cur = cur; nth_element(nodes + L, nodes + M, nodes + R + 1);

	build(lc[o], L, M - 1, cur ^ 1); build(rc[o], M + 1, R, cur ^ 1);

	return maintain(o);

}

bool all(int o) { return x.x[0] <= nodes[o].mn[0] && nodes[o].mx[0] <= y.x[0] && x.x[1] <= nodes[o].mn[1] && nodes[o].mx[1] <= y.x[1]; }

bool has(int o) { return !(nodes[o].mx[0] < x.x[0] || nodes[o].mn[0] > y.x[0] || nodes[o].mx[1] < x.x[1] || nodes[o].mn[1] > y.x[1]); }

int query(int o) {

	if(!o) return 0;

	int ans = 0;

	if(all(lc[o])) ans += nodes[lc[o]].sum;

	else if(has(lc[o])) ans += query(lc[o]);

	if(all(rc[o])) ans += nodes[rc[o]].sum;

	else if(has(rc[o])) ans += query(rc[o]);

	int nx = nodes[o].x[0], ny = nodes[o].x[1];

	if(x.x[0] <= nx && nx <= y.x[0] && x.x[1] <= ny && ny <= y.x[1]) ans += nodes[o].val;

	return ans;

}

int main() {

	nodes[0].mx[0] = nodes[0].mx[1] = -oo;

	nodes[0].mn[0] = nodes[0].mn[1] = oo;

	nodes[0].val = nodes[0].sum = 0;

	int n = read(), lastans = 0;

	n = read();

	while(n < 3) {

		if(n == 1) {

			x.x[0] = read() ^ lastans; x.x[1] = read() ^ lastans; x.val = read() ^ lastans;

			add(root, 1);

			if(ToT % 1000 == 0) build(root, 1, ToT, 1);

		}

		if(n == 2) {

			x.x[0] = read() ^ lastans; x.x[1] = read() ^ lastans; y.x[0] = read() ^ lastans; y.x[1] = read() ^ lastans;

			printf("%d\n", lastans = query(root));

		}

		n = read();

	}

	return 0;

}





其实这道题我第一眼看上去是树套树删边题，结果一看内存限制20M。。。还好没写树套树= =
												


											
[BZOJ2683][BZOJ4066]简单题的更多相关文章
	

    
								
  1. 【BZOJ1176】[Balkan2007]Mokia/【BZOJ2683】简单题 cdq分治
		
    [BZOJ1176][Balkan2007]Mokia Description 维护一个W*W的矩阵,初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=1600 ...
    
		
    2. 
						
  2. bzoj4066: 简单题 K-Dtree
		
    bzoj4066: 简单题 链接 bzoj 思路 强制在线.k-dtree. 卡常啊.空间开1e6就T了. 代码 #include <bits/stdc++.h> #define my_m ...
    
		
    3. 
						
  3. Bzoj4066 简单题
		
    Time Limit: 50 Sec  Memory Limit: 20 MBSubmit: 2185  Solved: 581 Description 你有一个N*N的棋盘,每个格子内有一个整数,初 ...
    
		
    4. 
						
  4. 【BZOJ2683】简单题 [分治][树状数组]
		
    简单题 Time Limit: 50 Sec  Memory Limit: 128 MB[Submit][Status][Discuss] Description 你有一个N*N的棋盘,每个格子内有一 ...
    
		
    5. 
						
  5. BZOJ4066 简单题（KD-Tree）
		
    板子题. #include<iostream> #include<cstdio> #include<cmath> #include<cstdlib> # ...
    
		
    6. 
						
  6. 【bzoj1176】[Balkan2007]Mokia/【bzoj2683】简单题  CDQ分治+树状数组
		
    bzoj1176 题目描述 维护一个W*W的矩阵,初始值均为S(题目描述有误,这里的S没有任何作用!).每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数 ...
    
		
    7. 
						
  7. 【BZOJ2683】简单题
		
    cdq分治妙啊 (被改过题面的)原题: dydxh所出的题目是这样的:有一个N*N矩阵,给出一系列的修改和询问,修改是这样的:将(x,y)中的数字加上k,而询问是这样的:求(x1,y1)到(x2,y2 ...
    
		
    8. 
						
  8. 【kd-tree】bzoj4066 简单题
		
    同p1176. #include<cstdio> #include<cmath> #include<algorithm> using namespace std;  ...
    
		
    9. 
						
  9. BZOJ4066:简单题(K-D Tree)
		
    Description 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作:   命令 参数限制 内容 1 x y A 1<=x,y<=N,A是正整数  ...
    
		
    10. 
		

	


随机推荐
	

    
									
  1. ModelProxy 前端接口配置建模框架
			
    ModelProxy    轻量级的接口配置建模框架(1) 先看一下这个博客说明为什么需要用ModelProxy的前端轻量级的框架吧:  http://developer.51cto.com/art/ ...
    
			
    2. 
						
  2. 四则运算 Day1
			
    需求分析 1.N(Need,需求) 该软件用户只有一类角色:小学生 要求能根据用户选择的算法(+-/),产生随机数(在0--10之间的整数)进行该算法的计算,程序能接收用户输入的答案,如果用户输入的为 ...
    
			
    3. 
						
  3. WebView加载本地html、js文件常见问题及解决办法
			
    声明:基于android studio平台,php语言搭建服务器 目录: 一.JavaScript脚本语言没有反应 二.alert无法弹出 三.html页面之间不能跳转 四.屏幕缩放没有达到预期效果  ...
    
			
    4. 
						
  4. DELL R710服务器做RAID5磁盘阵列图文教程
			
    本文转载于:http://www.jb51.net/article/53707.htm,只为做笔记使用 同时我们还可以选择一块硬盘做热备盘,就是说当配的raid5中有一个硬盘坏了的时候,它立马顶上,然 ...
    
			
    5. 
						
  5. [Asp.net mvc] 在Asp.net mvc 中使用MiniProfiler
			
    MiniProfiler是Stack Overflow团队设计的一款性能分析的小程序.可以对一个页面本身,及该页面通过直接引用.Ajax.Iframe形式访问的其它页面进行监控,监控内容包括数据库内容 ...
    
			
    6. 
						
  6. ACM算法总结及刷题参考
			
    参考:http://bbs.byr.cn/#!article/ACM_ICPC/11777 OJ上的一些水题(可用来练手和增加自信)(poj3299,poj2159,poj2739,poj1083,p ...
    
			
    7. 
						
  7. Cocos2d-X3.0 刨根问底（五）----- Node类及显示对象列表源码分析
			
    上一章 我们分析了Cocos2d-x的内存管理,主要解剖了 Ref.PoolManager.AutoreleasePool这三个类,了解了对象是如何自动释放的机制.之前有一个类 Node经常出现在各种 ...
    
			
    8. 
						
  8. BZOJ-1045 糖果传递   数学+递推
			
    1045: [HAOI2008] 糖果传递 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 2975 Solved: 1327 [Submit][Sta ...
    
			
    9. 
						
  9. RegexBuddy正则表达式工具
			
    RegexBuddy非常的好用,而且还能生成.net的代码. 我们在使用正则匹配时,毕竟.net提供的方法中,对于多行匹配就不能用单纯的正则去实现,而我们需要把它转换成相应的类库方法进行实现. 那么R ...
    
			
    10. 
						
  10. Java输入、输入、IO流 类层次关系梳理
			
    本文主要关注在Java编程中涉及到的IO相关的类库.方法.以及对各个层次(抽线.接口继承)的流之间的关系进行梳理 相关学习资料 http://baike.baidu.com/view/1007958. ...
    
			
    11.