【BZOJ1007】[HNOI2008]水平可见直线【半平面交】

【题目链接】

学习一发计算几何基础。

/* Telekinetic Forest Guard */
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;

typedef double DB;

const int maxn = 50005;
const DB eps = 1e-6;

struct _line {
	DB k, b;
	int id;
} l[maxn];

int n, sta[maxn], top;

inline bool cmp(_line x, _line y) {
	return fabs(x.k - y.k) < eps ? x.b < y.b : x.k < y.k;
}

inline bool cmp2(int x, int y) {
	return l[x].id < l[y].id;
}

inline DB cosp(_line x, _line y) {
	return (x.b - y.b) / (y.k - x.k);
}

inline void insert(int id) {
	while(top) {
		if(fabs(l[sta[top]].k - l[id].k) < eps) top--;
		else if(top > 1 && cosp(l[id], l[sta[top - 1]]) <= cosp(l[sta[top]], l[sta[top - 1]])) top--;
		else break;
	}
	sta[++top] = id;
}

int main() {
	scanf("%d", &n);
	for(int i = 1; i <= n; i++) scanf("%lf%lf", &l[i].k, &l[i].b), l[i].id = i;

	sort(l + 1, l + 1 + n, cmp);

	top = 0;
	for(int i = 1; i <= n; i++) insert(i);
	sort(sta + 1, sta + 1 + top, cmp2);
	for(int i = 1; i <= top; i++) printf("%d ", l[sta[i]].id);
	return 0;
}