凸包模板

凸包

描述
给定n个二维平面上的点，求他们的凸包。
输入
第一行包含一个正整数n。
接下来n行，每行包含两个整数x,y，表示一个点的坐标。
输出
令所有在凸包极边上的点依次为p1,p2,…,pm（序号），其中m表示点的个数，请输出以下整数：
(p1 × p2 × … × pm × m) mod (n + 1)
样例1输入

10
7 9
-8 -1
-3 -1
1 4
-3 9
6 -4
7 5
6 6
-6 10
0 8

样例1输出

7

样例1解释

所以答案为(9 × 2 × 6 × 7 × 1 × 5) % (10 + 1) = 7

限制
3 ≤ n ≤ 10^5
所有点的坐标均为范围(-10⁵, 10⁵)内的整数，且没有重合点。每个点在(-10⁵, 10⁵) × (-10⁵, 10⁵)范围内均匀随机选取
极边上的所有点均被视作极点，故在输出时亦不得遗漏
时间：4 sec
空间：512 MB

注意这里代码这里是逆时针顺序的凸包，需要先求出下凸壳，也就是<0然后再求出上凸壳。

//先求下凸壳      for(; m>1 && ((b[m-1]-b[m-2])^(a[i]-b[m-2])) < 0; --m);

完整代码如下所示

#include <vector>
#include <stdio.h>
#include <iostream>
#include <cstdio>
#include <algorithm>
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的叉积（外积）
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 && ((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), ans = m;
    for (int i = 0; i < m; ++i)
        ans = ((ll)ans * b[i].i) % (n + 1);
    printf("%d\n", ans);
    return 0;
}