bzoj1007: [HNOI2008]水平可见直线 关于一类半半平面交题目的解法研究

bzoj1007: [HNOI2008]水平可见直线

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

知识点：半半平面交

这是一半的半平面交，只需要维护下凸壳。而这类型的题目也有一个比较水的写法。
当然是先极角排序。考虑下面这种情况

你瞧绿色和红色的直线盖住了蓝色的。
你瞧绿色和蓝色的交点在红色和蓝色的交点前面。
这个时候蓝色的线就没用了，如果这个时候蓝色的直线在栈顶，可以直接弹掉。
否则的话就可以直接入栈了。
也就是维护下凸壳的交点横坐标单调。
可以自己画画图，一定是正确的（我不会证，水过去水过去）

代码

/**************************************************************
    Problem: 1007
    User: 2014lvzelong
    Language: C++
    Result: Accepted
    Time:364 ms
    Memory:2584 kb
****************************************************************/

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
const int maxn = 55000;
const double eps = 1e-8;
using namespace std;
int n;
int read()
{
    char ch = getchar(); int x = 0, f = 1;
    while(ch < '0' || ch > '9') {if(ch == '-') f = -1; ch = getchar();}
    while(ch >= '0' && ch <= '9') {x = x * 10 + ch - '0'; ch = getchar();}
    return x * f;
}

int dblcmp(double a) {
    if(abs(a) < eps) return 0;
    return a > 0 ? 1 : -1;
}

struct Line {
    double a, b; int id;
}line[maxn];
int st[maxn], top;
bool cmp1(Line a, Line b) {return !dblcmp(a.a - b.a) ? dblcmp(a.b - b.b) > 0 : a.a > b.a;}
double inter(Line a, Line b) {return (b.b - a.b) / (a.a - b.a);}

int main()
{
    n = read();
    for(int i = 1;i <= n; ++i) {
        scanf("%lf%lf", &line[i].a, &line[i].b);
        line[i].id = i;
    }
    sort(line + 1, line + n + 1, cmp1); st[top = 1] = 1;
    for(int i = 2;i <= n; ++i) 
    {
        if(dblcmp(line[i].a - line[st[top]].a)) {
            while(top >= 2) {
                double x1 = inter(line[st[top]], line[i]);
                double x2 = inter(line[st[top - 1]], line[i]);
                int p = dblcmp(x1 - x2);
                if(p >= 0) --top;
                else break;
            }
            st[++top] = i;
        }
    }
    for(int i = 1;i <= top; ++i) st[i] = line[st[i]].id;
    sort(st + 1, st + top + 1);
    for(int i = 1;i <= top; ++i) printf("%d ", st[i]);
    puts("");
    return 0;
}