bzoj。 1007: [HNOI2008]水平可见直线

1007: [HNOI2008]水平可见直线

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 8043  Solved: 3060
[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

类半平面交的题。。。先极角排序，因为k给出了直接用k排序，然后根据斜率，画图看一看。。。

后序的形成的点如果在前序点的左边，就会把之前的线给覆盖掉。。。

#include <iostream>#include <stdio.h>#include <string.h>#include <queue>#include <vector>#include <cmath>#include <algorithm>using namespace std;const int maxn = 51550;const double esp = 1e-9;int n, pn,  top;vector<int> ans;struct Line {    int k, b, id;    bool operator < (const Line& l) const {        return k < l.k || (k == l.k && b > l.b);    }    bool operator == (const Line& l) const {        return k == l.k;    }}l[maxn], dq[maxn];double get_x(const Line& l1, const Line& l2) {    return (l2.b - l1.b) * 1.0/ (l1.k - l2.k);}void solve() {    sort(l + 1, l + n + 1);    n = int(unique(l + 1, l + n + 1) - l - 1);    top = 1;    for (int i = 1; i <= n; ++i) {        while (top > 2 && (get_x(dq[top - 1], dq[top - 2]) > (get_x(dq[top - 1], l[i])) - esp))            top--;        dq[top++] = l[i];    }       for (int i = 1; i < top; ++i) {        ans.push_back(dq[i].id);    }    sort(ans.begin(), ans.end());    for (int i = 0; i < ans.size(); ++i) {        cout << ans[i] << " ";    }    //cout << endl;}int main() {    //freopen("in.txt", "r", stdin);    cin >> n;    int k, b;    for (int i = 1; i <= n; ++i) {        cin >> k >> b;        l[i] = Line{k, b, i};    }    solve();}

啊蠢了一晚上。