1. 程式人生 > >【HNOI 2008】水平可見直線

【HNOI 2008】水平可見直線

【題目】

傳送門

題目描述:

x o y xoy 直角座標平面上有 n n 條直線 l

1 , l 2 , . . . , l
n
l_1,l_2,...,l_n ,若在y值為正無窮大處往下看,能見到Li的某個子線段,則稱Li為可見的,否則Li為被覆蓋的。

例如,對於直線 l 1 l_1

y = x y=x l 2 l_2 y = x y=-x l 3 l_3 y = 0 y=0 ,則 l 1 l_1 l 2 l_2 是可見的, l 3 l_3 是被覆蓋的。給出 n n 條直線,表示成 y = A x + B y=Ax+B 的形式( A , B |A|,|B| 500000 500000 ),且 n n 條直線兩兩不重合。求出所有可見的直線.

輸入格式:

第一行為 n n ( 0 < n < 50000 ) (0 < n < 50000) ,接下來的 n n 行輸入 A i , B i A_i,B_i

輸出格式:

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

樣例資料:

輸入
3
-1 0
1 0
0 0

輸出
1 2

【分析】

其實看這道題的題意,應該還是很容易可以想到半平面交吧

還是按照斜率大小排序,每次加新線的時候就比較一下交點位置,然後判斷是否需要彈掉棧頂

不過這道題比半平面交模板不同的是不用彈掉隊首和隊尾(因為它們肯定能被看見)

注意這道題有平行線的情況,我們只用保留 B 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;
}