题意:
有
2000
∗
2000
2000*2000
2000∗2000 的二维平面,给你
100000
100000
100000 个矩形,有
1
e
5
1e5
1e5个询问,每次询问 给你
l
和
r
l和r
l和r ,问去掉这个区间的矩形之后,平面中被矩形覆盖的点数。
思路1:
用扫描线思想,并且用 优先队列的线段树,处理出 二维平面上 所有点,被覆盖的最左边的矩形 和 最右边的矩形。(听说用set会T,可能优先队列的常熟小吧)
因为题目的时限只有1s,考虑离线的算法,用树状数组来维护。将所有的询问离线,并按照
r
r
r升序排序。
当
i
i
i等于
r
r
r的时候,将对应的
l
l
l加一,然后区间询问的时候,直接查树状数组对应的区间,即此时覆盖这个点的矩形都包含在
l
l
l和
r
r
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;
}