是一道比較神仙的數據結構（我自己肯定想不出來）
這個查詢其實是比較難處理的，對於一段詢問區間[l, r]，每一個矩形和其他的都是可能相互聯系的
但是可以注意到區間的范圍非常小，是[0, 2000]之間，所以可以考慮對每一個1×1大小的方格進行處理，可以發現方格的數量是2000*2000個，是一個可以接受的范圍，那么我們現在來考慮對於每次詢問的[l, r]，每個小方格是否存在一個不在[l, r]的區間的矩形能包含它。
由於每次詢問的區間都是連續的，所以只有當前這個詢問區間完全包含所有出現了這個小方格的矩形的時候，這個小方格不被任何一個矩形包含，所以我們只需要知道每個小方格出現在的矩形中編號最小的和編號最大的就可以了。
求解這個的部分可以用掃描線+線段樹維護優先隊列解決，對於所有矩形，掃描線按照x從小到大掃，對每一個矩形的左邊，區間修改（就是往優先隊列里丟這個編號），如果掃到了矩形的右邊，就打上標記。
然后查詢所有的點，得到答案，其實在這里，這個優先隊列相當於lazy標記，所以在查詢的時候需要下傳，而我們之所以要用一個優先隊列就是為了把打上標記的點去掉，所以取出數的時候要找到第一個沒有標記的點。
現在我們得到了每個點的最小出現編號和最大出現編號，問題就又轉換成了給定\(C^2\)條線段，每次詢問\([l, r]\)，求多少個線段不被\([l, r]\)完全包含。
這個方法就比較多了，在線的話可以用什么分塊，線段樹，主席樹來維護，當然這道題是離線，所以可以用掃描線+樹狀數組簡單維護，大概就是先按照y的坐標排序，然后按照從小到大把x加進樹狀數組里，加的過程中順便算出來每個區間的答案
代碼參考了一下claris，寫的太優雅了%%%
#include <bits/stdc++.h>
using namespace std;
const int N = 100010;
bool vis[N];
struct node {
int x, y1, y2, id, flag;
inline bool operator < (const node & o) const {
return x < o.x;
}
}v[N << 1];
priority_queue < int > maxn[N], minn[N];
void update(int o, int l, int r, int L, int R, int x) {
if (r < L || l > R) return;
if (L <= l && r <= R) {
maxn[o].push(x);
minn[o].push(-x);
return ;
}
int mid = (l + r) >> 1;
update(o << 1, l, mid, L, R, x);
update(o << 1 | 1, mid + 1, r, L, R, x);
}
vector < pair<int, int> > seg;
int tot;
void query(int o, int l, int r, int x, int y) {
while (!maxn[o].empty() && vis[maxn[o].top()]) maxn[o].pop();
while (!minn[o].empty() && vis[-minn[o].top()]) minn[o].pop();
if (!maxn[o].empty()) x = max(x, maxn[o].top()), y = max(y, minn[o].top());
if (l == r) {
if (x) seg.emplace_back(x, -y), tot++;
return ;
}
int mid = (l + r) >> 1;
query(o << 1, l, mid, x, y);
query(o << 1 | 1, mid + 1, r, x, y);
}
int n, q, c[N];
inline int lowbit(int x) { return x & (-x); }
inline void add(int x) { while (x <= n) {c[x]++; x += lowbit(x);} }
inline int get(int x) { int ans = 0; while (x) { ans += c[x]; x -= lowbit(x);} return ans; }
vector < tuple<int, int, int> > qs;
int ans[N];
int main() {
scanf("%d%d", &n, &q);
for (int i = 1; i <= n; i++) {
int xa, ya, xb, yb;
scanf("%d%d%d%d", &xa, &ya, &xb, &yb);
xa++; ya++; xb++; yb++;
v[i] = node{xa, ya, yb, i, 1};
v[i + n] = node{xb, ya, yb, i, -1};
}
sort(v + 1, v + 2 * n + 1);
int p = 1;
for (int i = 1; i <= 2000; i++) {
while (v[p].x == i) {
if (v[p].flag == 1)
update(1, 1, 2000, v[p].y1, v[p].y2 - 1, v[p].id);
else vis[v[p].id] = true;
p++;
}
query(1, 1, 2000, 0, -n);
}
sort(seg.begin(), seg.end());
for (int i = 1; i <= q; i++) {
int x, y;
scanf("%d%d", &x, &y);
qs.emplace_back(y, x, i);
}
sort(qs.begin(), qs.end());
int cnt = 0; p = 0;
for (auto [y, x, id]:qs) {
while (seg[p].first <= y && p < tot) add(seg[p++].second), cnt++;
ans[id] = tot - cnt + get(x - 1);
}
for (int i = 1; i <= q; i++) printf("%d\n", ans[i]);
return 0;
}