1790: [Ahoi2008]Rectangle 矩形藏宝地
Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 393 Solved: 129
[Submit][Status][Discuss]
Description
欢乐岛上最著名的游戏是一个寻宝游戏,小可可来到宝藏的埋藏地,这是一块开阔地,宝藏被分散的埋藏在这块地下,现在要做的是一件件的把宝藏挖出来。为了提示宝藏的埋藏点,游戏的主办方把这块开阔地当作第一象限,将所有可能埋藏宝藏的地方划成一个个矩形的土地,并把这些矩形土地的坐标都告诉了参赛者。挖宝的提示很简单,只要某一个矩阵土地至少被另外一个矩阵土地所包含,那么这个矩阵土地里肯定埋有宝藏。其实这些宝藏都是一些精美的纪念品,如果谁挖出来了纪念品就归谁了,小可可很想为这次旅程画上完美的句号,有你的帮助他信心十足,你只要告诉他.有多少个矩形土地里肯定埋有宝藏就行了。胜利就在眼前,加油吧!!
Input
第一行包含一个整数N(N≤200000),表示矩形的个数。接下来N行,每行用4个整数x1,y1,x2,y2,描述了一个矩形。其中(x1,y1)表示这个矩形左下角的坐标,(x2,y2)表示右上角的坐标,一个xi值或yi值最多出现一次.
Output
只包含一个整数,表示肯定埋有宝藏的矩形土地的个数。
Sample Input
3
0 0 5 5
1 2 3 4
2 1 4 3
Sample Output
2
HINT
100%的数据中,N<=200000
70%的数据申,N<=50000
30%的数据中,N<=5000
所有数据中,一个x值或Y值最多出现一次
每一个二维的平面可以用两个点来表示。
kd
树不能处理这种包含的问题。那就可以吧这个平面转化一下。
两个点表示的平面需要四个坐标,那么我们随便拿出三个来,比如说选
(x1,y1,x2)(x1<x2,y1<y2)
令
y2
表示这个点权值,那么每次就是相当于查询
(0 x1−1,0 y1−1,x2 max)
的最大值。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=200010;
int n,root,D,ans,lim,now,va;
struct S{
int l,r,d[3],mi[3],ma[3],v,maxn;
int &operator [](int x){
return d[x];
}
bool operator < (const S &x)const{
return d[D]<x.d[D];
}
}tr[N],p[N];
inline int IN(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){
if(ch=='-') f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x*f;
}
inline void update(int k){
int l=tr[k].l,r=tr[k].r,i;
tr[k].maxn=max(tr[k].v,max(tr[l].maxn,tr[r].maxn));
for(i=0;i<=2;++i){
tr[k].mi[i]=tr[k].ma[i]=tr[k][i];
if(l){
tr[k].mi[i]=min(tr[k].mi[i],tr[l].mi[i]);
tr[k].ma[i]=max(tr[k].ma[i],tr[l].ma[i]);
}
if(r){
tr[k].mi[i]=min(tr[k].mi[i],tr[r].mi[i]);
tr[k].ma[i]=max(tr[k].ma[i],tr[r].ma[i]);
}
}
}
#define mid (l+r)/2
inline int build(int l,int r,int flag){
if(l>r) return 0;
D=flag;nth_element(p+l,p+mid,p+r+1);
tr[mid]=p[mid];
tr[mid].l=build(l,mid-1,(flag+1)%3);
tr[mid].r=build(mid+1,r,(flag+1)%3);
update(mid);
return mid;
}
inline bool in(int k,int x1,int y1,int x2,int y2,int x3,int y3){
return tr[k].mi[0]>=x1&&tr[k].ma[0]<=y1&&tr[k].mi[1]>=x2&&tr[k].ma[1]<=y2&&tr[k].mi[2]>=x3&&tr[k].ma[2]<=y3;
}
inline bool out(int k,int x1,int y1,int x2,int y2,int x3,int y3){
return tr[k].mi[0]>y1||tr[k].ma[0]<x1||tr[k].mi[1]>y2||tr[k].ma[1]<x2||tr[k].mi[2]>y3||tr[k].ma[2]<x3;
}
inline bool in1(int k,int x1,int y1,int x2,int y2,int x3,int y3){
return tr[k][0]>=x1&&tr[k][0]<=y1&&tr[k][1]>=x2&&tr[k][1]<=y2&&tr[k][2]>=x3&&tr[k][2]<=y3;
}
inline void query(int k,int x1,int y1,int x2,int y2,int x3,int y3){
if(now>va) return ;
if(x1>y1||x2>y2||x3>y3||!k) return ;
if(in(k,x1,y1,x2,y2,x3,y3)){
now=max(now,tr[k].maxn);
return ;
}
if(out(k,x1,y1,x2,y2,x3,y3)) return ;
if(in1(k,x1,y1,x2,y2,x3,y3)) now=max(now,tr[k].v);
if(tr[tr[k].l].maxn>va) query(tr[k].l,x1,y1,x2,y2,x3,y3);
if(tr[tr[k].r].maxn>va) query(tr[k].r,x1,y1,x2,y2,x3,y3);
}
int main(){
int i,x1,x2,y1,y2,j;
n=IN();
for(i=1;i<=n;++i){
x1=IN();y1=IN();x2=IN();y2=IN();
lim=max(lim,x2);
p[i][0]=x1;p[i][1]=y1,p[i][2]=x2;
p[i].v=p[i].maxn=y2;
for(j=0;j<=2;++j)
p[i].mi[j]=p[i].ma[j]=p[i][j];
}
root=build(1,n,0);
for(i=1;i<=n;++i){
now=0;va=p[i].v;
query(root,0,p[i][0]-1,0,p[i][1]-1,p[i][2],lim);
ans+=(now>p[i].v);
}
printf("%d\n",ans);
}