hdu 5033 ( Building ) 单调栈

题意： 给出建筑物位置与高度，以及人的位置，询问人看到天空的最大角度

思路：

对于每个位置，人的视野受到左侧与右侧楼与这个点相连的直线  ， 左侧斜率最小的 以及右侧斜率最大的 两条直线限制，如果暴力枚举显然会T

可以处理为一个栈的形式： 

以左侧楼的限制为例：右侧同理

 对于每个位置 如果这个位置是楼，就扔进先除去所有栈中比它高度低的， 再以当前点为基准，所描栈中的点，确保栈中点的到当前点的斜率是递减的（左侧栈是个 凸形），这样就会保证了当前点是卡住人的视野最大的。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<stack>
#include<cmath>
using namespace std;

const int N = 3e5 + 50;
const double pi = acos(-1.0);
typedef long long ll;
struct p{
    ll x, h, id;
}a[N], ansl[N], ansr[N];
bool cmp(p a, p b){
    return a.x < b.x;
}
stack<int> s;

int main(){
    int T;
    scanf("%d", &T);
    for(int kase = 1; kase <= T; kase++){
        memset(a, 0, sizeof a);
        memset(ansl, 0, sizeof ansl);
        memset(ansr, 0, sizeof ansr);
        int n, m;
        scanf("%d", &n);
        for(int i = 1; i <= n; i++) 
            scanf("%lld %lld", &a[i].x, &a[i].h);
        scanf("%d", &m);
        for(int i = n + 1; i <= n + m; i++) {
            scanf("%lld", &a[i].x);
            a[i].h = -1, a[i].id = i - n;
        }
        int tot = n + m;
        sort(a + 1, a + tot + 1, cmp);
        while(!s.empty()) 
            s.pop();
        for(int i = 1; i <= tot; i++){
            if(a[i].h != -1){
                while(!s.empty() && a[i].h >= a[s.top()].h) s.pop();
                while(s.size() > 1){
                    ll x1 = a[i].x - a[s.top()].x, h1 = a[s.top()].h - a[i].h;
                    int t = s.top();
                    s.pop();
                    ll x2 = a[t].x - a[s.top()].x, h2 = a[s.top()].h - a[t].h;
                    if(h1 * x2 > h2 * x1) {
                        s.push(t);
                        break;
                    }
                }
                s.push(i);
            }else {
                while(s.size() > 1){
                    ll x1 = a[i].x - a[s.top()].x, h1 = a[s.top()].h;
                    int t = s.top();
                    s.pop();
                    ll x2 = a[i].x - a[s.top()].x, h2 = a[s.top()].h;
                    if(h1 * x2 > x1 * h2) {
                        s.push(t);
                        break;
                    }
                }    
                ansl[a[i].id].h = s.top();
                ansl[a[i].id].x = a[i].x;    
            }
        }
        while(!s.empty()) 
            s.pop();
        for(int i = tot; i >= 1; i--){
            if(a[i].h != -1){
                while(!s.empty() && a[s.top()].h <= a[i].h) s.pop();
                while(s.size() > 1){
                    ll x1 = a[s.top()].x - a[i].x, h1 = a[s.top()].h - a[i].h;
                    int t = s.top();
                    s.pop();
                    ll x2 = a[s.top()].x - a[t].x, h2 = a[s.top()].h - a[t].h;
                    if(h1 * x2 > h2 * x1) {
                        s.push(t);
                        break;
                    }
                }
                s.push(i);
            }else {
                while(s.size() > 1){
                    ll x1 = a[s.top()].x - a[i].x, h1 = a[s.top()].h;
                    int t = s.top();
                    s.pop();
                    ll x2 = a[s.top()].x - a[i].x, h2 = a[s.top()].h;
                    if(h1 * x2 > x1 * h2) {
                        s.push(t);
                        break;
                    }
                }    
                ansr[a[i].id].h = s.top();
                ansr[a[i].id].x = a[i].x;    
            }
        }
        printf("Case #%d:\n", kase);
        for(int i = 1; i <= m; i++){
            double j1 = 180 * 1.0 / pi * atan2(1.0 * a[ansl[i].h].h, (ansl[i].x - a[ansl[i].h].x));
            double j2 = 180 * 1.0 / pi * atan2(1.0 * a[ansr[i].h].h, (a[ansr[i].h].x - ansr[i].x));
            printf("%0.10lf\n", 180 - j1 - j2);
        }
    }
    return 0;
}