【HNOI 2008】水平可见直线

【题目】

传送门

题目描述：

在 $xoy$ 直角坐标平面上有 $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=Ax+B$ 的形式（$|A|,|B|$ ≤ $500000$），且 $n$ 条直线两两不重合。求出所有可见的直线.

输入格式：

第一行为 $n$ $(0<n<50000)$，接下来的 $n$ 行输入 $A_i,B_i$

输出格式：

从小到大输出可见直线的编号，两两中间用空格隔开，最后一个数字后面也必须有个空格

样例数据：

输入
3
-1 0
1 0
0 0

输出
1 2

【分析】

其实看这道题的题意，应该还是很容易可以想到半平面交吧

还是按照斜率大小排序，每次加新线的时候就比较一下交点位置，然后判断是否需要弹掉栈顶

不过这道题比半平面交模板不同的是不用弹掉队首和队尾（因为它们肯定能被看见）

注意这道题有平行线的情况，我们只用保留 $B$ 值最大的那条

【代码】

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 50005
#define eps 1e-8
using namespace std;
struct line
{
	int id;
	double k,b;
}l[N],sta[N];
bool comp1(const line &p,const line &q)
{
	if(fabs(p.k-q.k)<eps)
	  return p.b<q.b;
	return p.k<q.k;
}
bool comp2(const line &p,const line &q)
{
	return p.id<q.id;
}
double interx(line a,line b)
{
	return (b.b-a.b)/(a.k-b.k);
}
int main()
{
	int n,i;
	scanf("%d",&n);
	for(i=1;i<=n;++i)
	{
		l[i].id=i;
		scanf("%lf%lf",&l[i].k,&l[i].b);
	}
	sort(l+1,l+n+1,comp1);
	int num=0,top=0;
	sta[++top]=l[1];
	for(i=2;i<=n;++i)
	{
		if(fabs(l[i].k-l[top].k)<eps)  top--;
		while(top>1&&interx(l[i],sta[top])<=interx(sta[top],sta[top-1]))  top--;
		sta[++top]=l[i];
	}
	sort(sta+1,sta+top+1,comp2);
	for(i=1;i<=top;++i)
	  printf("%d ",sta[i].id);
	return 0;
}