bzoj千题计划126:bzoj1038: [ZJOI2008]瞭望塔
Posted 日拱一卒 功不唐捐
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj千题计划126:bzoj1038: [ZJOI2008]瞭望塔相关的知识,希望对你有一定的参考价值。
http://www.lydsy.com/JudgeOnline/problem.php?id=1038
本题可以使用三分法
将点按横坐标排好序后
对于任意相邻两个点连成的线段,瞭望塔的高度 是单峰函数,而且是下凸函数
感性理解单峰就是
瞭望塔建的靠左,为了能看到右边的,要高一点
瞭望塔建的靠右,为了能看到左边的,要高一点
所以
枚举所有线段,三分线段上建造瞭望塔的位置,所有线段上的瞭望塔高度取最小
#include<cmath> #include<cstdio> #include<algorithm> using namespace std; #define N 301 const double eps=1e-7; int n; struct POINT { double x,y; }e[N]; struct LINE { double k,b; bool exit; }l[N]; bool cmp(POINT p,POINT q) { return p.x<q.x; } double cal(double x,double y) { double tmp=0; for(int i=1;i<n;++i) { if(!l[i].exit) continue; tmp=fmax(tmp,l[i].k*x+l[i].b-y); } return tmp; } int main() { scanf("%d",&n); for(int i=1;i<=n;++i) scanf("%lf",&e[i].x); for(int i=1;i<=n;++i) scanf("%lf",&e[i].y); sort(e+1,e+n+1,cmp); for(int i=1;i<n;++i) { if(fabs(e[i].x-e[i+1].x)<eps) continue; l[i].k=(e[i].y-e[i+1].y)/(e[i].x-e[i+1].x); l[i].b=e[i].y-e[i].x*l[i].k; l[i].exit=true; } double L,R,mid1,mid2; int T; double ans=1e20; for(int i=1;i<n;++i) { if(!l[i].exit) continue; T=100; L=e[i].x; R=e[i+1].x; while(T--) { mid1=L+(R-L)/3; mid2=R-(R-L)/3; if(cal(mid1,l[i].k*mid1+l[i].b)<cal(mid2,l[i].k*mid2+l[i].b)) R=mid2; else L=mid1; } ans=fmin(ans,cal(mid1,l[i].k*mid1+l[i].b)); } printf("%.3lf",ans); }
1038: [ZJOI2008]瞭望塔
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2658 Solved: 1084
[Submit][Status][Discuss]
Description
致力于建设全国示范和谐小村庄的H村村长dadzhi,决定在村中建立一个瞭望塔,以此加强村中的治安。我们
将H村抽象为一维的轮廓。如下图所示 我们可以用一条山的上方轮廓折线(x1, y1), (x2, y2), …. (xn, yn)来描
述H村的形状,这里x1 < x2 < …< xn。瞭望塔可以建造在[x1, xn]间的任意位置, 但必须满足从瞭望塔的顶端可
以看到H村的任意位置。可见在不同的位置建造瞭望塔,所需要建造的高度是不同的。为了节省开支,dadzhi村长
希望建造的塔高度尽可能小。请你写一个程序,帮助dadzhi村长计算塔的最小高度。
Input
第一行包含一个整数n,表示轮廓折线的节点数目。接下来第一行n个整数, 为x1 ~ xn. 第三行n个整数,为y1
~ yn。
Output
仅包含一个实数,为塔的最小高度,精确到小数点后三位。
Sample Input
6
1 2 4 5 6 7
1 2 2 4 2 1
【输入样例二】
4
10 20 49 59
0 10 10 0
Sample Output
1.000
【输出样例二】
14.500
HINT
N ≤ 300,输入坐标绝对值不超过106,注意考虑实数误差带来的问题。
以上是关于bzoj千题计划126:bzoj1038: [ZJOI2008]瞭望塔的主要内容,如果未能解决你的问题,请参考以下文章
bzoj千题计划118:bzoj1028: [JSOI2007]麻将
bzoj千题计划144:bzoj1176: [Balkan2007]Mokia