1. 程式人生 > >[HNOI2008]水平可見直線

[HNOI2008]水平可見直線

個數字 printf 空格 描述 const include 行為 ... 半平面

題目描述

在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條直線兩兩不重合.求出所有可見的直線.

輸入輸出格式

輸入格式:

第一行為N(0 < N < 50000),接下來的N行輸入Ai,Bi

輸出格式:

從小到大輸出可見直線的編號,兩兩中間用空格隔開,最後一個數字後面也必須有個空格

輸入輸出樣例

輸入樣例#1:
3
-1 0
1 0
0 0
輸出樣例#1:
1 2


半平面交板題。
#include<bits/stdc++.h>
#define ll long long
#define maxn 50005
using namespace std;
const double eps=1e-9;
inline int zt(double x){
	if(fabs(x)<eps) return 0;
	return x>0?1:-1;
}
struct lines{
	int num;
	double k,b;
	bool operator <(const lines &U)const{
		return zt(k-U.k)?zt(k-U.k)<0:b>U.b;
	}
}l[maxn];
int n,s[maxn],tp;
bool v[maxn];

inline double X(lines x,lines y){
	return (x.b-y.b)/(y.k-x.k);
}

inline void getplain(){
	for(int i=1;i<=n;i++){
		if(tp&&!zt(l[i].k-l[s[tp]].k)) continue;
		while(tp>1&&zt(X(l[i],l[s[tp]])-X(l[i],l[s[tp-1]]))<=0) tp--;
		s[++tp]=i;
	}
} 

int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%lf%lf",&l[i].k,&l[i].b),l[i].num=i;
	sort(l+1,l+n+1);
	
	getplain();
	
	for(int i=1;i<=tp;i++) v[l[s[i]].num]=1;
	for(int i=1;i<=n;i++) if(v[i]) printf("%d ",i);
	
	return 0;
}

  



[HNOI2008]水平可見直線