hdu 5033 buiding(单调栈)

某年某月某天,马特去了一个小镇。这个小镇如此狭窄,以至于他可以把小镇当作一个枢纽。在镇上有一些摩天大楼,其中一栋位于xi,高度为hi。所有的摩天大楼位于不同的地方。为了简化题目,假设摩天大楼没有宽度。由于摩天大楼如此之高,马特几乎看不到天空。对于马特所在的位置,他想知道他能看到天空的角度范围有多大。假设马特的身高是0。可以保证,对于每个查询,马特的左右两边至少有一座建筑物,而且他的位置上没有建筑物。建筑物的数量n<1e5,查询次数Q<1e5。

我用常规思路,想不出什么好方法。既然这道题的普遍解法是单调栈,那么我们就来看一下单调栈怎么实现。首先,单调队列和单调栈的前提是,加入的元素满足一种有序的关系。所以我们自然而然地就想到,将所有查询排序,从左向右和从右向左分别进行查询的计算。以从左向右为例,如何做到在(均摊)O1的时间内找到一个查询的答案呢?我们先来找些性质。

性质1:如果一个建筑物i比它左侧的j高,那么当i和j都在马特左侧时,马特一定看不到j。也就是说,如果维护一个关于建筑物下标的单调栈,新加进来一个建筑时,可以把它左侧所有更矮的建筑删掉。那么,现在的单调栈就是单调递减的。这就去除了第一种冗余状态。

性质2:从左向右而言,如果有这样的情况:

,也就是单调栈中的顶上两个建筑和新加进来的建筑组成一个凹包(不好意思,下凸包),那么可以直接把单调栈顶的建筑删了,一直循环删下去。因为无论马特怎么站,都不会被它挡住。第二种冗余状态也被我们去除了。那么现在,单调栈存储的建筑就组成了一个上凸包。

性质3:然而,优化力度还不够大。这样子算法依然不是O(n)的。那怎么办呢?

由于马特一直往右跑,那么他与凸包相切的点只会往左移动:

所以,在加入新的建筑之后,判断一下马特与当前凸包的交点,然后将交点右侧的建筑全删了(弹出)即可。栈顶的建筑就是最优建筑(切线嘛)。这就去除了第三种冗余状态。然而这样的时间复杂度是均摊O(1)的吗?额,显然的。

这可能是目前为止我写的最详细的一篇博文。。

#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std; const int maxn=1e5+5;
const double pi=3.1415926535898;
struct stack{
int t, a[maxn];
void push(int x){ a[++t]=x; }
void pop(){ if (--t==-1) ++t; }
void reset(){ t=0; }
inline int top(){ return a[t]; }
inline int top2(){ return a[t-1]; }
}s;
struct query{
int id, ans1, ans2; double x;
}q[maxn], q2[maxn];
struct building{
double x, h;
}b[maxn]; int n, T, Q; bool cmp1(query &a, query &b){
return a.x<b.x; }
bool cmp2(query &a, query &b){
return a.id<b.id; }
bool cmpb(building &a, building &b){
return a.x<b.x; } double angle(double x, double y){
return atan(x/y)*180/pi;
} int main(){
scanf("%d", &T); int tmp=0;
while (T--){
++tmp; s.reset();
scanf("%d", &n);
for (int i=1; i<=n; ++i)
scanf("%lf%lf", &b[i].x, &b[i].h);
scanf("%d", &Q);
for (int i=1; i<=Q; ++i){
scanf("%lf", &q[i].x);
q[i].id=i;
}
sort(b+1, b+n+1, cmpb);
sort(q+1, q+Q+1, cmp1);
int nowbui=1; s.push(1);
for (int i=1; i<=Q; ++i){
while (b[nowbui+1].x<q[i].x){ //从左到右加入建筑物
++nowbui;
//第一种冗余状态:左边的楼比右边的矮
while (s.t>0&&b[s.top()].h<=b[nowbui].h) s.pop();
//所以现在这个建筑物序列严格下降
//第二种冗余状态:是个下凸函数
while (s.t>1&&(b[nowbui].h-b[s.top()].h)/
(b[nowbui].x-b[s.top()].x)>=
(b[s.top()].h-b[s.top2()].h)/
(b[s.top()].x-b[s.top2()].x)) s.pop();
//所以现在这个建筑物序列是凸包了
s.push(nowbui);
}
//第三种冗余状态:若一栋楼在视线所切的那栋楼的右边,则可以删去
//可以用均摊,证明这里的复杂度是常数
while (s.t>1&&(b[s.top2()].h)/(q[i].x-b[s.top2()].x)>
(b[s.top()].h)/(q[i].x-b[s.top()].x)) s.pop();
//现在终于不仅排除了冗余状态,并且找到了最优解。
q[i].ans1=s.top();
} //从右往左也一样
nowbui=n; s.reset(); s.push(n);
for (int i=Q; i>=1; --i){
while (b[nowbui-1].x>q[i].x){
--nowbui;
//右边的楼比左边的矮
while (s.t>0&&b[s.top()].h<=b[nowbui].h) s.pop();
//下面有图 (其实这里封装一个关于斜率的函数会更好)
while (s.t>1&&(b[nowbui].h-b[s.top()].h)/
(b[nowbui].x-b[s.top()].x)<=
(b[s.top()].h-b[s.top2()].h)/
(b[s.top()].x-b[s.top2()].x)) s.pop();
s.push(nowbui);
}
//下面有图
while (s.t>1&&b[s.top2()].h/(b[s.top2()].x-q[i].x)>
b[s.top()].h/(b[s.top()].x-q[i].x)) s.pop();
q[i].ans2=s.top();
}
sort(q+1, q+Q+1, cmp2);
printf("Case #%d:\n", tmp);
for (int i=1; i<=Q; ++i){
printf("%.5lf\n", 180-angle(b[q[i].ans1].h,
q[i].x-b[q[i].ans1].x)-angle(b[q[i].ans2].h,
b[q[i].ans2].x-q[i].x));
}
}
// printf("%.6lf", angle(1, 1));
return 0;
}

