该算法是在gift-warping算法的基础上通过排序来预处理数据然后降低时间复杂度的
#include <iostream>
#include <memory>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
#define MAXN 100001
using namespace std;
/* 向量 */
struct Vertex {
double x, y;
Vertex(const double & _x = double(), const double & _y = double()) : x(_x), y(_y) {}
};
/* 点 */
struct Point {
double x, y;
Point(const double & _x = double(), const double & _y = double()) : x(_x), y(_y) {}
} points[MAXN];
/* 计算叉积 */
inline double cross(const Vertex & a, const Vertex & b) {
return a.x * b.y - a.y * b.x;
}
/* 计算叉积<a, b> * <a, c> */
inline double cross(const Point & a, const Point & b, const Point & c) {
Vertex vab = Vertex(b.x - a.x, b.y - a.y);
Vertex vac = Vertex(c.x - a.x, c.y - a.y);
return cross(vab, vac);
}
/* 计算两点间的距离 */
inline d