题目描述
JOI村有一片荒地,上面竖着N个稻草人,村民们每年多次在稻草人们的周围举行祭典。
有一次,JOI村的村长听到了稻草人们的启示,计划在荒地中开垦一片田地。和启示中的一样,田地需要满足以下条件:
田地的形状是边平行于坐标轴的长方形;
左下角和右上角各有一个稻草人;
田地的内部(不包括边界)没有稻草人。
给出每个稻草人的坐标,请你求出有多少遵从启示的田地的个数
输入
第一行一个正整数N,代表稻草人的个数
接下来N行,第i行(1<=i<=N)包含2个由空格分隔的整数Xi和Yi,表示第i个稻草人的坐标
输出
输出一行一个正整数,代表遵从启示的田地的个数
样例输入
4
0 0
2 2
3 4
4 3
样例输出
3
题解
自己yy出来的分治+单调栈+二分
如果直接对于每个点求以它为顶点的满足条件的矩形数目比较难求,所以考虑分治处理。
按x分治,处理完左右区间后处理左边对右边的影响。
此时左半部分的x严格小于右半部分的x,如果再按x排序则没有意义,所以按照y排序。
按照y从小到大排序后,遍历区间内所有的点。对于每个右半部分的点,在左半部分里寻找答案。
我们思考:如果在左半部分扫到了这样两个点a、b:ay>by且ax>bx,那么在按y从小到大遍历时,剩下的点的y都比这两个点大,画图可知b不可能再形成矩形,所以弹掉。
即对左半部分维护一个从栈底到栈顶x值递减的单调栈,碰到不满足条件的则弹出。
再考虑右半部分:如果扫到了这样两个点a、b:ay>by且ax<bx,那么这两个点是互不影响的,直接寻找答案即可。所以b没有意义,所以弹掉。而当ay>by且ax>bx时,b限制了a形成的矩形的范围,所以应当保留,并且在遍历到a时在左半部分中二分。
即对右半部分维护一个从栈底到栈顶x值递增的单调栈,碰到不满足条件的则弹出。
综上,我们维护两个单调栈,每次遍历到一个点,就把它压到对应的单调栈中,如果这个点是右半部分的点,就在左半部分的单调栈中二分求出比栈顶元素y值大的点的个数,并累加到答案中。注意左右部分的单调栈是不同的,因为它们的意义是不同的。
这样做的时间复杂度是$O(nlog^2n)$,亲测使用归并排序,二分的常数极小,可以使时间减到5s左右。
#include <cstdio> #include <algorithm> #define N 200010 using namespace std; struct data { int x , y; }a[N] , tmp[N]; int s1[N] , t1 , s2[N] , t2; long long ans; bool cmp(data a , data b) { return a.x < b.x; } int getnum(int t) { int l = 1 , r = t1 , mid , tmp = t1 + 1; while(l <= r) { mid = (l + r) >> 1; if(a[s1[mid]].y >= t) tmp = mid , r = mid - 1; else l = mid + 1; } return t1 - tmp + 1; } void solve(int l , int r) { if(l >= r) return; int mid = (l + r) >> 1 , tx = a[mid].x , i , p1 = l , p2 = mid + 1; solve(l , mid); solve(mid + 1 , r); for(i = l ; i <= r ; i ++ ) { if(p2 > r || (p1 <= mid && a[p1].y < a[p2].y)) tmp[i] = a[p1 ++ ]; else tmp[i] = a[p2 ++ ]; } t1 = t2 = 0; for(i = l ; i <= r ; i ++ ) { a[i] = tmp[i]; if(a[i].x <= tx) { while(t1 && a[i].x > a[s1[t1]].x) t1 -- ; s1[++t1] = i; } else { while(t2 && a[i].x < a[s2[t2]].x) t2 -- ; s2[++t2] = i; ans += getnum(a[s2[t2 - 1]].y); } } } int main() { int n , i; scanf("%d" , &n); for(i = 1 ; i <= n ; i ++ ) scanf("%d%d" , &a[i].x , &a[i].y); sort(a + 1 , a + n + 1 , cmp); solve(1 , n); printf("%lld " , ans); return 0; }