P4557 [JSOI2018]战争 凸包的闵可夫斯基和

https://www.luogu.com.cn/problem/P4557

  • 题目描述很简单,给定两个凸包,我们设为 a , b a,b a,b,现在问给 b b b凸包一个整体偏移,偏移量为 k → \overrightarrow{k} k ,给定 k → = ( x , y ) \overrightarrow{k}=(x,y) k =(x,y),问这两个凸包有没有公共部分
  • 换句话说,就是问给定的 k → \overrightarrow{k} k 是否满足 b + k → = a b+\overrightarrow{k}=a b+k =a,即 k → = a − b \overrightarrow{k}=a-b k =ab

首先介绍一下闵可夫斯基和
定义:两个欧几里得空间点集的和,也称作两个空间的膨胀集,点集 A A A B B B的闵可夫斯基和定义为 A + B = { a + b , a ∈ A , b ∈ B } A+B=\{a+b,a\in A,b\in B\} A+B={a+b,aA,bB}

  • 求取方法是把这两个空间内部所有的向量按照极角排序,然后找到一个起点,按照极角从小到大依次放边,最后构成的形状就是两个空间的闵可夫斯基和

那么这道题我们其实要求的是 a a a − b -b b的闵可夫斯基和,只要 k → \overrightarrow{k} k 终点在两个空间的膨胀集内部,或者边界,就说明会发生冲突,否则不会

  • 什么叫做凸集?凸集的几何意义是如果一个集合任意两个元素连线上的点也在集合内部,那么这个集合就是一个凸集,如下图,左侧是一个凸集但是右侧却不是一个凸集;代数定义是对于 S S S的子集内的所有 x x x y y y,并且在区间 [ 0 , 1 ] [0,1] [0,1]内所有的值 t t t,如果点 ( 1 − t ) x + t y (1-t)x+ty (1t)x+ty也属于 S S S,那么 S S S就是一个凸集
    在这里插入图片描述
  • 那么这里产生了一个疑问,为什么两个凸集的闵可夫斯基和仍然是一个凸集?我们可以进行如下证明,设 A , B A,B A,B是两个凸集,设 C = A + B , α 1 , α 2 ∈ A , β 1 , β 2 ∈ B C=A+B,\alpha_1,\alpha_2\in A,\beta_1,\beta_2\in B C=A+B,α1,α2A,β1,β2B
  • 由于 ( 1 − t ) α 1 + t α 2 ∈ A (1-t)\alpha_1+t\alpha_2\in A (1t)α1+tα2A ( 1 − t ) β 1 + t β 2 ∈ B (1-t)\beta_1+t\beta_2\in B (1t)β1+tβ2B,所以根据凸集和的定义有 ( 1 − t ) α 1 + t α 2 + ( 1 − t ) β 1 + t β 2 ∈ C (1-t)\alpha_1+t\alpha_2+(1-t)\beta_1+t\beta_2\in C (1t)α1+tα2+(1t)β1+tβ2C
  • 所以 ( 1 − t ) ( α 1 + β 1 ) + t ( α 2 + β 2 ) ∈ C (1-t)(\alpha_1+\beta_1)+t(\alpha_2+\beta_2)\in C (1t)(α1+β1)+t(α2+β2)C又因为根据凸集和的定义有 α 1 + β 1 ∈ C , α 2 + β 2 ∈ C \alpha_1+\beta_1\in C,\alpha_2+\beta_2\in C α1+β1C,α2+β2C,这就证明了C是一个凸集

接下来回到这道题,梳理一下思路,我们首先要求出 A A A − B -B B两个凸集的闵可夫斯基和,这个结果是一个凸包,求闵可夫斯基和我们可以使用固定一个初始点然后进行归并,始终把外面的线加进去,这里使用向量叉积来判断线是否在外部, − B -B B的意思就是坐标都取相反数,算出来要再跑一次 A n d r e w Andrew Andrew防止出现共线,然后进行点是否在凸包内部的二分判定,程序如下

#include <bits/stdc++.h>

