凸包基础

凸包入门

用途

从一堆结点中，获取构成面积最大的图形所需的顶点。

结点定义

将结点定义为一个变量，以及可以进行相减运算。

代码实现

// p[] 输入的点
// s[] 外围的点（构成凸包的点）
struct point {
    double x, y;

    //对减法符号进行重载
    point friend operator-(point a, point b) {
        return {a.x - b.x, a.y - b.y};
    }
} p[maxn], s[maxn];

叉乘

用途

判断结点相对位置

示例

$\overrightarrow{p1}$×$\overrightarrow{p2}$（假设两向量均以原点出为起点）

1. >0：$\overrightarrow{p1}$在$\overrightarrow{p2}$的顺时针方向
2. <0：$\overrightarrow{p1}$在$\overrightarrow{p2}$的逆时针方向
3. =0：共线

判断一个点是否在一条线段的左边还是右边
例如：相交的两条直线AB、CD

1. C、D在AB两侧，则 ($\overrightarrow{AB}$×$\overrightarrow{AC}$)* ($\overrightarrow{AB}$×$\overrightarrow{AD}$) <=0
2. A、B在CD两侧，则 ($\overrightarrow{CD}$×$\overrightarrow{CA}$)* ($\overrightarrow{CD}$×$\overrightarrow{CB}$) <=0
3. <0 表示在两边，=0表示在直线上

代码实现

// 叉积
double X(point a, point b) {
    return a.x * b.y - a.y * b.x;
}

// 携带极点的叉积
double multi(point p1, point p2, point p3) {
    return X(p1 - p3, p2 - p3);
}

Graham

操作

将输入的点从 1~n 进行排序，以 p[1]（所有点中最左下角的点，如若仅有非负数，则可以直接以原点 p[0] 作为极点）为极点，先按极角由小到大排序，极角如若相同，则按距离由小到大排序排序。

代码实现

// 距离公式
double dis(point a, point b) {
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}

// 极角排序
bool cmp(point a, point b) {
    double use = multi(a, b, p[1]);
    if (use > 0) {
        return true;
    } else if (use == 0) {
        return dis(a, p[1]) < dis(b, p[1]);
    } else {
        return false;
    }
}

// 获取左下的点作为极点
void makePool() {
    int k = 1;
    for (int i = 2; i <= n; i++) {
        if (p[i].y < p[k].y || (p[k].y == p[i].y && p[i].x < p[k].x)) {
            k = i;
        }
    }
    swap(p[1], p[k]);
}

void graham() {
    makePool();
    sort(p + 2, p + n + 1, cmp);
    s[1] = p[1];
    s[2] = p[2];
    top = 2;
    for (int i = 3; i <= n; i++) {
        while (top >= 2 && multi(s[top - 1], s[top], p[i]) <= 0) {
            top--;
        }
        s[++top] = p[i];
    }
//    for (int i = 1; i <= top; i++) {
//        cout << s[i].x << " " << s[i].y << endl;
//    }
}

例题

仅仅是求凸包的周长。

护城河

题目描述

为了防止口渴的食蚁兽进入他的农场，Farmer John决定在他的农场周围挖一条护城河。农场里一共有N(8<=N<=5,000)股泉水，并且，护城河总是笔直地连接在河道上的相邻的两股泉水。护城河必须能保护所有的泉水，也就是说，能包围所有的泉水。泉水一定在护城河的内部，或者恰好在河道上。当然，护城河构成一个封闭的环。
挖护城河是一项昂贵的工程，于是，节约的FJ希望护城河的总长度尽量小。请你写个程序计算一下，在满足需求的条件下，护城河的总长最小是多少。
所有泉水的坐标都在范围为(1…10,000,000,1…10,000,000)的整点上，一股泉水对应着一个唯一确定的坐标。并且，任意三股泉水都不在一条直线上。
以下是一幅包含20股泉水的地图，泉水用"*"表示：

图中的直线，为护城河的最优挖掘方案，即能围住所有泉水的最短路线。
路线从左上角起，经过泉水的坐标依次是：(18,0),(6,-6),(0,-5),(-3,-3),(-17,0),(-7,7),(0,4),(3,3)。绕行一周的路径总长为70.8700576850888(…)。答案只需要保留两位小数，于是输出是70.87。

输入描述：

第1行: 一个整数，N
第2…N+1行: 每行包含2个用空格隔开的整数，x[i]和y[i]，即第i股泉水的位置坐标。

输出描述：

第1行: 输出一个数字，表示满足条件的护城河的最短长度。保留两位小数。

输入

20
2 10
3 7
22 15
12 11
20 3
28 9
1 12
9 3
14 14
25 6
8 1
25 1
28 4
24 12
4 15
13 5
26 5
21 11
24 4
1 8

输出

70.87

完整代码

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>

using namespace std;

const int maxn = 1e5 + 10;
int top, n;
double ans;

struct point {
    double x, y;

    //对减法符号进行重载
    point friend operator-(point a, point b) {
        return {a.x - b.x, a.y - b.y};
    }
} p[maxn], s[maxn];

// 距离公式
double dis(point a, point b) {
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}

// 叉积
double X(point a, point b) {
    return a.x * b.y - a.y * b.x;
}

// 携带极点的叉积
double multi(point p1, point p2, point p3) {
    return X(p1 - p3, p2 - p3);
}

// 极角排序
bool cmp(point a, point b) {
    double use = multi(a, b, p[1]);
    if (use > 0) {
        return true;
    } else if (use == 0) {
        return dis(a, p[1]) < dis(b, p[1]);
    } else {
        return false;
    }
}

// 获取左下的点作为极点
void makePool() {
    int k = 1;
    for (int i = 2; i <= n; i++) {
        if (p[i].y < p[k].y || (p[k].y == p[i].y && p[i].x < p[k].x)) {
            k = i;
        }
    }
    swap(p[1], p[k]);
}

void graham() {
    makePool();
    sort(p + 2, p + n + 1, cmp);
    s[1] = p[1];
    s[2] = p[2];
    top = 2;
    for (int i = 3; i <= n; i++) {
        while (top >= 2 && multi(s[top - 1], s[top], p[i]) <= 0) {
            top--;
        }
        s[++top] = p[i];
    }
}

int main() {
    cin >> n;
    for (int i = 1; i <= n; i++) {
        cin >> p[i].x >> p[i].y;
    }
    graham();
    for (int i = 1; i < top; i++) {
        ans += dis(s[i], s[i + 1]);
    }
    ans += dis(s[1], s[top]);
    printf("%.2lf\n", ans);
    return 0;
}