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
HINT

..打的时候各种小错误。。我的正确率啊！！

。

分析：肯定是先按照斜率从小到大排序，然后按照截距从大到小，因为斜率相同，截距小的肯定被覆盖了。
然后先把前两条线加入队列，每次枚举的线i，如果i与s[top]的交点在s[top-1]和s[top]交点的左边或者重合，s[top]肯定被覆盖了，然后把栈顶元素弹出。如此往复，最后栈中的元素就是最终的答案。

为什么是这么做？因为s[top]与s[top-1]的交点左边，s[top]已经被s[top-1]覆盖了，如果i和s[top]的交点还在这个点左边，那么这个点右边的s[top]都会被遮住。

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
int n,m;
const int N=1e6+6;
struct node
{
    double k,b;
    int id;
}a[N];
int s[N];
bool cmp(node a,node b)
{
    if (a.k==b.k)return a.b>b.b;
    else return a.k<b.k;
}
inline double getpos(node a,node b)
{
    return (b.b-a.b)/(a.k-b.k);
}
bool ans[N];
int main()
{
    scanf("%d",&n);
    fo(i,1,n)
    {
        scanf("%lf%lf",&a[i].k,&a[i].b);
        a[i].id=i;
    }
    sort(a+1,a+1+n,cmp);
    int top=1;
    s[1]=1;
    fo(i,2,n)
    {
        if (a[i].k-a[i-1].k< 1e-8)continue;
        while(top>1&&getpos(a[i],a[s[top]])<=getpos(a[s[top]],a[s[top-1]]))top--;
        s[++top]=i;
    }
    fo(i,1,top)ans[a[s[i]].id]=1;
    fo(i,1,n)
    if (ans[i])printf("%d ",i);
    return 0;

}