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

[bzoj1007] [HNOI2008]水平可見直線

pac output const clu \n 水平 n) bit 一個數

Description

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

Input

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

Output

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

Sample Input

3
-1 0
1 0
0 0

Sample Output

1 2

Solution

用棧維護下凸殼即可。

#include<bits/stdc++.h>
using namespace std;

void read(int &x) {
    x=0;int f=1;char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-f;
    for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';x*=f;
}

#define write(x) printf("%d\n",x)

const int maxn = 2e5+10;

struct line {int k,b,id;}l[maxn],sta[maxn];

int cmp1(line a,line b) {return a.k==b.k?a.b<b.b:a.k<b.k;}

int cmp2(line a,line b) {return a.id<b.id;}

double inter(line a,line b) {return (double)(a.b-b.b)/(b.k-a.k);}

int main() {
    int n;read(n);
    for(int i=1;i<=n;i++) read(l[i].k),read(l[i].b),l[i].id=i;
    sort(l+1,l+n+1,cmp1);int top=0;
    for(int i=1;i<=n;i++) {
        while(top) {
            if(sta[top].k==l[i].k) top--;
            else if(top>1&&inter(l[i],sta[top-1])<=inter(sta[top],sta[top-1])) top--;
            else break; 
        }sta[++top]=l[i];
    }
    sort(sta+1,sta+top+1,cmp2);
    for(int i=1;i<=top;i++) printf("%d ",sta[i].id);puts("");
    return 0;
}

[bzoj1007] [HNOI2008]水平可見直線