BZOJ 1007: [HNOI2008]水平可见直线

Posted czy020202

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ 1007: [HNOI2008]水平可见直线相关的知识,希望对你有一定的参考价值。

1007: [HNOI2008]水平可见直线

Time Limit: 1 Sec  Memory Limit: 162 MB

Submit: 7769  Solved: 2973

[Submit][Status][Discuss]

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

题解

按斜率排序后,从小到大枚举,用栈维护当前可见的直线,如果当前考虑的直线和栈顶直线的交点在栈顶直线与栈顶第二条直线的交点的左侧,那么栈顶直线一定看不见,弹出。

因为要用栈顶第二条直线,用手写栈即可。

注意处理精度问题。

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
using namespace std;
const int N=500005;
const double eps=1e-8;
int n,top;
struct line{
	int id;
	double a,b;
}l[N],st[N];
bool cmp1(line a,line b){
	if(fabs(a.a-b.a)<=eps)return a.b<b.b;
	return a.a<b.a;
}
bool cmp2(line a,line b){
	return a.id<b.id;
}
double cross(line a,line b){
	return (b.b-a.b)/(a.a-b.a);
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		scanf("%lf%lf",&l[i].a,&l[i].b);
		l[i].id=i;
	}
	sort(l+1,l+n+1,cmp1);
	for(int i=1;i<=n;i++){
		while(top){
			if(fabs(l[i].a-st[top].a)<=eps)top--;
			else if(top>1&&cross(l[i],st[top-1])<=cross(st[top],st[top-1]))top--;
			else break;
		}
		st[++top]=l[i];
	} 
	sort(st+1,st+top+1,cmp2);
	for(int i=1;i<=top;i++){
		printf("%d ",st[i].id);
	}
	printf("\n");
	return 0;
}

以上是关于BZOJ 1007: [HNOI2008]水平可见直线的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ 1007: [HNOI2008]水平可见直线

bzoj 1007 [HNOI2008]水平可见直线(单调栈)

BZOJ 1007: [HNOI2008]水平可见直线

bzoj 1007: [HNOI2008]水平可见直线

bzoj 1007: [HNOI2008]水平可见直线

bzoj1007: [HNOI2008]水平可见直线