算法导论学习笔记（五）计算几何之最远点对

利用凸包的定义，可以确定最远点对一定包含于凸包点集中，所有先利用Graham’s scan算法求出凸包点集，然后利用旋转卡壳，求出最远距离。
旋转卡壳的变种：以凸包中某一条边为起点，先利用叉积求出离这一条边最远的点，而当边顺次选择时，最远点也会顺次转移，使得时间复杂度降为O(n)。

伪代码

ROTATE_CALIPERS( A )
    A[size] = A[0]
    let q = 1
    let dis = -
    for i=0 to A.size-1
        while GET_CROSS( A[i+1]-A[i], A[q]-A[i]) > GET_CROSS( A[i+1]-A[i], A[q+1]-A[i]) 
            q = (q+1) % A.size
        dis = MAX( dis, GET_DIS(A[i],A[q]), GET_DIS(A[i+1], A[q+1]) )
    return dis

GET_FARTHEST_DISTANCE( A )
    let B is a convex hull of A
    return ROTATE_CALIPERS( B )

代码实现

#include <cmath>
#include <algorithm>
const double ACCURACY = 0.0000001;
typedef struct Point
{
    int x,y;
}VectorLine;
VectorLine operator- ( Point a, Point b) 
{
    VectorLine c;
    c.