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
;
int n;
int S
,top,ans
,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;
}