Cod sursa(job #220415)

Utilizator savimSerban Andrei Stan savim Data 10 noiembrie 2008 19:45:47
Problema Zoo Scor 80
Compilator cpp Status done
Runda Arhiva de probleme Marime 3.09 kb
#include <stdio.h>
#include <string.h>
#include <algorithm>

using namespace std;

#define MAX_N 100000
#define MAX_L 417000
#define inf 2147000000

int n, m, i, j, k, x1, y1, x2, y2, p, q, sum, semn;
int aib[MAX_L];
int intreb[MAX_L][2], v[MAX_L];
int ordon[MAX_L], vord[MAX_L];
int qry[MAX_N][5];
char s[55];

struct punct {
	int x;
	int y;
};

punct a[MAX_L], b[MAX_L];

void cit() {
	freopen("zoo.in", "r", stdin);
	freopen("zoo.out", "w", stdout);
	
	//citirea parsata
	scanf("%d\n",&n);
	for (i = 1; i <= n; v[i] = i, i++) 
		scanf("%d %d",&a[i].x, &a[i].y);
	
	scanf("%d\n", &m);
	for (i = 1; i <= m; i++) {
		fgets(s,50,stdin);
		
		j = 0;
		semn = 1;x1 = 0;
		if (s[j] == '-') {semn = -1; j++;}
		for (j = j; s[j] != ' '; j++)
			x1 = x1 * 10 + s[j] - '0';
		x1 *= semn;
		
		j++;
		semn = 1; y1 = 0;
		if (s[j] == '-') {semn = -1; j++;}
		for (j = j; s[j] != ' '; j++)
			y1 = y1 * 10 + s[j] - '0';
		y1 *= semn;
		
		j++;
		semn = 1; x2 = 0;
		if (s[j] == '-') {semn = -1; j++;}
		for (j = j; s[j] != ' '; j++)
			x2 = x2 * 10 + s[j] - '0';
		x2 *= semn;
		
		j++;
		semn = 1; y2 = 0;
		if (s[j] == '-') {semn = -1; j++;}
		for (j = j; s[j] != '\n' && s[j] != '\0'; j++)
			y2 = y2 * 10 + s[j] - '0';
		y2 *= semn;

		a[++n].x = x1 - 1; a[n].y = y1 - 1; v[n] = n; intreb[n][0] = i;
		a[++n].x = x1 - 1; a[n].y = y2; v[n] = n; intreb[n][0] = i;
		a[++n].x = x2; a[n].y = y1 - 1; v[n] = n; intreb[n][0] = i;
		a[++n].x = x2; a[n].y = y2; v[n] = n; intreb[n][0] = i;
	}		
}

inline int cmp2(int a, int b) {
	return ordon[a] < ordon[b];
}

void normalizare () {
	p = 1; q = 1;
	for (i = 1; i <= n; i++) {
		b[i].x = p;
		if (a[v[i]].x > a[v[i - 1]].x) b[i].x = ++p;

		ordon[i]  = a[v[i]].y;
		vord[i] = i;
	}

	sort(vord + 1, vord + n + 1, cmp2);

	for (i = 1; i <= n; i++) {
		b[vord[i]].y = q;
		if (ordon[vord[i]] > ordon[vord[i - 1]]) b[vord[i]].y = ++q;
	}
}

void update(int k) {
	while (k <= q) {
	    aib[k]++;
		k += k ^ (k & (k - 1));
	}
}

int query(int k) {
	sum = 0;
	while (k) {
		sum += aib[k];
		k -= k ^ (k & (k - 1));
	}
	return sum;
}

void parcurg() {
	for (i = 1; i <= n; i++) {
		//update
		if (intreb[v[i]][0] == 0) update(b[i].y);
		
		//query
		if (intreb[v[i]][0] != 0) intreb[v[i]][1] = query(b[i].y);
	}
}

inline int cmp(int p, int q) {
    if ((a[p].x < a[q].x) || (a[p].x == a[q].x && a[p].y < a[q].y) ||   
        (a[p].x == a[q].x && a[p].y == a[q].y && intreb[p] < intreb[q]))  
       return true;  
    else return false;   
}

void solve() {
	//sortarea elementelor
	sort(v + 1, v + n + 1, cmp);
	
	//normalizarea elementelor
	normalizare();
	
	//parcurgerea elementelor si query - urile respective
	parcurg();
	
	//updatez valorile intrebarilor
	for (i = 1; i <= n; i++) 
		if (intreb[v[i]][0] != 0) {
			p = intreb[v[i]][0];
			qry[p][++qry[p][0]] = intreb[v[i]][1];
		}
}

void write() {
	for (i = 1; i <= m; i++) {
		sum = qry[i][4] - qry[i][2] - qry[i][3] + qry[i][1];
		printf("%d\n",sum);
	}
}

int main() {

	cit();
	solve();
	write();
	
	return 0;
}