// hdu5033 Building 单调队列
//
// 题目大意:
//
// n栋大楼,有一个高度h和位置x,现在有一个人高度为0,有q个询问
// 每个询问有一个位置x,求在位置x能看到天空的最大的角度。
//
// 解题思路:
//
// 首先得想到将q个询问的位置作为一栋大楼放在整个大楼中考虑,这样
// 问题就比较一致,也比较容易处理啦。
//
// 想象一下,每次从左往右走,对于当前位置的左半边的90度内,看天的角度是
// 逐渐减小的,这样,如果当前位置的左边比当前位置要低,那么当前位置就挡住了
// 之前的更矮的,也就是说之后的位置是不可能看到比当前更矮的大楼了,这样,就可
// 以看成是一个单调的序列,序列中元素高度单调递减,即所谓的单调队列。但仅仅只是
// 这样是不够的,也会存在比当前元素(设为h)高设为h1,在队列中h1前面的元素是h2,虽然
// h1是比h大,但是如果顺着h1的高度看过去,会被h2挡住,这时h1是没有作用的,也要将
// 这个元素从队列中移除。
//
// 当遇到的大楼要查询的位置的时候,此时队列最后的一个元素并不一定是最优的值,因为
// 这样的值可能是不合法的,比如上面的h1会被h2挡住,并不会看到天空。照样移除这样的元素
// 最后,队列的最后一个元素就是解。
//
// 右边的问题完全可以转化为左边的问题,只是从右往左处理而已。
//
//
// 感悟:
//
// 这道题是14年北京区的网赛的一道题目,当时的我虽然看懂了题目的意思,但是真的完全不会做。
// 直到现在伟大的MW大咖,说是单调队列可以做,我就做了,但是卡了一天半,还是想不出来具体怎么解
// 首先将询问的位置当作大楼这一点我想到了,维持高度递减的单调队列我想到了,就是最后一种情况没有
// 想到,更不知道要怎么处理,在伟大的MW大咖的敦敦教诲还有耐心的提示下,终于终于想到了解决的办法
// 不容易啊不容易。在此真诚感谢MW大咖~
//
// 过程无疑会有很多的疑惑,一个地方不慎,满盘皆输,找到一处的错误,心里十分的欣喜,特别是在ac
// 之后,那种心情,实在是难以描述,感觉到自己的付出,真的是有回报,以前不懂的东西,自己认真学
// 总会有收获,哪怕收获是那么一点点,微不足道,但收获就是收获,没有这个能让人更加欣喜啦,痛苦并快乐着
// 十分的享受,还是那句话,继续练吧~~~
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <iostream>
#include <cstring>
using namespace std;
const int MAX_N = 2e5 + 9;
const double PI = acos(-1.0);
int n,q;
double deqh[MAX_N];
double deqx[MAX_N];
double angle[MAX_N];
struct node {
double pos;
double h;
int id;
node(){
}
node(double pos,double h,int id):pos(pos),h(h),id(id){
}
};
node sky[MAX_N];
bool cmp(node a,node b){
return a.pos < b.pos;
}
void input(){
scanf("%d",&n);
double x,h;
for (int i = 1;i <= n;i++){
scanf("%lf%lf",&x,&h);
sky[i] = node(x,h,0);
}
scanf("%d",&q);
for (int i = 1;i <= q;i++){
scanf("%lf",&x);
sky[i+n] = node(x,0.0,i);
angle[i] = 0.0;
}
n += q;
sort(sky+1,sky+n+1,cmp);
}
double getk(double a,double b){
return fabs(a/b);
}
void getleft(){
int head,tail;
head = tail = 0;
for (int i=1;i<=n;i++){
if (sky[i].id){
while(head + 1 < tail){
double k1 = getk(deqh[tail-1],deqx[tail-1]-sky[i].pos);
double k2 = getk(deqh[tail-2],deqx[tail-2]-sky[i].pos);
if (k1 <= k2)
tail--;
else break;
}
angle[sky[i].id] += atan(deqh[tail-1]/fabs(deqx[tail-1]-sky[i].pos));
}else {
while(head < tail && deqh[tail-1] < sky[i].h)
tail--;
while(head + 1 < tail){
double k1 = getk(deqh[tail-1]-sky[i].h,deqx[tail-1]-sky[i].pos);
double k2 = getk(deqh[tail-2]-sky[i].h,deqx[tail-2]-sky[i].pos);
if (k1 <= k2)
tail--;
else
break;
}
deqh[tail] = sky[i].h;
deqx[tail++] = sky[i].pos;
}
}
}
void getright(){
int head,tail;
head = tail = 0;
for (int i = n;i >= 1;i--){
if (sky[i].id){
while(head + 1 < tail){
double k1 = getk(deqh[tail-1],deqx[tail-1]-sky[i].pos);
double k2 = getk(deqh[tail-2],deqx[tail-2]-sky[i].pos);
if (k1 <= k2)
tail--;
else
break;
}
angle[sky[i].id] += atan(deqh[tail-1]/fabs(deqx[tail-1]-sky[i].pos));
}else {
while(head < tail && deqh[tail-1] < sky[i].h)
tail--;
while(head + 1 < tail){
double k1 = getk(deqh[tail-1]-sky[i].h,deqx[tail-1] - sky[i].pos);
double k2 = getk(deqh[tail-2]-sky[i].h,deqx[tail-2] - sky[i].pos);
if (k1 <= k2){
tail--;
}else
break;
}
deqh[tail] = sky[i].h;
deqx[tail++] = sky[i].pos;
}
}
}
void solve(){
getleft();
getright();
for (int i=1;i<=q;i++){
printf("%.10lf\n",(PI - angle[i]) * 180.0 / PI);
}
}
int main(){
int t;
//freopen("1.txt","r",stdin);
scanf("%d",&t);
int kase = 1;
while(t--){
input();
printf("Case #%d:\n",kase++);
solve();
}
}
hdu5033 Building 单调队列
最新推荐文章于 2019-04-17 23:26:01 发布