BZOJ1007【HNOI2008】水平可见直线

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1007

【分析】

先按直线的斜率排好序。斜率相同的直线只有截距最大的那个可见。

然后记录一个栈，以及一个数组D[i]记录位于栈中的第i条直线与栈中第i-1条直线交点的横坐标。（以斜率从小到大排序为例，那么栈中第i条直线可见部分仅在x>D[i]处）

如果当前加入的直线与栈顶直线的交点横坐标小于D[top]（即当前加入直线覆盖了栈顶直线的可见部分），则栈顶直线不可见，如此依次弹出栈顶直线，再加入当前直线。

【代码】

#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn=50010;
const double eps=1e-10;
 
struct Line
{
    double a,b;
    int ID;
    bool operator < (const Line &p) const {return a<p.a || (abs(a-p.a)<eps && b>p.b);}
}L[maxn];
typedef Line Node;
double GetX(Line x,Line y) {return (y.b-x.b)/(x.a-y.a);}
 
int Q[maxn],top;
int I[maxn];
double D[maxn];
int n;
 
int main()
{
    scanf("%d",&n);
    for (int i=0;i<n;i++) scanf("%lf%lf",&L[i].a,&L[i].b),L[i].ID=i+1;
    sort(L,L+n);
    Q[top]=0;D[top]=-1e80;
    for (int i=1;i<n;i++)
    {
        if (abs(L[i].a-L[Q[top]].a)<eps) continue;
        D[++top]=GetX(L[i],L[Q[top-1]]);
        while (D[top]-D[top-1]<eps) D[--top]=GetX(L[i],L[Q[top-1]]);
        Q[top]=i;
    }
    for (int i=0;i<=top;i++) I[L[Q[i]].ID]=1;
    for (int i=1;i<=n;i++) if (I[i]) printf("%d ",i);
    return 0;
}