稻草人(bzoj 4237)

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)互不相同。
/*
  网上说这是一道cdq分治
  我们对于y坐标进行分治,对于某一段先按照x坐标排序,考虑上半部分对于下半部分的影响。 
  考虑上半部分某个点i,如果它能作为右上的点与下半部分的某个点形成长方形,那么必须保证在他们之间没有某个点j,使得x[j]<x[i]&&y[j]<y[i],画图可知,我们可以在上半部分维护一个y递增的单调栈,下半部分维护一个y递减的单调栈,这样能保证一一对应,然后二分求出第一个满足条件的坐标。 
*/
#include<cstdio>
#include<iostream>
#include<algorithm>
#define N 200010
#define lon long long
using namespace std;
int n,p[N],q[N];lon ans;
struct node{int x,y;}a[N],b[N];
bool cmpy(const node&s1,const node&s2){return s1.y<s2.y;}
int find(int x,int l,int r){
    while(l+1<r){
        int mid=l+r>>1;
        if(a[q[mid]].x<x) l=mid;
        else r=mid;
    }
    return l;
}
void solve(int l,int r){
    if(l==r)return;
    int mid=l+r>>1;
    solve(l,mid);solve(mid+1,r);
    int tp1=0,tp2=0,i,j=l,k;
    for(i=mid+1;i<=r;i++){
        while(tp1&&a[i].y<a[p[tp1]].y)tp1--;
        p[++tp1]=i;
        for(;a[j].x<a[i].x&&j<=mid;j++){
            while(tp2&&a[j].y>a[q[tp2]].y)tp2--;
            q[++tp2]=j;
        }
        ans+=tp2-find(a[p[tp1-1]].x,0,tp2+1);
    }
    j=l;k=mid+1;
    for(i=l;i<=r;i++)
        b[i]=((j<=mid&&a[j].x<a[k].x)||k>r)?a[j++]:a[k++];
    for(i=l;i<=r;i++)a[i]=b[i];
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d%d",&a[i].x,&a[i].y);
    sort(a+1,a+n+1,cmpy);
    solve(1,n);
    cout<<ans;
    return 0;
}
原文地址:https://www.cnblogs.com/harden/p/6421784.html