BZOJ1007/洛谷P3194 [HNOI2008]水平可见直线(凸包(单调栈)求n条直线可见部分)

题目

给出 n 条直线，表示成 y=Ax+B 的形式(∣A∣,∣B∣≤500000)，且 n 条直线两两不重合。

求出在y轴的无穷大处往下看时，所有可见的直线。

思路来源

夏老师

https://www.luogu.com.cn/blog/LiM-817/solution-p3194

题解

要不是夏老师说eduF是这个原题，我都不知道直线能求凸包……

其实是求了个下凸壳，

考虑按斜率增序，斜率相同时按截距增序排序，加入的时候，

如果栈顶小于两条直线，直接加入

如果和栈顶直线平行，说明完全把上一条压在底下，

如果最后一条直线和栈顶直线交点x1，栈顶直线和次栈顶直线交点x2，x1<x2的话栈顶直线根本不会露出来，弹栈

单调栈维护一下，发现这个形状就是下凸壳

代码

#include<bits/stdc++.h>
using namespace std;
typedef double db;
const db eps=1e-8;
const int N=5e4+10;
struct line{
    db k,b;
    int id;
}l[N],stk[N];
bool ok[N];
int n,c;
bool cmp(line a,line b){
    return a.k<b.k||(fabs(a.k-b.k)<eps && a.b<b.b);
}
db crossx(line a,line b){
    return (b.b-a.b)/(a.k-b.k);
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;++i){
        scanf("%lf%lf",&l[i].k,&l[i].b);
        l[i].id=i;
    }
    sort(l+1,l+n+1,cmp);
    for(int i=1;i<=n;++i){
        while(c){
            //这里可以统一写作b1*k2<=b2*k1 注意分母符号 但要求分开考虑c>1和c=1
            if(fabs(stk[c].k-l[i].k)<eps)c--;
            else if(c>1 && crossx(stk[c],l[i])<=crossx(stk[c-1],stk[c]))c--;
            else break;
        }
        stk[++c]=l[i];
    }
    for(int i=1;i<=c;++i){
        ok[stk[i].id]=1;
    }
    for(int i=1;i<=n;++i){
        if(ok[i]){
            printf("%d ",i);
        }
    }
    return 0;
}