CF535E 计算几何

题目大意：

$n$个运动员，有两段赛道，第$i$个运动员在第一段赛道跑的速度是$v_{i1}$，第二段赛道跑的速度是$v_{i2}$，现在不知道每种赛道的长度，询问有哪些运动员可能最后赢得比赛

解题思路：

• 可以先将运动员速度从小到大排序，第一关键字为第一种速度

• 通过排序，可以先将一些必不可能成为赢家的筛选掉（其两种速度都小于某一个人的两种速度），并且对于两种速度都相同的去重处理

• 对于三个人，中间那个人可能成为赢家的条件是：（已经排序去重）（$s$为第一种速度，$t$为第二种速度，x为第一段长度，y为第二段长度）
  $$\begin{gathered} \frac{x}{s_2}-\frac{x}{s_1}+\frac{y}{t_2}-\frac{y}{t_1}\le 0,a_1=\frac{1}{s_2}-\frac{1}{s_1}<0,b_1=\frac{1}{t_2}-\frac{1}{t_1}>0 \\ \frac{x}{s_2}-\frac{x}{s_3}+\frac{y}{t_2}-\frac{y}{t_3}\le 0,a_2=\frac{1}{s_2}-\frac{1}{s_3}>0,b_1=\frac{1}{t_2}-\frac{1}{t_3}<0 \\ \Downarrow \\ {b}_{1}y\le -a_{1}x\Rightarrow y\le -\frac{a_1}{b_1}x \\ {b}_{2}y\le -a_{2}x\Rightarrow y\ge -\frac{a_2}{b_2}x \\ 成为赢家的条件就是要都满足上述条件：则-\frac{a_1}{b_1}\ge -\frac{a_2}{b_2}\Rightarrow \frac{a_1}{b_1}\le \frac{a_2}{b_2} \\ 经过化简就变成了：s_3{t}_1(s_1-s_2)(t_3-t_2)\ge s_1{t}_3(s_3-s_2)(t_1-t_2) \end{gathered}$$

• 然后根据这个式子维护一下栈里元素即可

更好的想法：

• 对于式子$\frac{A}{a_i}+\frac{B}{b_i}=T$，题目就转为对于点$(\frac{1}{a_i},\frac{1}{b_i})$，如果存在$(A,B)$使得在该点对最小，则输出
• 原式可以看成：$Ax+By=T\to y=-\frac{A}{B}x+\frac{T}{B}$
• 由于$(x,y)$给定，而斜率恒为负数，所以只要维护一个左下凸包即可，且最终维护凸包的式子就上面的的式子

AC代码：

#include <bits/stdc++.h>
#define ft first
#define sd second
#define pb push_back
#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0) //不能跟puts混用
#define seteps(N) fixed << setprecision(N)
#define endl "\n"
const int maxn = 2e5 + 10;
using namespace std;
typedef long long ll;
typedef double db;
typedef pair<int, int> pii;
const ll mod = 1e9 + 7;
int n, m;
struct Speed {
    int v1, v2, id;
    bool operator < (Speed p) {
        if (v1 != p.v1) return v1 < p.v1;
        else if (v2 != p.v2) return v2 < p.v2;
        else return id < p.id;
    }
} sp[maxn], ps[maxn];
int st[maxn], s;
bool vis[maxn];
db slope(int s1, int s2, int t1, int t2) {
    return (1.0 / s2 - 1.0 / s1) / (1.0 / t2 - 1.0 / t1);
}
int main() {
    cin >> n;
    for (int i = 1; i <= n; i++) cin >> sp[i].v1 >> sp[i].v2, sp[i].id = i;
    sort (sp + 1, sp + n + 1);
    m = 1;
    ps[1] = sp[1];
    for (int i = 2; i <= n; i++) {
        while (m >= 1 && ps[m].v2 <= sp[i].v2) m--;
        ps[++m] = sp[i];
    }
    for (int i = 1; i <= m; i++) {
        while (s >= 2 && 1ll * ps[i].v1 * ps[st[s - 1]].v2 * (ps[st[s - 1]].v1 - ps[st[s]].v1) * (ps[i].v2 - ps[st[s]].v2) < 1ll * ps[st[s - 1]].v1 * ps[i].v2 * (ps[i].v1 - ps[st[s]].v1) * (ps[st[s - 1]].v2 - ps[st[s]].v2))
            s--;
        st[++s] = i;
    }
    for (int i = 1; i <= s; i++) vis[ps[st[i]].id] = true;
    for (int i = n; i > 1; i--) {
        if (sp[i].v1 == sp[i - 1].v1 && sp[i].v2 == sp[i - 1].v2)
            vis[sp[i - 1].id] = vis[sp[i].id];
    }
    for (int i = 1; i <= n; i++)
        if (vis[i])
            cout << i << " ";
    cout << endl;
    return 0;
}