曼哈顿距离几何意义 矩形面积并

这篇博客探讨了如何利用曼哈顿距离的几何意义解决矩形面积并问题。通过二分答案,构建以坐标点为中心,对角线长度为2*mid的正方形,并讨论了将正方形顺时针旋转45度以简化求解过程。文章提到了在旋转后如何维护正方形的边界,并指出在实际计算中,边长乘以sqrt(2)并不影响求交的结果。虽然数据集比较简单,但未深入讨论交集中是否存在整点的情况。

11.7

思路:
二分ans,把所有长度>mid的(lf,rg)都搞出来,
然后我们要找出一个点对(u,v),使得|x - u| + |y - v| <= mid。
考虑几何意义,曼哈顿距离的图像限定。
把一个区间(lf,rg)转换成一个点(x,y)。
对于每个点(x,y),我们构建出一个以(x,y)为几何中心,
对角线长2*mid,且对角线平行于坐标轴的正方形。
满足条件的(u,v)就在这个正方形当中(可在边上)。
那么要让所有的(x,y)满足条件,就是看所有正方形是否有交。
我们发现斜45度的正方形交不好求,考虑把正方形扳正。
也就是对于所有的点,都围绕原点顺时针旋转45度。
(当然我们维护一条对角线两端的两个点就好了)
这里维护扳正后正方形的左下角右上角,也就是原正方形的上下两个顶点。
我们发现旋转之后的坐标也不好维护,但是如果我们把正方形边长*sqrt(2),
(x,y)就会转到(y+x,y-x),很好维护。
而且所有的边都*sqrt(2)对求交并没有影响。
其实还应该讨论一下交中有没有整点,不过数据水还是过了。

#include <cstdio>
#include <iostream> 
#include <algorithm>
#define N 2000010
using namespace std;

int n, m;

struct Edge{
    int lf, rg;
}ed[N];

struct Matrix{
    int x1, x2, y1, y2;//1为左下角 2为右上角 
};

bool cmp(Edge a, Edge b){ return a.rg - a.lf > b.rg - b.lf;}

Matrix tra(int i, int K){//把一段区间转成一个矩阵 
    int x = ed[i].lf, y = ed[i].rg;
    Matrix cc;
    cc.x1 = x + y - K, cc.y1 = y - K - x;
    cc.x2 = x + y + K, cc.y2 = y + K - x;
    return cc;
}

bool unionn(Matrix &A, Matrix B){//合并维护交集 
    if(A.x2 < B.x1 || A.x1 > B.x2 || A.y1 > B.y2 || A.y2 < B.y1) return 0;//
    A.x1 = max(A.x1, B.x1), A.y1 = max(A.y1, B.y1);
    A.x2 = min(A.x2, B.x2), A.y2 = min(A.y2, B.y2);
    return 1;
}

bool check(int x){
    int flag = 1, cnt = 0;
    Matrix cc;
    for(int i=1; i<=m; i++)
        if(ed[i].rg - ed[i].lf > x){
            if( !cnt ) cc = tra(i, x), cnt = 1;
            else flag = unionn(cc, tra(i, x));
            if(!flag) return 0;
        }
        else break;
    return 1;
}

int main(){
    freopen ("b.in", "r", stdin);
    freopen ("b.out", "w", stdout);
    int lf = 0, rg = 0, ans = 0;
    scanf("%d%d", &n, &m);
    for(int i=1; i<=m; i++){
        scanf("%d%d", &ed[i].lf, &ed[i].rg);
        rg = max(rg, ed[i].rg - ed[i].lf);
    }
    sort(ed+1, ed+m+1, cmp);
    while(lf <= rg){
        int mid = (lf + rg) >> 1;
        if( check( mid ) )
            rg = mid - 1, ans = mid;
        else lf = mid + 1;
    }
    printf("%d\n", ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值