Pagini recente » Cod sursa (job #1017945) | Cod sursa (job #664938) | Cod sursa (job #760088) | Cod sursa (job #782212) | Cod sursa (job #680210)
Cod sursa(job #680210)
#include <cstdio>
#include <cassert>
#include <algorithm>
#include <vector>
using namespace std;
const int N = 16500, M = 100500, INF = 2000000500;
struct point {
int x, y;
};
struct rectangle {
int x1, y1, x2, y2;
};
rectangle q[M];
point v[N], t[4 * M];
int n, m, newm, ans[4 * M], vy[3 * M], cory[3 * M], aib[3 * M];
void read() {
assert(scanf("%d", &n));
for (int i = 1; i <= n; ++i)
assert(scanf("%d%d", &v[i].x, &v[i].y));
assert(scanf("%d", &m));
for (int i = 1; i <= m; ++i)
assert(scanf("%d%d%d%d", &q[i].x1, &q[i].y1, &q[i].x2, &q[i].y2));
}
int cautb(int x, int *v) {
int i = 0, pas = 1 << 18;
for (i = 0; pas; pas >>= 1)
if (i + pas <= v[0] && v[i + pas] < x)
i += pas;
return i + 1;
}
void normalization() {
vy[++vy[0]] = -INF;
for (int i = 1; i <= n; ++i)
vy[++vy[0]] = v[i].y;
for (int i = 1; i <= m; ++i) {
vy[++vy[0]] = q[i].y1;
vy[++vy[0]] = q[i].y2;
}
vy[++vy[0]] = INF;
sort(vy + 1, vy + vy[0] + 1);
for (int i = 1; i <= vy[0]; ++i) {
cory[i] = cory[i - 1];
if (vy[i] != vy[i - 1])
++cory[i];
}
for (int i = 1; i <= n; ++i)
v[i].y = cautb(v[i].y, vy);
for (int i = 1; i <= m; ++i) {
q[i].y1 = cautb(q[i].y1, vy);
q[i].y2 = cautb(q[i].y2, vy);
}
}
bool comp(point A, point B) {
if (A.x < B.x)
return 1;
if (A.x > B.x)
return 0;
return A.y < B.y;
}
void init() {
sort(v + 1, v + n + 1, comp);
newm = 0;
for (int i = 1; i <= m; ++i) {
++newm;
t[newm].x = q[i].x2;
t[newm].y = q[i].y2;
++newm;
t[newm].x = q[i].x2;
t[newm].y = q[i].y1 - 1;
++newm;
t[newm].x = q[i].x1 - 1;
t[newm].y = q[i].y2;
++newm;
t[newm].x = q[i].x1 - 1;
t[newm].y = q[i].y1 - 1;
}
sort(t + 1, t + newm + 1, comp);
}
inline int bit(int x) {
return x & (-x);
}
inline int query(int poz) {
int rez = 0;
while (poz) {
rez += aib[poz];
poz -= bit(poz);
}
return rez;
}
inline void insert(int poz) {
while (poz < 3 * M) {
++aib[poz];
poz += bit(poz);
}
}
int cb(int x, int y) {
int pas = 1 << 18, i = 0;
for (i = 0; pas; pas >>= 1)
if (i + pas <= newm && t[i + pas].x <= x && t[i + pas].y < y)
i += pas;
return i + 1;
}
void solve() {
for (int i = 0, j = 0; i <= n || j <= newm;)
if (i > n)
ans[j] = query(t[j].y);
else if (j > newm)
insert(v[i].y);
else if (v[i].x <= t[j].x)
insert(v[i].y);
else
ans[j] = query(t[j].y);
for (int i = 1; i <= m; ++i)
printf("%d\n", ans[cb(q[i].x2, q[i].y2)] - ans[cb(q[i].x2, q[i].y1 - 1)] - ans[cb(q[i].x1 - 1, q[i].y2)] + ans[cb(q[i].x1 - 1, q[i].y1 - 1)]);
}
int main() {
assert(freopen("zoo.in", "r", stdin));
assert(freopen("zoo.out", "w", stdout));
read();
normalization();
init();
solve();
return 0;
}