bzoj 4237 稻草人

http://www.elijahqi.win/archives/3461
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)互不相同。

Source
JOI 2013~2014 春季training合宿 竞技3 By PoPoQQQ

cdq分治

考虑没有中间那个不能有稻草人的限制直接cdq求偏序即可

考虑有了这样的限制之后我们仍然将y从小到大排序然后每次做的时候先c递归处理左右两边 然后考虑合并两边的答案

先将两边各自按照x排序 然后针对左边按照y维护递降单调队列

右边维护y递升单调队列 然后每次相当于是做右边的时候取单调队列里前一个的x然后去左边单调队列二分找到第一个大的即可 那么左边单调队列那个位置到队尾都是可以算答案的 具体画图感受

#include<cstdio>
#include<cctype>
#include<algorithm>
#define ll long long
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int N=2e5+10;
struct node{
    int x,y;
}a[N];
inline bool cmp1(const node &a,const node &b){return a.y<b.y;}
inline bool cmp2(const node &a,const node &b){return a.x<b.x;}
int n,q1[N],q2[N];ll ans;
inline void cdq(int l,int r){
    if(l==r) return;int mid=l+r>>1;
    cdq(l,mid);cdq(mid+1,r);int top1=0,top2=0,j=l;
    sort(a+l,a+mid+1,cmp2);sort(a+mid+1,a+r+1,cmp2);
    for (int i=mid+1;i<=r;++i){
        while(top1&&a[q1[top1]].y>a[i].y) --top1;
        q1[++top1]=i;
        while(a[j].x<a[i].x&&j<=mid){
            while(top2&&a[q2[top2]].y<a[j].y) --top2;
            q2[++top2]=j;++j;
        }int x=a[q1[top1-1]].x,l1=1,r1=top2,md=0;
        while(l1<=r1){
            int mid=l1+r1>>1;
            if (a[q2[mid]].x>x) md=mid,r1=mid-1;else l1=mid+1;
        }if(md) ans+=top2-md+1;
    }
}
int main(){
    freopen("bzoj4237.in","r",stdin);
    n=read();a[0].x=-1;
    for (int i=1;i<=n;++i) a[i].x=read(),a[i].y=read();
    sort(a+1,a+n+1,cmp1);cdq(1,n);
    printf("%lld\n",ans);
    return 0;
}