HDU5033 - Building(暴力，DP思想)

最新推荐文章于 2020-03-17 01:18:21 发布

dawxy

最新推荐文章于 2020-03-17 01:18:21 发布

阅读量589

点赞数

分类专栏：竞赛题解文章标签： DP

本文链接：https://blog.csdn.net/u012760629/article/details/39493371

版权

竞赛题解专栏收录该内容

3 篇文章 0 订阅

订阅专栏

题目链接 HDU5033

【题意】在直线上给出n幢房子的坐标和高度，求出人站在任意没有房子的位置上时能看到的天空的角度范围；房子坐标不会重复，并且保证人站的位置至少左右各有一座房子，注意坐标和高度是实数。

【分析】只要DP思想那样记录每个点左右两边斜率最高的那个点，然后查询的时候二分找查询点在哪两个点之间，再计算一遍左右两边最大的斜率即可求出答案，因为用了这个思想DP查询时间复杂度很低的，就算有极端数据斜率全部递减也不会很慢。比赛的时候竟然想着计算每个点高度为0时的左右最高斜率，一直WA到死，后来才发现一个是要减去当前高度后的斜率。

【AC CODE】546ms

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define toj(a) (a*180.0/PI)
#define MAXN 100010
struct NODE{
    double x, h;
    bool operator<(const NODE& a) const{
        return x < a.x;
    }
}p[MAXN];
struct A{
    int x;
    double val;
}left[MAXN], right[MAXN];
const double PI = 4.0 * atan(1.0);
int n;

inline double clc_k(const NODE& a, const NODE& b){
    return (a.h-b.h)/fabs(a.x-b.x);
}
int find(double num)
{
    int l = 0, r = n-1, mid;
    while(l <= r)
    {
        mid = (l+r)>>1;
        if(p[mid].x < num) 
            l = mid+1;
        else 
            r = mid-1;
    }
    return l;
}
int main()
{
#ifdef SHY
    freopen("e:\\1.txt","r",stdin);
#endif
    int t, count = 0;
    scanf("%d%*c", &t);
    while(t--)
    {
        scanf("%d%*c", &n);
        for(int i = 0; i < n; i++)
            scanf("%lf %lf%*c", &p[i].x, &p[i].h);
        sort(p,p+n);
        left[0].x = -1, left[0].val = 0;
        for(int i = 1; i < n; i++)
        {
            int a = i-1;
            while(~left[a].x && left[a].val > clc_k(p[a], p[i]))
                a = left[a].x;
            left[i].x = a;
            left[i].val = clc_k(p[a], p[i]);
        }
        right[n-1].x = -1, right[n-1].val = 0;
        for(int i = n-2; i >= 0; i--)
        {
            int a = i+1;
            while(~right[a].x && right[a].val > clc_k(p[a], p[i]))
                a = right[a].x;
            right[i].x = a;
            right[i].val = clc_k(p[a], p[i]);
        }
        int q;
        NODE a;
        a.h = 0;
        scanf("%d%*c", &q);
        printf("Case #%d:\n", ++count);
        while(q--)
        {
            scanf("%lf%*c", &a.x);
            int rt = find(a.x);
            int lt = rt-1;
            while(~left[lt].x && left[lt].val > clc_k(p[lt], a))
                lt = left[lt].x;
            while(~right[rt].x && right[rt].val > clc_k(p[rt], a))
                rt = right[rt].x;
            printf("%.10lf\n", 180-toj(atan(clc_k(p[lt],a)))-toj(atan(clc_k(p[rt],a))));
        }
    }
    return 0;
}