Cod sursa(job #680210)

Utilizator dushmiMihai-Alexandru Dusmanu dushmi Data 14 februarie 2012 12:25:39
Problema Zoo Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.8 kb
#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;
}