bzoj1007/luogu3194 水平可见直线 (单调栈)

先按斜率从小到大排序,然后如果排在后面的点B和前面的点A的交点是P,那B会把A在P的右半段覆盖掉,A会把B在P的左半段覆盖掉。

然后如果我们现在又进来了一条线,它跟上一条的交点还在上一条和上上条的左边,这就说明上一条完全被覆盖了

这样的话,维护一个单调栈做一做就可以了

(要先处理一下,斜率相同的只留下B最大的,而且会有重合的线,都要输出)

 1 #include<bits/stdc++.h>
 2 #define pa pair<int,int>
 3 #define ll long long
 4 using namespace std;
 5 const int maxn=50050;
 6
 7 ll rd(){
 8     ll x=0;char c=getchar();int neg=1;
 9     while(c<‘0‘||c>‘9‘){if(c==‘-‘) neg=-1;c=getchar();}
10     while(c>=‘0‘&&c<=‘9‘) x=x*10+c-‘0‘,c=getchar();
11     return x*neg;
12 }
13
14 struct Line{
15     int a,b,i;
16 }p[maxn],q[maxn];
17 int N,NN,stk[maxn],sh;
18 bool flag[maxn];
19
20 inline bool cmp(Line a,Line b){
21     return a.a==b.a?a.b>b.b:a.a<b.a;
22 }
23
24 inline double getx(Line a,Line b){
25     return (double)(b.b-a.b)/(a.a-b.a);
26 }
27 inline bool nsame(int a,int b){
28     return q[a].a!=q[b].a||q[a].b!=q[b].b;
29 }
30
31 int main(){
32     int i,j,k;
33     N=rd();for(i=1;i<=N;i++){
34         p[i].a=rd(),p[i].b=rd();p[i].i=i;
35     }
36     sort(p+1,p+N+1,cmp);
37     for(i=1,j=0;i<=N;i++){
38         if(p[i].a==p[i-1].a&&p[i].b<p[i-1].b) continue;
39         q[++j]=p[i];
40     }NN=j;
41     stk[sh=1]=1;
42     for(i=2;i<=NN;i++){
43         while(sh>=2&&nsame(stk[sh],i)&&(getx(q[i],q[stk[sh]])<=getx(q[stk[sh]],q[stk[sh-1]]))) sh--;
44         stk[++sh]=i;
45     }for(i=sh;i;i--) flag[q[stk[i]].i]=1;
46     for(i=1;i<=N;i++) if(flag[i]) printf("%d ",i);
47     return 0;
48 }

原文地址:https://www.cnblogs.com/Ressed/p/9634188.html

时间: 2024-11-09 10:14:56

bzoj1007/luogu3194 水平可见直线 (单调栈)的相关文章

bzoj1007: [HNOI2008]水平可见直线(单调栈)

1007: [HNOI2008]水平可见直线 题目:传送门 题解: 蒟蒻在bzoj上做的第一道计算几何 其实这道题并不难...(所以我A了) 仔细想想不难发现,其实我们只需要维护一个下凸的图形... 只有在这个图形上的直线才不会被覆盖,也就是可以被上帝直线看到的孩子 为什么呢...自己画个图模拟吧. 那么具体的做法我们就要采用单调栈啦! 把给出的直线按照斜率从小到大排序(如果斜率相同的话就按照b来排) 然后一个一个放入我们强大的单调栈中,在加入的同时我们当然还要进行维护: 如果栈顶的直线与新加直

【BZOJ1007】水平可见直线(单调栈)

[BZOJ1007]水平可见直线(单调栈) 题解 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 < 50

BZOJ 1007 HNOI 2008 水平可见直线 计算几何+栈

题目大意:给出一些笛卡尔系中的一些直线,问从(0,+∞)向下看时能看到哪些直线. 思路:半平面交可做,但是显然用不上.类似于求凸包的思想,维护一个栈.先将所有直线按照k值排序,然后挨个压进去,遇到有前一个交点被挡住的话就先弹栈. 比较闹心的是去重.我的方法是压栈之前先去重,然后在处理. CODE: #include <cmath> #include <cstdio> #include <cstring> #include <iostream> #includ

[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

[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 从

[BZOJ1007][HNOI2008]水平可见直线 计算几何

题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1007 可以发现题目求的就是一个下凸包,把直线按斜率排序,再来维护凸包就好了.可以发现下凸包上的拐点横坐标单增.同时注意处理斜率相同的直线的情况. 1 #include<cstdio> 2 #include<cstring> 3 #include<algorithm> 4 using namespace std; 5 const double eps=1e-14;

[日常摸鱼]bzoj1007[HNOI2008]水平可见直线-半平面交(对偶转凸包)

不会写半平面交-然后发现可以转成对偶凸包问题 具体见这里:http://trinkle.blog.uoj.ac/blog/235 相关的原理我好像还是不太懂-orz #include<cstdio> #include<algorithm> const int N=50005; inline int read() { int s=0,f=1;char c=getchar(); while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}

【bzoj1007】[HNOI2008]水平可见直线 半平面交/单调栈

题目描述 在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 输出 从小到大输出可见直线的编号,两两中间

BZOJ_1007_ [HNOI2008]_水平可见直线_(单调栈+凸包)

描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1007 给出一些直线,沿着y轴从上往下看,能看到多少条直线. 分析 由于直线相交,会遮挡住一些直线. 自己画画图就可以发现,最后能看见的直线,也就是在最上面的那些直线一定构成一个凸包的下凸壳(没错一定是凸的). 接下来就是如何求这个下凸壳了. 先按照斜率为第一关键字,截距为第二关键字,将直线从小到大排序.用一个斜率单调递增的栈来维护凸壳. 我们按照排序后的顺序添加直线,画画图会发现: 1.斜率