利用凸包的定义,可以确定最远点对一定包含于凸包点集中,所有先利用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.