18th 省赛 B. Restore Atlantis

题意：
有$2000\ast 2000$ 的二维平面，给你 $100000$ 个矩形，有 $1e5$个询问，每次询问 给你 $l和r$ ，问去掉这个区间的矩形之后，平面中被矩形覆盖的点数。

思路1：
用扫描线思想，并且用 优先队列的线段树，处理出 二维平面上 所有点，被覆盖的最左边的矩形 和 最右边的矩形。（听说用set会T，可能优先队列的常熟小吧）
因为题目的时限只有1s，考虑离线的算法，用树状数组来维护。将所有的询问离线，并按照 $r$升序排序。
当$i$等于$r$的时候，将对应的$l$加一，然后区间询问的时候，直接查树状数组对应的区间，即此时覆盖这个点的矩形都包含在 $l$和$r$ 中。感性理解即可。

代码1：

#include<bits/stdc++.h>
#define rep(i,n,m) for(ll i=n;i<=m;i++)
#define repp(i,n,m) for(ll i=n;i>=m;i--)
#define mem(a) memset(a,0,sizeof(a))
#define debug1(a) cout<<a<<endl
#define debug2(a,b) cout<<a<<" "<<b<<endl
#define debug3(a,b,c) cout<<a<<" "<<b<<" "<<c<<endl;
#define whlie while
#define itn int
#define icn cin
#define inc cin
#define mian main
typedef int ll;
using namespace std;

int n,m;
struct node{
	priority_queue<int>q1[2];
	priority_queue<int,vector<int>,greater<int>>q2[2];
}tr[8060];
int L[2005][2005],R[2005][2005];
struct no{
	int op,l,r,id;
};
vector<no>v[2005];
void update(int p,int l,int r,int x,int y,int op,int w){
	if(x<=l&&r<=y){
		tr[p].q1[op].push(w);
		tr[p].q2[op].push(w);
		return;
	}
	int mid=l+r>>1;
	if(x<=mid)update(2*p,l,mid,x,y,op,w);
	if(mid<y)update(2*p+1,mid+1,r,x,y,op,w);
}
void dfs(int p,int l,int r,int mx,int mi,int id){
	while(!tr[p].q1[0].empty()&&!tr[p].q1[1].empty()&&tr[p].q1[0].top()==tr[p].q1[1].top()){
		tr[p].q1[0].pop();tr[p].q1[1].pop();
	}
	while(!tr[p].q2[0].empty()&&!tr[p].q2[1].empty()&&tr[p].q2[0].top()==tr[p].q2[1].top()){
		tr[p].q2[0].pop();tr[p].q2[1].pop();
	}
	if(!tr[p].q1[1].empty())mx=max(mx,tr[p].q1[1].top());
	if(!tr[p].q2[1].empty())mi=min(mi,tr[p].q2[1].top());
	if(l==r){
		L[id][l]=mi;
		R[id][l]=mx;
		return;
	}
	int mid=l+r>>1;
	dfs(2*p,l,mid,mx,mi,id); 
	dfs(2*p+1,mid+1,r,mx,mi,id); 
}
vector<int>vv[100005];
struct nod{
	int x,id;
};
int ans[100050];
vector<nod>v1[100005];
int sum[100005];
int lowbit(int x){
	return x&(-x);
}
void add(int x,int w){
	while(x<=100001){
		sum[x]+=w;
		x+=lowbit(x);
	}
}
int ask(int x){
	int ans=0;
	while(x){
		ans+=sum[x];
		x-=lowbit(x);
	}
	return ans;
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		int a,b,c,d;
		scanf("%d%d%d%d",&a,&b,&c,&d);
		a++,b++,c++,d++;
		v[a].push_back({1,b,d-1,i});
		v[c].push_back({0,b,d-1,i});
	}
	int tot=0;
	for(int i=1;i<=2001;i++){
		for(auto to:v[i]){
			update(1,1,2001,to.l,to.r,to.op,to.id);
		}
		dfs(1,1,2001,0,n+1,i);
		for(int j=1;j<=2001;j++){
			if(R[i][j]){
				vv[R[i][j]].push_back(L[i][j]);
				tot++;
			}
		}
	}
	for(int i=1;i<=m;i++){
		int a,b;
		scanf("%d%d",&a,&b);
		v1[b].push_back({a,i});
	}
	for(int i=1;i<=n;i++){
		for(auto to:vv[i]){
			add(to,1);
		}
		for(auto to:v1[i]){
			ans[to.id]=tot-ask(i)+ask(to.x-1);
		}
	}
	for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
	return 0;
}