BZOJ 1007: [HNOI2008]水平可见直线 栈/计算几何

1007: [HNOI2008]水平可见直线

Time Limit: 1 Sec  Memory Limit: 162 MB

题目连接

http://www.lydsy.com/JudgeOnline/problem.php?id=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

题解：

首先我们按照斜率从大到小排序，然后我们对一个堆进行优化

如果要插入一条直线的话，K最大的那条线和K最小的肯定不会背盖住，只会有中间的那条直线被压住

所以我们就判定一下，然后不停的更新就好了

具体判定是看 斜率小的直线与斜率中间的直线，斜率小的直线和斜率大的直线，这两个x的坐标大小进行比较，然后就可以啦~

代码：

//qscqesze
#include <cstdio>
#include <cmath>
#include <cstring>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <typeinfo>
#include <fstream>
#include <map>
typedef long long ll;
using namespace std;
//freopen("D.in","r",stdin);
//freopen("D.out","w",stdout);
#define sspeed ios_base::sync_with_stdio(0);cin.tie(0)
#define maxn 50001
#define mod 10007
#define eps 1e-9
//const int inf=0x7fffffff;   //无限大
const int inf=0x3f3f3f3f;
/*

*/
//**************************************************************************************
struct node
{
    double x,y;
    int id;
};
bool cmp(node c,node d)
{
    return c.x>d.x;
}
double kiss(node a,node b)
{
    return (b.y-a.y+0.0)/(a.x-b.x+0.0);
}
node a[maxn],aa[maxn];
int s[maxn];
inline ll read()
{
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
bool cmp2(int c,int d)
{
    return a[c].id<a[d].id;
}
int main()
{
    int n;
    cin>>n;
    for(int i=;i<=n;i++)
    {
        cin>>aa[i].x>>aa[i].y;
        aa[i].id=i;
    }
    sort(aa+,aa+n+,cmp);
    int M=;
    for(int i=;i<=n;i++)
    {
        if(aa[i].x!=a[i-].x)
            a[++M]=aa[i];
        else if(aa[i].y>a[M].y)
            a[M].y=aa[i].y,a[M].id=aa[i].id;
    }
    int top=;
    s[]=;top=;
    for(int i=;i<=M;i++)
    {
        while(top>=)
        {
            double x1=kiss(a[s[top-]],a[i]);
            double x2=kiss(a[s[top]],a[i]);
            if(x1<=x2+1e-)
                top--;
            else
                break;
        }
        s[++top]=i;
    }
    sort(s+,s+top+,cmp2);
    for(int i=;i<=top;i++)
        cout<<a[s[i]].id<<" ";
    return ;
}