hdu 5033 Building 2014 ACM/ICPC Asia Regional Beijing Online

题目链接:hdu 5033

        给定一系列的建筑的高度的位置,然后给定m个位置求能看到的最大仰角

        将建筑与人按位置排序,对每个一点记录与前一个点构成的角度,维护一个角度递增的的序列。然后每个人最大的仰角的位置就是在维护的序列中角度最大且小于该点与人构成的角度所在位置。然后利用三角函数计算即可

/******************************************************
 * File Name:   5033.cpp
 * Author:      kojimai
 * Creater Time:2014年09月23日 星期二 22时11分06秒
******************************************************/

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
using namespace std;
#define FFF 100005
int num,now;
const double pi=atan(1.0)*4;
struct node
{
	double h,ang,loc;
	int num;
}p[FFF*2],stack[FFF];
bool cmp(node a,node b)
{
	return a.loc<b.loc;
}
double ans[FFF];
int qfind(node x,double ju)
{
	int l=0,r=now-1,mid;
	while(l<=r)
	{
		mid=(l+r)/2;
		if(stack[mid].ang<=(stack[mid].h-x.h)/(x.loc-stack[mid].loc)*ju)
		{
			l=mid+1;
		}
		else
		{
			r=mid-1;
		}
	}
	return l;
}
void solve(node x,double ju)
{
	int t=qfind(x,ju);
	now = t+1;
	stack[t]=x;
	if(t==0)
		stack[t].ang=0;
	else
	{
		stack[t].ang=(stack[t-1].h-stack[t].h)/(stack[t].loc-stack[t-1].loc)*ju;
	}
	return;
}
void add(int num,node x,double ju)
{
	int t=qfind(x,ju);
	t--;
	double tt=1.0/(stack[t].h/(x.loc-stack[t].loc)*ju);
	//cout<<endl<<"t="<<t<<" tan="<<tt<<endl;
	//cout<<"x.loc="<<x.loc<<" stack[t].loc="<<stack[t].loc<<" stack[t].h="<<stack[t].h<<endl;
	//cout<<"atan="<<atan(tt)*180.0/pi<<endl<<endl;;
	ans[num]+=atan(tt)*180.0/pi;
}
void out()
{
	cout<<" now="<<now<<endl;
	for(int i=0;i<now;i++)
	{
		cout<<i<<" ang="<<stack[i].ang<<endl;
	}
}

int main()
{
	int m,keng,n,ca=1;
	scanf("%d",&keng);
	while(keng--)
	{
		scanf("%d",&n);
		for(int i=0;i<n;i++)
		{
			scanf("%lf%lf",&p[i].loc,&p[i].h);
			p[i].num=-1;
		}
		scanf("%d",&m);
		for(int i=0;i<m;i++)
		{
			scanf("%lf",&p[i+n].loc);
			p[i+n].h=0;
			p[i+n].num=i;
		}
		sort(p,p+n+m,cmp);
		p[0].ang=0;
		stack[0]=p[0];
		now=1;
		memset(ans,0,sizeof(ans));
		for(int i=1;i<n+m;i++)
		{
			if(p[i].num==-1)
			{
				solve(p[i],1.0);
			}
			else
			{
				add(p[i].num,p[i],1.0);
			}
		}
		stack[0]=p[n+m-1];
		now=1;
		for(int i=n+m-1;i>=0;i--)
		{
			if(p[i].num==-1)
				solve(p[i],-1.0);
			else
				add(p[i].num,p[i],-1.0);
		}
		printf("Case #%d:\n",ca++);
		for(int i=0;i<m;i++)
			printf("%.10f\n",ans[i]);
	}
	return 0;
}