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

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

到现在才做稻草人的我...，其实这个题也不算太难，今天考试就是要用这个东西优化dp；

对于平面点对计数，我们考虑分治，那么我们假设按照y来分治，那么图分为上半部分和下半部分；

我们只考虑跨过分界线的，同侧的递归处理(注意如果是用cdq处理dp的话，因为右区间要先被左区间更新，再去更新同侧的，所以顺序有点不一样，还要sort)；

那么右上角在上半部分，左下角在下半部分，于是我们枚举右上角；

对于右上角的限制，只要有一个横坐标小于他且纵坐标小于他的点就不行，于是我们先按照按照横坐标排序，这样我们就能用递增单调栈求出左边第一个小于他的位置；

然后我们考虑左下角的限制，如果存在一个横坐标大于他且纵坐标大于他的就不合法，我们依然是按照横坐标排序，然后我们相当于是要维护一个递减的单调栈，每个点后面的点就都合法；

我们需要在下方找到横坐标比右上角限制点大的第一个点，然后从这个点开始往后到栈顶就都是合法的了，这个我们在单调栈中二分，然后直接统计即可；

具体实现方法就是两个单调栈一起建，然后更新；

//MADE BY QT666

#include<iostream>

#include<cstdio>

#include<algorithm>

#include<cstring>

using namespace std;

typedef long long ll;

const int N=300050;

struct data{

    int x,y;

}a[N],b[N];

int n,tp1,tp2,p[N],q[N];

ll ans;

bool cmp(const data &a,const data &b){

    return a.y<b.y;

};

int find(int x,int l,int r){

    int ret=r+1;

    while(l<=r){

	int mid=(l+r)>>1;

	if(a[q[mid]].x>=x) r=mid-1,ret=mid;

	else l=mid+1;

    }

    return ret;

}

void solve(int l,int r){

    if(l==r) return;

    int mid=(l+r)>>1;

    solve(l,mid);solve(mid+1,r);

    tp1=0;tp2=0;int j=l;

    for(int 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,1,tp2)+1;

    }

    j=l;int k=mid+1;

    for(int i=l;i<=r;i++){

	if(j<=mid&&(a[j].x<a[k].x||k>r)) b[i]=a[j++];

	else b[i]=a[k++];

    }

    for(int 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+1+n,cmp);

    solve(1,n);printf("%lld\n",ans);

    return 0;

}