BZOJ 1007 水平可见直线

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

两种做法吧：

1. 直接半平面交暴搞，虽然我不会写；
2. 第二种主要是针对这个题目的： 利用题所隐藏的性质。稍微YY一下，所能看到的直线沿着X正方向斜率一定是单调递增的。

　　首先将斜率排序(从小到大)之后一次加入栈中。对于栈顶直线l1,次顶直线l2，以及所枚举到的直线i,如果l与l1的交点在l1与l2交点左方，栈顶的弹出(证明：画画图就知道了)。

　　最后栈中的直线即为答案。(这个好像就是半平面交)

代码如下：

 #include<cstdlib>
 #include<cstdio>
 #include<vector>
 #include<algorithm>
 using namespace std;
 
 #define esp (1e-6)
 #define maxn 50010
 int n,S[maxn],top;
 vector <int> vec;
 struct Line{double A,B; int ord;}line[maxn];
 
 inline bool cmp(Line a,Line b)
 {
     if (a.A == b.A) return a.B > b.B;
     return a.A < b.A;
 }
 
 inline double calc(Line a,Line b) {return (double)(b.B-a.B)/(double)(a.A-b.A);}
 
 inline bool okay(double a,double b)
 {
     if (a + esp > b &&a-esp<b) return true;
     if (a > b) return true;
     return false;
 }
 
 inline void work()
 {
     int i;
     S[++top] = ; i = ;
     while (i <= n&&line[i].A == line[i-].A) i++;
     if (i <= n)
     {
         S[++top] = i; ++i;
         while (i <= n&&line[i].A == line[i-].A) i++;
         for (;i <= n;++i)
         {
             if (line[S[top]].A == line[i].A)
                 continue;
             else
             {
                 while (top >  && okay(calc(line[S[top]],line[S[top-]]),calc(line[S[top]],line[i])))
                     --top;
                 S[++top] = i;
             }
         }
     }
     for (i = ;i <= top;++i) vec.push_back(line[S[i]].ord);
     sort(vec.begin(),vec.end());
 }
 
 int main()
 {
     freopen("1007.in","r",stdin);
     freopen("1007.out","w",stdout);
     scanf("%d",&n); int i;
     for (i = ;i <= n;++i)
     {
         scanf("%lf %lf",&line[i].A,&line[i].B);
         line[i].ord = i;
     }
     sort(line+,line+n+,cmp);
     work();
     int nn = vec.size();
     for (i = ;i < nn;++i) printf("%d ",vec[i]);
     fclose(stdin); fclose(stdout);
     return ;
 }