Problem:
宝探し

太郎君はある広場にお宝を探しにやってきました。この広場にはたくさんのお宝が埋められていますが、太郎君は最新の機械を持っているので、どこにお宝が埋まっているかをすべて知っています。広場はとても広いので太郎君は領域を決めてお宝を探すことにしましたが、お宝はたくさんあるためどのお宝がその領域の中にあるかすぐにわかりません。そこで太郎君はその領域の中にあるお宝の数を数えることにしました。

Input

n m
x1 y1
x2 y2
...
xn yn
x11 y11 x12 y12
x21 y21 x22 y22
...
xm1 ym1 xm2 ym2

nは広場に埋まっているお宝の数を表す

mは調べる領域の数を表す

2行目からn+1行目はそれぞれのお宝が埋まっている座標を表す

n+2行目からn+m+1行目はそれぞれの調べる領域を表す

x軸の正の方向が東、y軸の正の方向が北を表す

それぞれの領域は長方形であり、xi1とyi1は長方形の南西の頂点の座標、xi2とyi2は長方形の北東の頂点の座標を表す

Constraints

1 ≤ n ≤ 5000
1 ≤ m ≤ 5×105
|xi|, |yi| ≤ 109 (1 ≤ i ≤ n)
|xi1|, |yi1|, |xi2|, |yi2| ≤ 109 (1 ≤ i ≤ m)
xi1 ≤ xi2, yi1 ≤ yi2 (1 ≤ i ≤ m)

すべての入力は整数で与えられる

Output

C1
C2
...
Cm

それぞれの領域に含まれるお宝の数を各行に出力せよ

Sample Input 1

3 1
1 1
2 4
5 3
0 0 5 5

Output for the Sample Input 1

3

領域の境界線上にあるお宝も領域に含まれるものとみなす

Sample Input 2

4 2
-1 1
0 3
4 0
2 1
-3 1 5 1
4 0 4 0

Output for the Sample Input 2

2
1

領域は直線や点のようになる場合もある

Sample Input 3

2 3
0 0
0 0
-1 -1 1 1
0 0 2 2
1 1 4 4

Output for the Sample Input 3

2
2
0

同じ場所に複数のお宝が存在する場合もある

Sample Input 4

5 5
10 5
-3 -8
2 11
6 0
-1 3
-3 1 3 13
-1 -1 9 5
-3 -8 10 11
0 0 5 5
-10 -9 15 10

Output for the Sample Input 4

2
2
5
0
4

Solution in C++ (judge verdict: Accepted)

#include <iostream>
#include <vector>
#include <algorithm>
#define inf (1e9+9)

using namespace std;

int n, m;
int x[5005], y[5005];
vector<int> compX, compY;
int dif[5005][5005], sum[5005][5005];

int main(void)
{
	cin >> n >> m;
	for(int i = 1; i <= n; i++) cin >> x[i] >> y[i];
	
	compX.push_back(-inf);
	for(int i = 1; i <= n; i++) compX.push_back(x[i]);
	sort(compX.begin(), compX.end());
	compX.erase(unique(compX.begin(), compX.end()), compX.end());
	for(int i = 1; i <= n; i++) x[i] = lower_bound(compX.begin(), compX.end(), x[i]) - compX.begin();
	
	compY.push_back(-inf);
	for(int i = 1; i <= n; i++) compY.push_back(y[i]);
	sort(compY.begin(), compY.end());
	compY.erase(unique(compY.begin(), compY.end()), compY.end());
	for(int i = 1; i <= n; i++) y[i] = lower_bound(compY.begin(), compY.end(), y[i]) - compY.begin();
	
	for(int i = 1; i <= n; i++) dif[x[i]][y[i]]++;
	for(int x = 1; x < 5005; x++){
		for(int y = 1; y < 5005; y++){
			sum[x][y] = sum[x-1][y] + sum[x][y-1] - sum[x-1][y-1] + dif[x][y];
		}
	}
	
	//for(int i = 0; i < compY.size(); i++) cout << compY[i] << " "; cout << endl;
	
	int sx, sy, tx, ty;
	for(int q = 0; q < m; q++){
		cin >> sx >> sy >> tx >> ty;
		sx = lower_bound(compX.begin(), compX.end(), sx) - compX.begin();
		sy = lower_bound(compY.begin(), compY.end(), sy) - compY.begin();
		tx = upper_bound(compX.begin(), compX.end(), tx) - compX.begin() - 1;
		ty = upper_bound(compY.begin(), compY.end(), ty) - compY.begin() - 1;
		//cout << sx << " " << sy << " " << tx << " " << ty << endl;
		cout << sum[tx][ty] - sum[sx-1][ty] - sum[tx][sy-1] + sum[sx-1][sy-1] << endl;
	}
	
	return 0;
}