hdu 5033 buiding(单调栈)的更多相关文章

  1. HDU 5033 Building(单调栈)

    HDU 5033 Building(单调栈) 题目链接http://acm.hdu.edu.cn/showproblem.php?pid=5033 Description Once upon a ti ...

  2. hdu 5033 Building (单调栈 或 暴力枚举 )

    Description Once upon a time Matt went to a small town. The town was so small and narrow that he can ...

  3. HDU - 5033 Building (单调栈+倍增)

    题意:有一排建筑,每座建筑有一定的高度,宽度可以忽略,求在某点的平地上能看到天空的最大角度. 网上的做法基本都是离线的...其实这道题是可以在线做的. 对于向右能看到的最大角度,从右往左倍增维护每个时 ...

  4. HDU 5033 Building(单调栈维护凸包)

    盗张图:来自http://blog.csdn.net/xuechelingxiao/article/details/39494433 题目大意:有一排建筑物坐落在一条直线上,每个建筑物都有一定的高度, ...

  5. hdu 5033 模拟+单调优化

    http://acm.hdu.edu.cn/showproblem.php?pid=5033 平面上有n个建筑,每个建筑由(xi,hi)表示,m组询问在某一个点能看到天空的视角范围大小. 维护一个凸包 ...

  6. Largest Rectangle in a Histogram HDU - 1506 (单调栈)

    A histogram is a polygon composed of a sequence of rectangles aligned at a common base line. The rec ...

  7. HDU 1506【单调栈】

    思路: 转化成对于某一位置为最小值求向两边最远>=他的位置,用单调栈就能轻易完成. 那么ans=(left+right)*h[i]; 维护单调递增还是递减呢? 我们能很快反应到,一旦碰到一个比他 ...

  8. HDU 3410【单调栈】

    思路: 单调栈. 鄙人的记忆:按当前为最大值的两边延伸就是维护单调递减栈. //#include <bits/stdc++.h> #include <iostream> #in ...

  9. Hdu 4923(单调栈)

    题目链接 Room and Moor Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Ot ...

随机推荐

  1. SIFT算法详解(转)

    原文地址 http://blog.csdn.net/pi9nc/article/details/23302075 尺度不变特征变换匹配算法详解 Scale Invariant Feature Tran ...

  2. 超限学习机 (Extreme Learning Machine, ELM) 学习笔记 (一)

    1. ELM 是什么 ELM的个人理解: 单隐层的前馈人工神经网络,特别之处在于训练权值的算法: 在单隐层的前馈神经网络中,输入层到隐藏层的权值根据某种分布随机赋予,当我们有了输入层到隐藏层的权值之后 ...

  3. log4j 2 入门实例(3)

    继承机制 所有logger都继承自root logger. 可以认为名为log4j2learn.Hello的logger继承自名为log4j2learn的logger. log4j会先查找名称是&qu ...

  4. python3使用pdfminer3k解析pdf文件

    安装pdfminer模块 pip3 install pdfminer3k 代码如下 #!/usr/bin/env python # coding:utf8 # author:Z time:2018/7 ...

  5. SpringMVC的API和Spring的官方说明文档的地址。

    SpringMVC的API和Spring的官方说明文档的地址. 1.SpringMVC的API的URL: http://docs.spring.io/spring/docs/current/javad ...

  6. 在js实现矩阵转置

    var arr=[[2,4,6,8],[8,9,0,-1],[9,6,2,1]]; //定义一个新的数组 var arr2=[]; for(var i=0;i<arr[0].length;i++ ...

  7. .net中后台c#数组与前台js数组交互

    第一步:定义cs数组  cs文件里后台程序中要有数组,这个数组要定义成公共的数组.  public string[] lat = null;  public string[] lng = null; ...

  8. 分享知识-快乐自己:spring_Boot 中文返回给浏览器乱码 解析成问号?? fastJson jackJson

    心路历程: 在Controller中return 对象的时候,对象中的属性值中文部分在浏览器中 显示为问号?? 然后结果是这样的:?? 尝试排查原因: 中文乱码常有以下三种: 1.request.re ...

  9. Linux下查看端口占用情况

    用启动服务的账号登录,然后运行命令: lsof -i:<端口号> 也可使用命令: netstat -apn|grep <端口号> 找到进程号以后,再使用以下命令查看详细信息: ...

  10. struts2框架xml验证

    struts2验证分为3步: 1.获取需要验证的信息,使用同名属性,提供getter,setter方法.然后框架使用反射将值自动注入. 2.对信息进行验证,成功失败作出对应的选择. xml验证和手动验 ...