模板 - 凸包

模板 - 凸包

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N = 300005;

// 存储二维平面点
struct ip {
    int x, y, i;
    ip(int x = 0, int y = 0) : x(x), y(y), i(0) { }
    void ri(int _i) {
        scanf("%d%d", &x, &y);
        i = _i;
    }
};

// iv表示一个向量类型，其存储方式和ip一样
typedef ip iv;

// 先比较x轴再比较y轴，
bool operator < (const ip &a, const ip &b) {
    return a.x == b.x ? a.y < b.y : a.x < b.x;
}

// 两点相减得到的向量
iv operator - (const ip &a, const ip &b) {
    return iv(a.x - b.x, a.y - b.y);
}

// 计算a和b的叉积（外积）
//可以判断方向性，例如向量a为(2,0)，向量b为(2,2) 返回正,即向量a在向量b逆时针偏转45度处
ll operator ^ (const iv &a, const iv &b) {
    return (ll)a.x * b.y - (ll)a.y * b.x;
}

/* 请在这里定义你需要的全局变量 */

// 计算二维点数组a的凸包，将凸包放入b数组中，下标均从0开始
// a, b：如上
// n：表示a中元素个数
// 返回凸包元素个数
int convex(ip *a, ip *b, int n) {
    /* 请在这里设计你的算法 */
    sort(a,a+n);
    //如有重点，必须去重
    //n = unique(a,a+n)-a;

    int m = 0;
    for (int i = 0;i < n;i ++){ //这里是向上凸的部分
        for (;m > 1 && ((b[m-1] - b[m-2])^(a[i] - b[m-2])) > 0;m --);
        b[m++] = a[i];
    }

    for (int i = n-2,t = m;i >= 0;i --){ //这里是向下凸的部分
        for (;m-t > 0 && ((b[m-1] - b[m-2])^(a[i] - b[m-2])) > 0;m --);
        b[m++] = a[i];
    }
    return m-1;
}


ip a[N], b[N];

int main() {
    int n;
    scanf("%d", &n);
    for (int i = 0; i < n; ++i)
        a[i].ri(i + 1);
    int m = convex(a, b, n);
    printf("%d\n", m); //凸包的点数
    return 0;
}