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
0 0
2 2
3 4
4 3
Sample Output
3
HINT
所有满足要求的田地由下图所示:
![](http://www.lydsy.com/JudgeOnline/upload/201508/aa.jpg)
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
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~CDQ分治+二分~
我们纵坐标在区间[l,r]内的点按照y排序分成上下两块,然后对于每个在下面一块中的点寻找在上面一块中的能与之匹配的点。
所以我们每次维护单调栈来记录目前的点,然后在其中二分查找符合要求的点更新答案。
ans要开long long!
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long
int n,p[200005],q[200005];
ll ans;
struct node{
int x,y;
bool operator < (const node&u) const {return y<u.y;}
}a[200005],b[200005];
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
while(ch>='0' && ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
int find(int x,int l,int r)
{
int mid;
while(l+1<r)
{
mid=l+r>>1;
if(a[q[mid]].x<x) l=mid;
else r=mid;
}
return l;
}
void sol(int l,int r)
{
if(l==r) return;
int mid=l+r>>1,top1=0,top2=0,l1=l,l2=mid+1;
sol(l,mid);sol(mid+1,r);
for(int i=mid+1,j=l;i<=r;i++)
{
while(top1 && a[p[top1]].y>a[i].y) top1--;
p[++top1]=i;
for(;j<=mid && a[j].x<a[i].x;j++)
{
while(top2 && a[q[top2]].y<a[j].y) top2--;
q[++top2]=j;
}
ans+=top2-find(a[p[top1-1]].x,0,top2+1);
}
for(int i=l;i<=r;i++) b[i]=(l1<=mid && a[l1].x<a[l2].x || l2>r) ? a[l1++]:a[l2++];
for(int i=l;i<=r;i++) a[i]=b[i];
}
int main()
{
n=read();
for(int i=1;i<=n;i++) a[i].x=read(),a[i].y=read();
sort(a+1,a+n+1);sol(1,n);
printf("%lld\n",ans);
return 0;
}