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)互不相同。

Source

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

一年前就写过这题。。结果各种细节不会处理就半途而废了。今天在bzoj上看到就想着赶紧把这题过了

我们考虑分治。对所有点的x坐标先排序。

这样当处理区间[l,r]的时候，

所有点就分成了[l,mid]和[mid+1,r]两块。左右的点x坐标互不影响。

然后我们再把两块分别按照y坐标从大到小排序。

维护一个单调栈，左边保证x单调减，右边保证x单调增

每加入左边的一个点k，就把右边y坐标大于该点的点加入

对于点k，可以取到的右边的点的y坐标在k.y到k-1.y之间

这里二分一下，再相减，就可以求出每个点可以形成多少个矩形

不断分治，即可得到ans

#include<cstdio>
#include<algorithm>
using namespace std;
struct point
{
	int x,y;
}a[200001];
inline bool cmp1(point x,point y)
{
	if(x.x<y.x)
		return true;
	return false;
}
inline bool cmp2(point x,point y)
{
	if(x.y>y.y)
		return true;
	return false;
}
point q1[200001],q2[200001];
inline int find(int x,int lim)
{
	int l=1,r=lim;
	int mid;
	while(l<=r)
	{
		mid=(l+r)/2;
		if(q1[mid].y<=x)
			r=mid-1;
		else
			l=mid+1;
	}
	return r;
}
long long ans;
inline void solve(int l,int r)
{
	if(l>=r)
		return ;
	int mid=(l+r)/2;
	solve(l,mid);
	sort(a+l,a+mid+1,cmp2);
	sort(a+mid+1,a+1+r,cmp2);
	int i;
	int top1=0,top2=0,p=mid+1;
	int ll,rr;
	q1[0].y=q2[0].y=2100000000;
	for(i=l;i<=mid;i++)
	{
		while(p<=r&&a[p].y>a[i].y)
		{
			while(top1!=0&&a[p].x<q1[top1].x)
				top1--;
			top1++;
			q1[top1]=a[p];
			p++;
		}
		while(top2!=0&&a[i].x>q2[top2].x)
			top2--;
		top2++;
		q2[top2]=a[i];
		ll=find(q2[top2-1].y,top1)+1;
		rr=find(q2[top2].y,top1);
		if(ll<=rr)
			ans+=(long long)(rr-ll+1);
	}
	sort(a+l,a+mid+1,cmp1);
	sort(a+mid+1,a+1+r,cmp1);
	solve(mid+1,r);
}
int main()
{
	//freopen("03-05.txt","r",stdin);
	//freopen("03-05.out","w",stdout);
	int n;
	scanf("%d",&n);
	int i;
	for(i=1;i<=n;i++)
		scanf("%d%d",&a[i].x,&a[i].y);
	sort(a+1,a+1+n,cmp1);
	solve(1,n);
	printf("%lld\n",ans);
	return 0;
}