Andrew’s Monotone Chain 图示
先以 X 座标排序,当 X 座标相同则以 Y 座标排序。从起点开始,按照顺序扫描,找到下半凸包。再从终点开始,按照相反顺序扫描,找到上半凸包。合起来就是完整的凸包。
解决了凸包有重叠的点、共线的点、退化成线段和点的问题(在Graham’s Scan可能遇到的问题)。
- 时间复杂度:O(NlogN),其中排序的时间 O(NlogN) ,加上包围的时间 O(N) 。
具体思想:
- 排序,根据某个坐标轴为主进行排序
- 从x最小的点开始
- 利用类似Graham’s Scan算法,利用栈去寻找下半包
- 从x最大的点开始
- 利用类似Graham’s Scan算法,利用栈去寻找上半包
- 两个半包结合,即是整个凸包
(上半包和下半包可以合并成一个while函数)
Andrew’s Monotone Chain 代码
// P为平面上撒的点。设定为10点。
// CH为凸包上的顶点。设定为逆时针方向排列。可以视作一个stack。
struct Point {int x, y;} P[10], CH[10+1];
// 向量OA叉积向量OB。大于零表示从OA到OB为逆时针旋转。
double cross(Point o, Point a, Point b)
{
return (a.x - o.x) * (b.y - o.y) - (a.y - o.y) * (b.x - o.x);
}
bool compare(Point a, Point b)
{
return (a.x < b.x) || (a.x == b.x && a.y < b.y);
}
void Andrew_monotone_chain()
{
//排序
sort(P, P+10, compare);
int m = 0; // m 是凸包的顶点数
// 包下半部
for (int i=0; i<10; ++i)
{
while (m >= 2 && cross(CH[m-2], CH[m-1], P[i]) <= 0) m--;
CH[m++] = P[i];
}
// 包上半部,不用再包入方才包过的终点,但会再包一次起点
for (int i=10-2, t=m+1; i>=0; --i)
{
while (m >= t && cross(CH[m-2], CH[m-1], P[i]) <= 0) m--;
CH[m++] = P[i];
}
m--; // 最后一个点是重复出现两次的起点,所以减1。
}