using namespace std;
#define db double
const db eps = 1e-10;
const int MAXN = 1e5 + 100;
int sgn(db x){
    if(fabs(x) < eps) return 0;
    return x < 0 ? -1 : 1;
}
struct Point{
    db x, y;
    Point(){}
    Point(db x, db y): x(x), y(y){}
    Point operator + (const Point &B)const{
        return Point(x + B.x, y + B.y);
    }
    Point operator - (const Point &B)const{
        return Point(x - B.x, y - B.y);
    }
    bool operator < (const Point &B)const{
        return sgn(x - B.x) < 0 || (sgn(x - B.x) == 0 && sgn(y - B.y) < 0);
    }
    bool operator == (const Point &B)const{
        return sgn(x - B.x) == 0 && sgn(y - B.y) == 0;
    }
}s1[MAXN], s2[MAXN], ch1[MAXN], ch2[MAXN];
Point M[MAXN], ans[MAXN];
typedef Point Vector;
db Cross(Vector A, Vector B){
    return A.x * B.y - A.y * B.x;
}
int Convex_hull(Point *s, Point *ch, int n){
    int v = 0;
    sort(s, s + n);
    n = unique(s, s + n) - s;
    for(int i=0;i<n;i++){
        while(v > 1 && sgn(Cross(ch[v - 1] - ch[v - 2], s[i] - ch[v - 2])) <= 0){
            v -= 1;
        }
        ch[v++] = s[i];
    }
    int j = v;
    for(int i=n-2;i>=0;i--){
        while(v > j && sgn(Cross(ch[v - 1] - ch[v - 2], s[i] - ch[v - 2])) <= 0){
            v -= 1;
        }
        ch[v++] = s[i];
    }
    if(n > 1) v -= 1;
    return v;
}
int check(Point A, Point *ch, int n){
    int l = 0;
    int r = n - 1;
    while(r - l > 1){
        int mid = ((r - l) >> 1) + l;
        db a1 = Cross(ch[mid] - ch[0], A - ch[0]);
        db a2 = Cross(ch[mid + 1] - ch[0], A - ch[0]);
        if(sgn(a1) >= 0 && sgn(a2) <= 0){
            if(sgn(Cross(ch[mid + 1] - ch[mid], A - ch[mid])) >= 0) return 1;
            return 0;
        }else if(sgn(a1) < 0){
            r = mid;
        }else{
            l = mid;
        }
    }
    return 0;
}
int Minkowski_sum(int n, int m){
    int tot;
    for(int i=1;i<n;i++){
        s1[i] = ch1[i] - ch1[i - 1];
    }
    s1[n] = ch1[0] - ch1[n - 1];
    for(int i=1;i<m;i++){
        s2[i] = ch2[i] - ch2[i - 1];
    }
    s2[m] = ch2[0] - ch2[m - 1];
    ans[tot = 0] = s1[0] + s2[0];
    int pt1, pt2;
    pt1 = pt2 = 1;
    while(pt1 <= n && pt2 <= m){
        tot += 1;
        ans[tot] = ans[tot - 1] + (sgn(Cross(s1[pt1], s2[pt2])) >= 0 ? s1[pt1++] : s2[pt2++]);
    }
    while(pt1 <= n){
        tot += 1;
        ans[tot] = ans[tot - 1] + s1[pt1++];
    }
    while(pt2 <= m){
        tot += 1;
        ans[tot] = ans[tot - 1] + s2[pt2++];
    }
    return tot;
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int n, m, q;
    db x, y;
    cin >> n >> m >> q;
    for(int i=0;i<n;i++){
        cin >> s1[i].x >> s1[i].y;
    }
    int a = Convex_hull(s1, ch1, n);
    for(int i=0;i<m;i++){
        cin >> s2[i].x >> s2[i].y;
        s2[i].x = -s2[i].x;
        s2[i].y = -s2[i].y;
    }
    int b = Convex_hull(s2, ch2, m);
    int c = Minkowski_sum(a, b);
    c = Convex_hull(ans, M, c);
    while(q--){
        cin >> x >> y;
        cout << check(Point(x, y), M, c) << '\n';
    }
    return 0;
}

热门文章

暂无图片
编程学习 ·

exe4j详细使用教程(附下载安装链接)

一、exe4j介绍 ​ exe4j是一个帮助你集成Java应用程序到Windows操作环境的java可执行文件生成工具&#xff0c;无论这些应用是用于服务器&#xff0c;还是图形用户界面&#xff08;GUI&#xff09;或命令行的应用程序。如果你想在任务管理器中及Windows XP分组的用户友好任务栏…
暂无图片
编程学习 ·

AUTOSAR从入门到精通100讲(126)-浅谈车载充电系统通信方案

01 引言 本文深入研究车载充电系统策略,设计出一套基于电动汽车电池管理系统与车载充电机的CAN通信协议,可供电动汽车设计人员参考借鉴。 02 电动汽车充电系统通讯网络 电动汽车整车控制系统中采用的是CAN总线通信方式,由一个整车内部高速CAN网络、内部低速CAN网络和一个充电…
暂无图片
编程学习 ·

CMake(九):生成器表达式

当运行CMake时&#xff0c;开发人员倾向于认为它是一个简单的步骤&#xff0c;需要读取项目的CMakeLists.txt文件&#xff0c;并生成相关的特定于生成器的项目文件集(例如Visual Studio解决方案和项目文件&#xff0c;Xcode项目&#xff0c;Unix Makefiles或Ninja输入文件)。然…
暂无图片
编程学习 ·

47.第十章 网络协议和管理配置 -- 网络配置(八)

