【51Nod 1559】 车和矩形 （线段树）

题意
有一个 $n*m$ 的棋盘，其中有 $k$ 个点有棋子， $q$ 次个矩形。若矩形中每行都有一个棋子或没列都有一个棋子，则为被保护矩形，输出 “YES” ；否则 “NO”

题解
可以把所有点、矩形排序，然后扫描线方法判断
线段树记录每列最靠右的点的横坐标
对于一次询问，如果横坐标所在区间纵坐标最小值大于矩形左边界，那么这个矩形能被保护

这样操作可以把每行都有一个的情况解决，每列呢？
其实可以把行列互换就行了
注意：互换完n、m也要换…wa的贼惨

代码

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 200010
#define inf 0x3f3f3f3f
int n,m,k,q,tree[N<<1],flag[N];
struct node1{int x,y;}a[N];
struct node2{int x1,x2,y1,y2,id;}b[N];
inline bool cmp1(node1 x,node1 y){return x.y<y.y;}
inline bool cmp2(node2 x,node2 y){return x.y2<y.y2;}
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
inline void insert(int v,int l,int r,int x,int y){
    if(l==r){
        tree[v]=y;
        return ;
    }int mid=l+r>>1;
    if(x<=mid) insert(v<<1,l,mid,x,y);
    else insert(v<<1|1,mid+1,r,x,y);
    tree[v]=min(tree[v<<1],tree[v<<1|1]);
}
inline int query(int v,int l,int r,int x,int y){
    if(x<=l && r<=y) return tree[v];
    int mid=l+r>>1,s=inf;
    if(x<=mid) s=min(s,query(v<<1,l,mid,x,y));
    if(mid<y) s=min(s,query(v<<1|1,mid+1,r,x,y));
    return s;
}
void solve(int x){
    sort(a+1,a+k+1,cmp1);
    sort(b+1,b+q+1,cmp2);
    memset(tree,0,sizeof(tree));
    int i=1,j=1;
    while(j<=q){
        while(i<=k && a[i].y<=b[j].y2) insert(1,1,x,a[i].x,a[i].y),i++;
        if(query(1,1,x,b[j].x1,b[j].x2)>=b[j].y1) flag[b[j].id]=true;
        j++;
    }
}
int main(){
    freopen("a.in","r",stdin);
    n=read();m=read();k=read();q=read();
    for(int i=1;i<=k;i++) a[i].x=read(),a[i].y=read();
    for(int i=1;i<=q;i++){
        b[i].x1=read();b[i].y1=read();b[i].x2=read();b[i].y2=read();
        b[i].id=i;
    }
    solve(n);
//  for(int i=1;i<=k;i++) printf("%d %d %d %d %d\n",b[i].id,b[i].x1,b[i].y1,b[i].x2,b[i].y2);

    for(int i=1;i<=k;i++) swap(a[i].x,a[i].y);
    for(int i=1;i<=q;i++){
        swap(b[i].x1,b[i].y1);
        swap(b[i].x2,b[i].y2);
    }
    solve(m);
//  for(int i=1;i<=k;i++) printf("%d %d %d %d %d\n",b[i].id,b[i].x1,b[i].y1,b[i].x2,b[i].y2);

    for(int i=1;i<=q;i++) if(flag[i]) puts("YES");else puts("NO");
    return 0;
}