水平可见直线

【题目描述】
在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条直线两两不重合.求出所有可见的直线。
【输入格式】
第一行为N(0 < N < 50000)
接下来的N行输入Ai,Bi
【输出格式】
从小到大输出可见直线的编号，两两中间用空格隔开,最后一个数字后面也必须有个空格
【样例输入】
3
-1 0
1 0
0 0
【样例输出】
1 2
【分析】
如果只有1条直线就太简单了，肯定可以看见。
如果有2条直线，也不难解决，如果两直线平行就遮住了一条，反之都能看见。
接下来的分析不考虑平行。
如果有3条直线，会组成一个三角形，不难发现三角形底边所在直线被遮住了。于是可以得出结论：若斜率最大的直线与斜率次大（或次次大）的直线的交点在斜率次大和斜率次次大的直线的交点的左边，则斜率次大的直线被覆盖。

于是可以得出算法：先把所有直线按斜率排序，然后枚举斜率次次大的直线，分别求出与斜率最大和次大的直线的交点，看交点的位置。
可以用栈来维护，有覆盖的就直接弹栈。最后栈里存放的直线就是没有被覆盖的。
至于求交点，可以列方程组：
A1x+B1=y
A2x+B2=y
然后解一下就可以。

#include<cstdio>#include<cmath>#include<algorithm>using namespace std;struct Line {    int k,b,i;//k表示斜率，b表示截距，i表示序号    bool operator < (const Line& rhs) const//重载运算符“<”    {        if (k==rhs.k) return b>rhs.b;        return k<rhs.k;    }} line[50001];int stack[50001]={0},ans[50001];double find_intersection_point(int i,int j) {    int k1=line[i].k, b1=line[i].b, k2=line[j].k, b2=line[j].b;    return (double)(b2-b1)/(k1-k2);//求交点，方程不解释}int main(){    int n,top=1;    scanf("%d",&n);    if (n==1) { printf("1 \n");return 0; }//特判    for (int i=0;i<n;i++) {        int k,b;        scanf("%d%d",&k,&b);        line[i]=(Line){k,b,i+1};    }    sort(line,line+n);//排序    //其实这里少了一步stack[0]=0，即斜率最大的直线先进栈，但是不影响结果    if (line[1].k!=line[0].k) stack[top++]=1;//若两直线不平行则斜率次大的进栈    for (int i=2;i<n;i++) {        if (line[i].k==line[i-1].k) continue;//平行肯定覆盖了        while (top>1 && find_intersection_point(i,stack[top-1])-find_intersection_point(stack[top-2],stack[top-1])<=1e-6) top--;//判断交点        stack[top++]=i;//入栈    }    for(int i=0;i<top;i++) ans[i]=line[stack[i]].i;//找到栈中的直线编号    sort(ans,ans+top);//栈里的直线不一定是按编号递增排列的    for(int i=0;i<top;i++) printf("%d ",ans[i]);    return 0;}