BZOJ4237(cdq分治+单调栈）

最新推荐文章于 2022-10-05 14:28:42 发布

creatorx

最新推荐文章于 2022-10-05 14:28:42 发布

阅读量320

点赞数 1

CC 4.0 BY-SA版权

分类专栏： cdq分治&整体二分

本文链接：https://blog.youkuaiyun.com/creatorx/article/details/76855163

cdq分治&整体二分专栏收录该内容

2 篇文章

订阅专栏

本文介绍了一道关于寻找特定矩形区域的问题——稻草人问题，并使用CDQ分治算法进行解答。该算法通过将点集按坐标排序并利用单调栈来高效计算满足条件的矩形数量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

4237: 稻草人

Time Limit: 40 Sec Memory Limit: 256 MB
Submit: 873 Solved: 382
[ Submit][ Status][ Discuss]

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

解题思路：用cdq分治，先把所有点按Y坐标从小到大排序，然后对于两个分治的区间求解答案时，先把两个区间各自按x坐标从小到大排序，然后左区间用一个单调栈维护一个y坐标递减的序列，右区间用一个单调栈维护一个y坐标单调递增的序列，然后在单调递增的单调栈中找一个次大的点（作为为标准，因为这个点对当前枚举的右上端点威胁最大，只要左区间中的点的横坐标大于这个点的横坐标，这个点就不在当前枚举点与左区间的点组成的长方形中间了），在单调递减的单调栈中用二分找出所以横坐标大于这个标准点的横坐标的点就是答案了。

觉得这个用单调栈有点难想，还是太菜了。。。。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 200000 + 10;
int N;
long long sum;
int Sta1[maxn];
int Sta2[maxn];//两个单调栈
int top1;
int top2;
struct Point{
    int x, y;
}P[maxn];
bool cmpx(Point p1, Point p2)
{
    return p1.x < p2.x;
}
bool cmpy(Point p1, Point p2)
{
    return p1.y < p2.y;
}
void cdq(int L, int R)
{
    if(L == R) return;
    int mid = (L + R)>>1;
    cdq(L, mid);
    cdq(mid + 1, R);
    sort(P + L, P + mid + 1, cmpx);
    sort(P + mid + 1, P + R + 1, cmpx);
    int ll = L;
    top1 = top2 = 0;
    for(int i = mid + 1; i <= R; i++)
    {
        while(top1 > 0 && P[i].y < P[Sta1[top1]].y)
        {
            top1--;
        }
        Sta1[++top1] = i;
        while(ll <= mid && P[ll].x < P[i].x)
        {
            while(top2 > 0 && P[ll].y > P[Sta2[top2]].y) top2--;
            Sta2[++top2] = ll++;
        }
        int l = 1, r = top2;
        int re = -1;
        int judge = P[Sta1[top1 - 1]].x;
        if(top1 <= 1) judge = -1;
        while(l <= r)
        {
            int m = (l + r)>>1;
            if(P[Sta2[m]].x > judge)
            {
                re = m;
                r = m - 1;
            }
            else l = m + 1;
        }
        if(re > 0) sum += top2 - re + 1;
    }

}
void init()
{
    sum = 0;
    top1 = top2 = 0;//栈顶指针
}
int main()
{
    scanf("%d", &N);
    init();
    for(int i = 1; i <= N; i++)
    {
        scanf("%d%d", &P[i].x, &P[i].y);
    }
    sort(P + 1, P + N + 1, cmpy);
    cdq(1, N);
    printf("%lld\n", sum);
    return 0;
}