2D凸包算法（三）：Andrew's Monotone Chain

Andrew’s Monotone Chain 图示

先以 X 座标排序，当 X 座标相同则以 Y 座标排序。从起点开始，按照顺序扫描，找到下半凸包。再从终点开始，按照相反顺序扫描，找到上半凸包。合起来就是完整的凸包。
解决了凸包有重叠的点、共线的点、退化成线段和点的问题（在Graham’s Scan可能遇到的问题）。

• 时间复杂度：O(NlogN)，其中排序的时间 O(NlogN) ，加上包围的时间 O(N) 。

具体思想：

1. 排序，根据某个坐标轴为主进行排序
2. 从x最小的点开始
3. 利用类似Graham’s Scan算法，利用栈去寻找下半包
4. 从x最大的点开始
5. 利用类似Graham’s Scan算法，利用栈去寻找上半包
6. 两个半包结合，即是整个凸包
   （上半包和下半包可以合并成一个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。
}