HYSBZ4066-简单题

简单题

Time Limit: 50 Sec  Memory Limit: 20 MB
Submit: 3360  Solved: 880
[ Submit][ Status][ Discuss]

Description

你有一个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	无	终止程序

Input

输入文件第一行一个正整数N。

接下来每行一个操作。每条命令除第一个数字之外，

均要异或上一次输出的答案last_ans，初始时last_ans=0。

Output

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

Sample Input

4
1 2 3 3
2 1 1 3 3
1 1 1 1
2 1 1 0 7
3

Sample Output

3
5

HINT

数据规模和约定

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

样例解释见OJ2683

新加数据一组，但未重测----2015.05.24

Source

By wjy1998

解题思路：kd-tree

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <cctype>
#include <set>
#include <stack>
#include <queue>
#include <vector>
#include <bitset>
#include <functional>

using namespace std;

#define LL long long
const int INF = 0x3f3f3f3f;
const int N = 200000 + 5;
const int demension = 2;//二维

struct node
{
	int pos[demension];
	int ma[demension], mi[demension];
	int l, r, val, sum;
}a[N], x;
int cmpDem;//以第cmpDem维作比较
int root, n, ans, op, m;
int x2, y2, x3, y3;

bool cmp(const node &a, const node&b)
{
	if (a.pos[cmpDem] != b.pos[cmpDem]) return a.pos[cmpDem] < b.pos[cmpDem];
	else return a.pos[!cmpDem] < b.pos[!cmpDem];
}

void Merge(int k)
{
	for (int i = 0; i < demension; i++)
	{
		if (a[k].l)
		{
			a[k].ma[i] = max(a[k].ma[i], a[a[k].l].ma[i]);
			a[k].mi[i] = min(a[k].mi[i], a[a[k].l].mi[i]);
		}
		if (a[k].r)
		{
			a[k].ma[i] = max(a[k].ma[i], a[a[k].r].ma[i]);
			a[k].mi[i] = min(a[k].mi[i], a[a[k].r].mi[i]);
		}
	}
	if (a[k].l) a[k].sum += a[a[k].l].sum;
	if (a[k].r) a[k].sum += a[a[k].r].sum;
}

void Merge(int x, int y)
{
	for (int i = 0; i < demension; i++)
	{
		a[x].ma[i] = max(a[x].ma[i], a[y].ma[i]);
		a[x].mi[i] = min(a[x].mi[i], a[y].mi[i]);
	}
	a[x].sum += a[y].sum;
}

int build(int l, int r, int k)
{
	if (l > r) return 0;
	int mid = (l + r) / 2;
	//以第mid个元素为中心排序
	cmpDem = k;
	nth_element(a + l, a + mid, a + r + 1, cmp);
	for (int i = 0; i < demension; i++) a[mid].ma[i] = a[mid].mi[i] = a[mid].pos[i];
	a[mid].sum = a[mid].val;
	//左右子树
	a[mid].l = build(l, mid - 1, k ^ 1);
	a[mid].r = build(mid + 1, r, k ^ 1);
	Merge(mid);
	return mid;
}

int update(int k, int p)
{
	if (!k) return n;
	Merge(k, n);
	if (a[k].pos[p] > a[n].pos[p]) a[k].l = update(a[k].l, p ^ 1);
	else a[k].r = update(a[k].r, p ^ 1);
	return k;
}

int check(int k)
{
	if (x3 >= a[k].ma[0] && x2 <= a[k].mi[0] && y2 <= a[k].mi[1] && y3 >= a[k].ma[1]) return 1;
	else if (x2 > a[k].ma[0] || a[k].mi[0] > x3 || a[k].ma[1]<y2 || a[k].mi[1]>y3) return -1;
	else return 0;
}

int check1(int k)
{
	if (a[k].pos[0] >= x2&&a[k].pos[0] <= x3&&a[k].pos[1] >= y2&&a[k].pos[1] <= y3) return 1;
	else return 0;
}

void query(int k)
{
	if (check1(k)) ans += a[k].val;
	if (a[k].l)
	{
		int temp = check(a[k].l);
		if (temp == 1) ans += a[a[k].l].sum;
		else if (!temp) query(a[k].l);
	}
	if (a[k].r)
	{
		int temp = check(a[k].r);
		if (temp == 1) ans += a[a[k].r].sum;
		else if (!temp)query(a[k].r);
	}
}

int main()
{
	scanf("%d", &m);
	while (1)
	{
		scanf("%d", &op);
		if (op == 3) break;
		else if (op == 1)
		{
			n++;
			scanf("%d%d%d", &a[n].pos[0], &a[n].pos[1], &a[n].val);
			a[n].pos[0] ^= ans, a[n].pos[1] ^= ans, a[n].val ^= ans;
			a[n].l = a[n].r = 0;
			a[n].ma[0] = a[n].mi[0] = a[n].pos[0];
			a[n].ma[1] = a[n].mi[1] = a[n].pos[1];
			a[n].sum = a[n].val;
			root = update(root, 0);
			if (n % 10000 == 0) root = build(1, n, 0);
		}
		else
		{
			scanf("%d%d%d%d", &x2, &y2, &x3, &y3);
			x2 ^= ans, y2 ^= ans, x3 ^= ans, y3 ^= ans;
			ans = 0;
			query(root);
			printf("%d\n", ans);
		}
	}
	return 0;
}