【HNOI 2008】水平可見直線
阿新 • • 發佈:2018-11-09
【題目】
題目描述:
在 直角座標平面上有 條直線 ,若在y值為正無窮大處往下看,能見到Li的某個子線段,則稱Li為可見的,否則Li為被覆蓋的。
例如,對於直線 : ; : ; : ,則 和 是可見的, 是被覆蓋的。給出 條直線,表示成 的形式( ≤ ),且 條直線兩兩不重合。求出所有可見的直線.
輸入格式:
第一行為 ,接下來的 行輸入
輸出格式:
從小到大輸出可見直線的編號,兩兩中間用空格隔開,最後一個數字後面也必須有個空格
樣例資料:
輸入
3
-1 0
1 0
0 0
輸出
1 2
【分析】
其實看這道題的題意,應該還是很容易可以想到半平面交吧
還是按照斜率大小排序,每次加新線的時候就比較一下交點位置,然後判斷是否需要彈掉棧頂
不過這道題比半平面交模板不同的是不用彈掉隊首和隊尾(因為它們肯定能被看見)
注意這道題有平行線的情況,我們只用保留
值最大的那條
【程式碼】
#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;
}