洛谷P3194 [HNOI2008]水平可见直线题解

题目链接：P3194 [HNOI2008]水平可见直线

题意：在$ x-y$ 直角坐标平面上有 $n$ 条直线 $L_1,L_2,…L_n$ ，若在 $y$ 值为正无穷大处往下看，能见到 $L_i$ 的某个子线段，则称 $L_i$ 为可见的，否则 $L_i$ 为被覆盖的。例如，对于直线: $L_1:y=x$ ; $L_2:y=-x$ ; $L_3:y=0$ ; 则 $L_1$ 和 $L_2$ 是可见的， $L_3$ 是被覆盖的。给出 $n$ 条直线，表示成 $y = A x + B$ 的形式( $\le 500000$ )，且 $n$ 条直线两两不重合，求出所有可见的直线。

一条直线C被另外两条直线A和B所覆盖

它的充分必要条件为直线C在A,B交点处的值更小一些

反之，直线C也可能对答案有贡献

那么交点的坐标是什么呢？
$A_1x+B_1=A_2x+B_2\\ x=\dfrac{B_2-B_1}{A_1-A_2}$
则
$A_1\dfrac{B_2-B_1}{A_1-A_2} + B_1 \ge A_3\dfrac{B_2-B_1}{A_1-A_2} + B_3\\ \color{red}\dfrac{B_2-B_1}{A_2-A_1}\ge\dfrac{B_3-B_1}{A_3-A_1}$
可以发现这个式子很像斜率

如果记 $P_i(A_i,B_i)$

那么答案就是点集 ${P_i\}$ 构成的上凸壳

为什么是上凸壳？显然我们选定的是斜率大且尽可能截距大

代码如下

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define INF 0x3f3f3f3f3f3f3f3f
#define N (int)(1e5+15)
struct vct
{
    double x,y;int id;
    vct operator-(const vct &o)const
    {
        return (vct){x-o.x,y-o.y};
    }
}p[N];
int n,stk[N],top;
double cross(vct a,vct b)
{
    return a.x*b.y-a.y*b.x;
}
signed main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);cout.tie(0);
    cin >> n;
    for(int i=1; i<=n; i++)
    {
        cin >> p[i].x >> p[i].y;
        p[i].id=i;
    }
    sort(p+1,p+1+n,[](vct a,vct b)
    {
        return a.x==b.x?a.y>b.y:a.x>b.x;
    });
    stk[++top]=1;
    for(int i=2; i<=n; i++)
    {
        if(p[i-1].x==p[i].x)continue;
        while(top>1&&cross(p[stk[top]]-p[stk[top-1]],p[i]-p[stk[top]])<=0)
            --top;
        stk[++top]=i;
    }
    sort(stk+1,stk+1+top,[](int a,int b)
    {
        return p[a].id<p[b].id;
    });
    for(int i=1; i<=top; i++)
        cout << p[stk[i]].id << " ";
    return 0;
}