[bzoj4237]稻草人

Posted orzzz

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[bzoj4237]稻草人相关的知识,希望对你有一定的参考价值。

1A了好高兴。。cdq分治的思想。感觉难点在单调栈。


Description 

JOI村有一片荒地,上面竖着$N$个稻草人,村民们每年多次在稻草人们的周围举行祭典。
有一次,JOI村的村长听到了稻草人们的启示,计划在荒地中开垦一片田地。和启示中的一样,田地需要满足以下条件:
田地的形状是边平行于坐标轴的长方形;
左下角和右上角各有一个稻草人;
田地的内部(不包括边界)没有稻草人。

给出每个稻草人的坐标,请你求出有多少遵从启示的田地的个数

Input

第一行一个正整数N,代表稻草人的个数
接下来N行,第i行($1<=i<=N$)包含2个由空格分隔的整数$Xi$和$Yi$,表示第$i$个稻草人的坐标

Output

输出一行一个正整数,代表遵从启示的田地的个数

Sample Input

4
0 0
2 2
3 4
4 3

Sample Output

3

HINT

所有满足要求的田地由下图所示:

 技术分享图片

$1<=N<=2*10^5$

$0<=Xi<=10^9(1<=i<=N)$

$0<=Yi<=10^9(1<=i<=N)$

$Xi(1<=i<=N)$互不相同。

$Yi(1<=i<=N)$互不相同。


%出题人!!

首先分治处理:在区间$[l,r]$内,计算左下角的$x$在$[l,mid]$,右上角的$x$在$[mid+1,r]$的数量,递归下去。不重不漏。

那么就先将所有点按$x$排序,没啥说的。

那么在区间$[l,r]$内,分为左右两边$[l,mid],[mid+1,r]$:

左右各维护一个单调栈,左边称为$s1$,是单调递减,右边叫$s2$,是单调递增的。为啥呢?

考虑枚举右边的点。每次按$y$从小到大将右边的点加入$s2$。

右边的一个点$p$与左边的一个点$q$不能遵从启示(误)的情况有三种:

1.$p_y>q_y$  2.$p,q$之间有左边的点  3.$p,q$之间有右边的点

于是我们知道对于左边的点,它左下方的点不会造成影响。类似的,对于右边的点,其右下方也不会造成影响。

所以按y从小到大将右边点顺序加入$s2$,同时将栈内比他的x大的弹出$s2$。因为这些点造成情况3时,当前点也可以造成。

然后将所有y小于当前点的左侧点加入$s1$。每加入一个点,将栈内比他的$x$小的弹出$s1$。因为这些点不可能再与右侧点形成矩形了——他们会被加入的这个点阻碍形成情况2。

于是每次查找比加入$s2$前,$s2$的栈顶的$y$大的,$s1$中的点有多少个。二分就行了。否则会出现情况3啊。

说的很麻烦,看代码吧qwq。跑得巨慢 =_=。下面我造的一组样例可以手玩一下。

#include<bits/stdc++.h>
using namespace std;
const int N=200010;
inline int read(){
    int r=0,c=getchar();
    while(!isdigit(c))c=getchar();
    while(isdigit(c))
    r=r*10+c-0,c=getchar();
    return r;
}
struct qwq{
    int x,y;
}a[N];
bool cmpx(qwq p,qwq q){
    return p.x<q.x;
}
bool cmpy(qwq p,qwq q){
    return p.y<q.y;
}
struct Stack{
    int t;qwq s[N];
    Stack(){
        t=0;memset(s,0,sizeof s);
    }
    int size(){return t;}
    bool empty(){return !t;}
    void pop(){t--;}
    void push(qwq x){s[++t]=x;}
    qwq top(){return s[t];}
    int find(int x){
        int l=0,r=t;
        while(r-l>1){
            int mid=l+r>>1;
            if(s[mid].y<x)l=mid;
            else r=mid;
        }
        if(s[r].y<x)return r;
        else return l;
    }
    void clear(){while(!empty())pop();}
}sl,sr;
long long ans;
void cdq(int l,int r){
    if(l==r)return;
    int mid=l+r>>1;
    sort(a+l,a+mid+1,cmpy);sort(a+mid+1,a+r+1,cmpy);
    sl.clear();sr.clear();
    for(int i=mid+1,j=l;i<=r;i++){
        while(!sr.empty()&&sr.top().x>a[i].x)sr.pop();
        while(j<=mid&&a[j].y<a[i].y){
            while(!sl.empty()&&sl.top().x<a[j].x)sl.pop();
            sl.push(a[j]);j++;
        }
        int t=sl.size()-sl.find(sr.top().y);
        ans+=1ll*t;
        sr.push(a[i]);
    }
    sort(a+l,a+r+1,cmpx);
    cdq(l,mid);cdq(mid+1,r);
}
int main(){
    int n=read();
    for(int i=1;i<=n;i++)
    a[i].x=read(),a[i].y=read();
    sort(a+1,a+n+1,cmpx);
    cdq(1,n);
    printf("%lld\n",ans);
}
/*
7
3 1
7 2
4 3
1 4
6 5
2 6
5 7
*/

 






以上是关于[bzoj4237]稻草人的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ4237 稻草人

BZOJ#4237. 稻草人

BZOJ4237 稻草人

bzoj4237稻草人

bzoj 4237: 稻草人 -- CDQ分治

bzoj 4237 稻草人