凸包大意:给一堆点,让你求一个最小的多边形使得其囊括所有点。
实际上就是拿条绳子在一堆钉子外面绕一圈,绕成的图形就是所求多边形,即绕绳法。
但这里并不讨论绕绳法,而是介绍一种叫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