1007: [HNOI2008]水平可见直线

Description

　　在xoy直角坐标平面上有n条直线L1,L2,...Ln,若在y值为正无穷大处往下看,能见到Li的某个子线段,则称Li为
可见的,否则Li为被覆盖的.
例如,对于直线:
L1:y=x; L2:y=-x; L3:y=0
则L1和L2是可见的,L3是被覆盖的.
给出n条直线,表示成y=Ax+B的形式(|A|,|B|<=500000),且n条直线两两不重合.求出所有可见的直线.

Input

　　第一行为N(0 < N < 50000),接下来的N行输入Ai,Bi

Output

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

Sample Input

3
-1 0
1 0
0 0

Sample Output

       1    2

从Y轴无穷远出看 应该看到的是一个凹形图，为了求出凹形图上有哪些直线的子线段，优先对斜率k升序排序，k相同对b降序排序。维护一个单调栈，如果当前直线能使栈顶

直线完全被覆盖则弹出，最后进栈。栈中元素就是答案。

#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;
const int maxm = 500005;
struct node
{
	double k, b;
	int id;
	bool operator<(const node &r)const
	{
		if (r.k != k)
			return k < r.k;
		return b > r.b;
	}
};
node p[maxm];
int s[maxm], len = 1, flag[maxm] = { 0 };
double judge(node x,node y);
int main()
{
	int n, i, j, k, sum;
	scanf("%d", &n);
	for (i = 1;i <= n;i++)
	{
		scanf("%lf%lf", &p[i].k, &p[i].b);
		p[i].id = i;
	}
	sort(p + 1, p + 1 + n);
	s[1] = 1; flag[p[1].id] = 1;
	for (i = 2;i <= n;i++)
	{
		if (p[i].k - p[i - 1].k < 1e-8)	continue;
		while (len > 1 && judge(p[i], p[s[len]]) <= judge(p[s[len]], p[s[len - 1]]))len--;
		s[++len] = i;
	}
	for (i = 1;i <= len;i++)
		flag[p[s[i]].id] = 1;
	for (i = 1;i <= n;i++)
	{
		if (flag[i])
			printf("%d ", i);
	}
	printf("\n");
	return 0;
}
double judge(node x, node y)
{
	return (y.b - x.b) / (x.k - y.k);
}