二维凸包详解

凸包大意:给一堆点,让你求一个最小的多边形使得其囊括所有点。实际上就是拿条绳子在一堆钉子外面绕一圈,绕成的图形就是所求多边形,即绕绳法。但这里并不讨论绕绳法,而是介绍一种叫Graham的算法。前置芝士:叉积(a×b=a.x·b.y+a.y·b.x)叉积不同于点积,叉积的结果是向量,而点积的结果是标量。设平面直角坐标系有三点A(x0,y0),B(x1,y1),C(x2,y2)。有共起点...
摘要由CSDN通过智能技术生成

凸包大意:给一堆点,让你求一个最小的多边形使得其囊括所有点。

实际上就是拿条绳子在一堆钉子外面绕一圈,绕成的图形就是所求多边形,即绕绳法。
但这里并不讨论绕绳法,而是介绍一种叫Graham的算法。


前置芝士:叉积(a×b=a.x·b.y+a.y·b.x)

叉积不同于点积,叉积的结果是向量,而点积的结果是标量。
设平面直角坐标系有三点A(x0,y0),B(x1,y1),C(x2,y2)。
有共起点p0的两个向量AB(x1-x0,y1-y0),AC(x2-x0,y2-y0),它们的叉积表示为AB×AC=(x1-x0)•(y2-y0)+(y1-y0)•(x2-x0)。

此时有三种情况:
  • 若p0p1×p0p2>0,那么向量p0p1在向量p0p2的逆时针方向;
  • 若p0p1×p0p2<0,那么向量p0p1在向量p0p2的顺时针方向;
  • 若p0p1×p0p2=0,那么向量p0p1与向量p0p2共线。

如下图,a为p0p1,b为p0p2。

至于为什么请读者自行百度思考(其实我也不知道)。

另一个用途:

若有向量AB,AC,则AB×AC=2S∆ABC。
这样就可以把求出来的多边形分割成一个个三角形,再统计这些三角形的总面积,就是多边形的面积。

在这里插入图片描述

算法:

Graham算法的实现步骤如下:

  • 设stk数组为答案栈,存的是答案点。使用时相邻两点连接就是答案多边形的一条边了。
  • 先取一个最左下的点x(下优先),易知x为所求多边形上一个点,可以直接入栈。以x对其余n-1个点做极角排序,这样可以保证是逆时针确定所求多边形的边。
    极角排序:按夹角的大小排序。
  • 接下来枚举剩下的n-1个点:设起点为stk[n-1],终点为stk[n] 的向量为a;起点为stk[n-1],终点为当前点的向量为b。用叉积检查a是否在b的逆时针方向(多边形的边尽量靠外即尽量靠顺时针才能包住所有点),如果a在b的逆时针方向(或共线)就将stk[n]出栈,重复执行该操作,直到不满足条件退出。最后加入当前点。(不管怎样都要加,大不了下次删就是了)
    tip:共线的情况视题目条件而定。
  • 简单地说,Graham就是一个检测上个答案点是否有用,即这个点连stk[n-1]是否比stk[n]连stk[n-1]更优,是一个类似于贪心的过程。
  • 枚举完后已经得到了答案数组stk,可以为所欲为了
大概长这个样:

(原作者不晓得,哪里都有这玩意,作者可以联系我删)
在这里插入图片描述





模板:

#include<bits/stdc++.h>
using namespace std;
int n,tot=1,xx,yy;
struct data{
   
    double x,y;
};
data nd[60060];
data stk[60060];
double ans=0,mid;
double cross
  • 1
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值