Hdu 5033(单调栈+离线处理)

Hdu 5033

思路:

这个题乍一看和可以用二分+查询处理,但是时间会超时。

所以需要先离线处理一下,将所有的查询视为高度为0的点,要求出一个人的最大视野就是找到最影响

他视线的点在哪里,用单调栈维护一个递减的斜率,当下一次斜率不递减时就说明这个高度是最影响视野的高度,

然后分别从左,从右遍历,求出左边的最影响的点,和右边最影响的点,然后更新就得到了这个点的最大视野,

就可以求出最大的角度了。

注意:斜率求解不能用绝对值,因为那样无法判断出反向的斜率,就是如果斜率由正变为负(带绝对值)就不能判断出正负。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 300300;
const double PI = acos(-1.0);
struct Node{
	double x,y;
	int id;
}cur[maxn];
double L[maxn],R[maxn];
int st[maxn];
bool cmp(Node a,Node b){
	return a.x<b.x;
}
double f(int i,int j,int fg){
	if(fg==1) return (cur[i].y-cur[j].y)/(cur[j].x-cur[i].x);
	return (cur[j].y-cur[i].y)/(cur[j].x-cur[i].x);  
}
int main(void)
{
	int n,T,i,j,pt=1,q;
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n);
		for(i=1;i<=n;i++) scanf("%lf%lf",&cur[i].x,&cur[i].y),cur[i].id = -1;
		scanf("%d",&q);
		for(j=1;j<=q;j++){
			scanf("%lf",&cur[j+n].x);
			cur[j+n].y = 0;
			cur[j+n].id = j;
		}
		n+=q;
		sort(cur+1,cur+1+n,cmp); //按照x坐标排序 
		memset(L,0,sizeof(L));
		memset(R,0,sizeof(R));
		int cnt = 0;
		//更新左边界角度 
		for(i=1;i<=n;i++){
			while(cnt>=2&&f(st[cnt-1],st[cnt],1)>=f(st[cnt],i,1)) cnt--;
			if(cur[i].id>=0){
				if(cnt==0) L[cur[i].id] = 0;
				else L[cur[i].id] = atan(f(st[cnt],i,1))/(2*PI)*360.0;
			}
			st[++cnt] = i;
		}
		cnt = 0;
		//更新右边界角度 
		for(i=n;i>=1;i--){
			while(cnt>=2&&f(st[cnt],st[cnt-1],2)>=f(i,st[cnt],2)) cnt--;
			if(cur[i].id>=0){
				if(cnt==0) R[cur[i].id] = 0;
				else R[cur[i].id] = atan(f(i,st[cnt],2))/(2*PI)*360.0;
			}
			st[++cnt] = i;
		}
		printf("Case #%d:\n",pt++);
		for(i=1;i<=q;i++){
			printf("%.10lf\n",180.0-(L[i]+R[i]));
			//for(i=1;i<=q;i++) printf("%lf %lf\n",L[i],R[i]);
		}
	}
	return 0;
}

猜你喜欢

转载自blog.csdn.net/qq_41829060/article/details/89367348