BZOJ1007 [HNOI2008]水平可见直线 贪心+栈

给出n条直线，求从y值无穷大处能看见的直线编号.n<=500000.

按照斜率k递增为第一关键字，与y轴截距b递减为第二关键字对直线排序.

如果两条直线斜率不同，那新直线要被看到，必须与之前直线的交点在当前交点的左侧。

如果斜率相同，则b值大的可见。

没有卡精度还是可以的

#include<bits/stdc++.h>
#define LL long long
#define clr(x,i) memset(x,i,sizeof(x))
using namespace std;
const int N=500005;
const double eps=1e-8;
struct line{
    double k,b;int id;
    friend bool operator <(line a,line b){
        return fabs(a.k-b.k)<eps ? a.b>b.b : a.k<b.k;
    }
}st[N];
int n;
int S[N],top,ans[N],cnt;
inline double cross(int x,int y)
{
    return (st[y].b-st[x].b)/(st[x].k-st[y].k);
}
void solve()
{
    S[1]=1;top=1;
    for(int i=2;i<=n;i++)
    {
        if(fabs(st[i].k-st[i-1].k)<eps)continue;
        while(top>1&&cross(S[top],S[top-1])+eps>cross(S[top],i))
          top--;
        S[++top]=i;
    }
    for(int i=1;i<=top;i++)
      ans[++cnt]=st[S[i]].id;
    sort(ans+1,ans+cnt+1);
    for(int i=1;i<=cnt;i++)
      printf("%d ",ans[i]);
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
      scanf("%lf%lf",&st[i].k,&st[i].b),st[i].id=i;
    sort(st+1,st+n+1);
    solve();
    return 0;
}