BZOJ1007: [HNOI2008]水平可见直线

BZOJ1007: [HNOI2008]水平可见直线

https://lydsy.com/JudgeOnline/problem.php?id=1007

分析:

• 这玩意是个半平面交，但是他可以用凸包来对偶搞一搞。
• 把直线按斜率排序，然后搞一搞就做完了

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <cmath>
using namespace std;
typedef double f2;
#define N 50050
int n,ans[N],S[N],tp;
struct A {
	f2 k,b;
	int id;
	bool operator < (const A &u) const {
		return abs(k-u.k)<1e-8 ? b<u.b : k<u.k;
	}
}a[N];
f2 X(A &x,A &y) {
	return -(x.b-y.b)/(x.k-y.k);
}
int main() {
	scanf("%d",&n);
	int i;
	for(i=1;i<=n;i++) {
		scanf("%lf%lf",&a[i].k,&a[i].b); a[i].id=i;
	}
	sort(a+1,a+n+1);

	for(i=1;i<=n;i++) {
		while(tp>1&&((abs(a[S[tp]].k-a[i].k)<1e-8)||X(a[S[tp-1]],a[S[tp]])>=X(a[S[tp]],a[i]))) tp--;
		if(tp==1&&(abs(a[S[tp]].k-a[i].k)<1e-8)) tp--;
		S[++tp]=i;
	}
	for(i=1;i<=tp;i++) ans[++ans[0]]=a[S[i]].id;
	sort(ans+1,ans+ans[0]+1);
	for(i=1;i<=ans[0];i++) printf("%d ",ans[i]);
}