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 x₁ y₁ x₂ y₂

1<=x₁<= x₂<=N

1<=y₁<= y₂<=N

输出x₁ y₁ x₂ y₂这个矩形内的数字和

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 longconst 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;}