4.3.3 route 命令 路由表管理命令 路由表主要构成: Destination: 目标网络ID,表示可以到达的目标网络ID,0.0.0.0/0 表示所有未知网络,又称为默认路由,优先级最低Genmask:目标网络对应的netmaskIface: 到达对应网络,应该从当前主机哪个网卡发送出来Gateway: 到达非直连的网络,…
暂无图片
编程学习 ·

元宇宙技术基础

请看图&#xff1a; 1、通过AR、VR等交互技术提升游戏的沉浸感 回顾游戏的发展历程&#xff0c;沉浸感的提升一直是技术突破的主要方向。从《愤怒的小鸟》到CSGO,游戏建模方式从2D到3D的提升使游戏中的物体呈现立体感。玩家在游戏中可以只有切换视角&#xff0c;进而提升沉浸…
暂无图片
编程学习 ·

flink的伪分布式搭建

一 flink的伪分布式搭建 1.1 执行架构图 1.Flink程序需要提交给 Job Client2.Job Client将作业提交给 Job Manager3.Job Manager负责协调资源分配和作业执行。 资源分配完成后&#xff0c;任务将提交给相应的 Task Manage。4.Task Manager启动一个线程以开始执行。Task Manage…
暂无图片
编程学习 ·

十进制正整数与二进制字符串的转换(C++)

Function one&#xff1a; //十进制数字转成二进制字符串 string Binary(int x) {string s "";while(x){if(x % 2 0) s 0 s;else s 1 s;x / 2;}return s; } Function two&#xff1a; //二进制字符串变为十进制数字 int Decimal(string s) {int num 0, …
暂无图片
编程学习 ·

[含lw+源码等]微信小程序校园辩论管理平台+后台管理系统[包运行成功]Java毕业设计计算机毕设

项目功能简介: 《微信小程序校园辩论管理平台后台管理系统》该项目含有源码、论文等资料、配套开发软件、软件安装教程、项目发布教程等 本系统包含微信小程序做的辩论管理前台和Java做的后台管理系统&#xff1a; 微信小程序——辩论管理前台涉及技术&#xff1a;WXML 和 WXS…
暂无图片
编程学习 ·

树莓派驱动DHT11温湿度传感器

1&#xff0c;直接使用python库 代码如下 import RPi.GPIO as GPIO import dht11 import time import datetimeGPIO.setwarnings(True) GPIO.setmode(GPIO.BCM)instance dht11.DHT11(pin14)try:while True:result instance.read()if result.is_valid():print(ok)print(&quo…
暂无图片
编程学习 ·

ELK简介

ELK简介 ELK是三个开源软件的缩写&#xff0c;Elasticsearch、Logstash、Kibana。它们都是开源软件。不过现在还新增了一个 Beats&#xff0c;它是一个轻量级的日志收集处理工具(Agent)&#xff0c;Beats 占用资源少&#xff0c;适合于在各个服务器上搜集日志后传输给 Logstas…
暂无图片
编程学习 ·

Linux 基础

通常大数据框架都部署在 Linux 服务器上&#xff0c;所以需要具备一定的 Linux 知识。Linux 书籍当中比较著名的是 《鸟哥私房菜》系列&#xff0c;这个系列很全面也很经典。但如果你希望能够快速地入门&#xff0c;这里推荐《Linux 就该这么学》&#xff0c;其网站上有免费的电…
暂无图片
编程学习 ·

Windows2022 无线网卡装不上驱动

想来 Windows2022 和 windows10/11 的驱动应该差不多通用的&#xff0c;但是死活装不上呢&#xff1f; 搜一下&#xff0c;有人提到 “默认安装时‘无线LAN服务’是关闭的&#xff0c;如果需要开启&#xff0c;只需要在“添加角色和功能”中&#xff0c;选择开启“无线LAN服务…
暂无图片
编程学习 ·

【嵌入式面试宝典】版本控制工具Git常用命令总结

目录 创建仓库 查看信息 版本回退 版本检出 远程库 Git 创建仓库 git initgit add <file> 可反复多次使用&#xff0c;添加多个文件git commit -m <message> 查看信息 git status 仓库当前的状态git diff 差异对比git log 历史记录&#xff0c;提交日志--pret…
暂无图片
编程学习 ·

用Postman生成测试报告

newman newman是一款基于nodejs开发的可以运行postman脚本的工具&#xff0c;使用Newman&#xff0c;可以直接从命令运行和测试postman集合。 安装nodejs 下载地址&#xff1a;https://nodejs.org/en/download/ 选择自己系统相对应的版本内容进行下载&#xff0c;然后傻瓜式安…
暂无图片
编程学习 ·

Java面向对象之多态、向上转型和向下转型

文章目录前言一、多态二、引用类型之间的转换Ⅰ.向上转型Ⅱ.向下转型总结前言 今天继续Java面向对象的学习&#xff0c;学习面向对象的第三大特征&#xff1a;多态&#xff0c;了解多态的意义&#xff0c;以及两种引用类型之间的转换&#xff1a;向上转型、向下转型。  希望能…