重叠矩形面积不会算？来这篇文章学扫描线算法就对了

最新推荐文章于 2024-09-13 10:21:40 发布

干货没人看

最新推荐文章于 2024-09-13 10:21:40 发布

阅读量1k

点赞数 13

文章标签： c语言 c++ 算法 python 经验分享

本文链接：https://blog.csdn.net/qq_36916968/article/details/141949990

版权

你有没有遇到过这样的问题：给你一大堆的矩形，让你找出在哪个地方重叠得最厉害？听起来有点像数学老师出的高难度题目对吧？很多人都会选择老办法，挨个比对，逐个判断，看起来很有道理，但实际上效率低到爆！今天，我要给你介绍一个高效又“扫”得干净的算法——扫描线法（Sweep Line Algorithm），用它解决矩形最大重叠问题，那叫一个清爽！

想象一下，一个简单的例子

我们先来点简单的，比如你有四个矩形，它们相互重叠在一起，想要知道在哪个区域重叠最多，怎么做？很多人的第一反应是逐个矩形拿出来比较，挨个算出重叠面积，最后找出最大值。但你有没有想过，当你面对几十上百个矩形时，这种方法就要让你哭了。因为这种暴力计算的时间复杂度可是(O(n^2)) 级别的！

别怕，有扫荡全场的“扫线法”

聪明人都知道，要解决这种问题，首先要找个高效的思路！你有没有注意到一个现象：矩形的重叠情况，只在“关键时刻”发生变化，比如在某个矩形的边进入或退出某个区域的时候。既然重叠变化只在这些关键时刻发生，那我们就只需要在这些时刻进行计算就好了。这就是扫描线法的核心思想！

扫描线法怎么玩？一条线扫出所有答案！

让我们一步一步来看这个神奇的算法到底怎么运作。

收集事件点：开门与关门

想象你站在一个原野上，有一大堆矩形在你的左边，你面前有一根横向的线，这根线会从左往右慢慢移动（我们称之为“扫描线”）。每当扫描线碰到一个矩形的左边（我们叫它“开门”事件），这个矩形就会被“激活”；而当扫描线碰到一个矩形的右边（我们叫它“关门”事件），这个矩形就会“关闭”。
按顺序处理事件：从左到右依次扫描

接下来，我们将所有的“开门”和“关门”事件按 x 坐标排序。扫描线从最左边开始，逐一处理这些事件。当遇到“开门”事件时，我们把这个矩形加入到当前的活动矩形集；遇到“关门”事件时，我们把相应的矩形移出活动集。
重叠高度怎么计算？一套数据结构搞定！

你可能会问：“在每个位置，重叠面积怎么快速计算？”这就要借助数据结构的力量了！聪明的程序员都会用一棵“有序集合”（比如 STL 中的 set）或者“平衡树”来存储和维护活动矩形集的纵坐标区间信息。通过动态插入和删除这些区间，我们就可以快速计算当前的重叠面积。
实时更新最大重叠面积

每次扫描线遇到一个事件点，我们就计算一下当前的重叠面积，并更新我们记录的最大重叠面积。就这样，一路扫过去，所有的重叠信息都能轻松搞定！

看代码，算法真理藏在细节中

让我们直接看看这个算法在代码中的应用，保证让你秒懂！

#include <stdio.h>
#include <stdlib.h>

typedef struct {
    int x;      // 事件的x坐标
    int y1, y2; // 矩形的y区间
    int type;   // 事件类型：1表示矩形左边缘（进入），-1表示右边缘（离开）
} Event;

// 按x坐标排序事件
int compareEvents(const void* a, const void* b) {
    return ((Event*)a)->x - ((Event*)b)->x;
}

#define MAX_EVENTS 1000
#define MAX_Y 10000

int sweepLine(Event events[], int n) {
    // 用于记录重叠的y区间
    int active[MAX_Y] = {0};
    int maxOverlap = 0;

    qsort(events, n, sizeof(Event), compareEvents); // 按x排序

    for (int i = 0; i < n; i++) {
        int y1 = events[i].y1;
        int y2 = events[i].y2;
        int type = events[i].type;

        // 更新活动区间
        for (int y = y1; y < y2; y++) {
            active[y] += type;
        }

        // 计算当前重叠
        int currentOverlap = 0;
        for (int y = 0; y < MAX_Y; y++) {
            if (active[y] > 0) {
                currentOverlap++;
            }
        }

        if (currentOverlap > maxOverlap) {
            maxOverlap = currentOverlap;
        }
    }

    return maxOverlap;
}

int main() {
    Event events[MAX_EVENTS] = {
        {1, 1, 5, 1},  {4, 1, 5, -1},  // 第一个矩形 (1,1) 到 (4,5)
        {2, 2, 6, 1},  {5, 2, 6, -1},  // 第二个矩形 (2,2) 到 (5,6)
        {3, 3, 7, 1},  {6, 3, 7, -1},  // 第三个矩形 (3,3) 到 (6,7)
    };
    int n = 6; // 6个事件点

    printf("最大重叠矩形数：%d\n", sweepLine(events, n));

    return 0;
}