bzoj 1007 [HNOI2008]水平可见直线（单调栈）

 

1007: [HNOI2008]水平可见直线

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 5120  Solved: 1899
[Submit][Status][Discuss]

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

HINT

Source

 

【思路】

       单调栈维护下凸包。

 

【代码】

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define FOR(a,b,c) for(int a=(b);a<=(c);a++)
using namespace std;

const int N = 50000+10;
const double eps = 1e-8;

struct Line {
    double a,b; int r;
    bool operator < (const Line& rhs) const {
        if(fabs(a-rhs.a)<eps) return b<rhs.b;
        else return a<rhs.a;
    }
}L[N],S[N];

double cross(Line x1,Line x2) {
    return (x2.b-x1.b)/(x1.a-x2.a);
}

int n,flag[N],top;

int main() {
    scanf("%d",&n);
    FOR(i,1,n) {
        scanf("%lf%lf",&L[i].a,&L[i].b);
        L[i].r=i;
    }
    sort(L+1,L+n+1);
    FOR(i,1,n) {
        while(top) {
            if(fabs(S[top].a-L[i].a)<eps) top--;
            else if(top>1 && cross(L[i],S[top-1])<=cross(S[top],S[top-1])) top--;
            else break;
        }
        S[++top]=L[i];
    }
    FOR(i,1,top) flag[S[i].r]=1;
    FOR(i,1,n) if(flag[i])
        printf("%d ",i);
    return 0;
}