LDUOJ——山(计算几何+二分+精度)

简介: LDUOJ——山(计算几何+二分+精度)

H. 山

Description

给一座山,如图所示


现在要在山上的某个部位装一盏灯,使得这座山的任何一个部位都能够被看到。给出最小的 y 坐标,如图的+号处 就是 y 坐标最小的安装灯的地方。


Input

第一行一个数 N,表示这座山有 N 个点构成,接下来 N 行从左到右给出了这座山的构造情况,每行两个数 Xi,Yi, 表示一个折点,保证 Xi>Xi−1(1<i≤N)。


Output

仅输出一行,为最小的 y 坐标,当你的答案与标准答案相差 0.01 时,则被认为是正确的。


Samples

Input 复制

6

0 0

10 0

11 1

15 1

16 0

25 0

Output

3.00

Hint

30%的数据,1≤N≤50;

100%的数据,1≤N≤5000,0≤Xi,Yi≤100000,保证答案不超过1000000。

思路:

说说我被这个题卡的全过程。

首先看到最低高度y,第一反应是二分,感觉有单调性,就打算先莽一发。

那么如何check呢。

我们要找到每条直线与当前枚举高度的交点,比如当直线的斜率小于0时,该点左边的所有点都可以看见。又因为要满足所有直线的要求,对所有斜率小于0的直线取左端点max。

同理,当该直线的斜率大于0时,该点右边的点都可以看见,取右端点的min。

当区间合法时,则该高度合法。

但是这样会漏点斜率为0的情况,特判一下就好了。具体原因画画图就出来了。

几个小tips:

1.注意斜率是y/x

2.注意计算过程中的精度和输出答案的精度是不一样的

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll,ll>PLL;
typedef pair<int,int>PII;
typedef pair<double,double>PDD;
#define I_int ll
inline ll read()
{
    ll x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
#define read read()
#define closeSync ios::sync_with_stdio(0);cin.tie(0);cout.tie(0)
#define multiCase int T;cin>>T;for(int t=1;t<=T;t++)
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repp(i,a,b) for(int i=(a);i<(b);i++)
#define per(i,a,b) for(int i=(a);i>=(b);i--)
#define perr(i,a,b) for(int i=(a);i>(b);i--)
ll ksm(ll a,ll b,ll p)
{
    ll res=1;
    while(b)
    {
        if(b&1)res=res*a%p;
        a=a*a%p;
        b>>=1;
    }
    return res;
}
#define PI acos(-1)
const double eps=1e-8;
const int maxn=1e6+100;
struct node{
    double x,y;
}a[maxn];///存点
struct node1{
    double k,b,x;
}b[maxn];///存线段
int n;
bool check(double y){
    double l=-1e9,r=1e9;
    for(int i=2;i<=n;i++){
        if(b[i].k<0){
            double t=(y-b[i].b)/b[i].k;
            l=max(l,t);
        }
        else if(b[i].k>0){
            double t=(y-b[i].b)/b[i].k;
            r=min(r,t);
        }
        else if(b[i].k==0&&y<b[i].b) return 0;
    }
    if(l<=r) return 1;
    return 0;
}
int main()
{
    n=read;
    rep(i,1,n){
        cin>>a[i].x>>a[i].y;
        if(i>1){
            b[i].k=(a[i].y-a[i-1].y)/(a[i].x-a[i-1].x);
            b[i].b=a[i].y-b[i].k*a[i].x;
            b[i].x=a[i].x;
        }
    }
    double l=0,r=1e9;
    double res=0;
    while((r-l)>eps){
        double mid=(l+r)/2;
        if(check(mid)) res=mid,r=mid-eps;
        else l=mid+eps;
        ///cout<<mid<<" "<<l<<" "<<r<<endl;
    }
    printf("%.2lf\n",res);
    return 0;
}
目录
相关文章
|
6月前
|
算法 安全 Java
非启发式算法——中国剩余定理
非启发式算法——中国剩余定理
119 0
|
人工智能 算法
算法提高:组合数学| 容斥原理常见应用
容斥原理常见的问题如下。 (1) 篮球、羽毛球、网球三种运动,至少会一种的有22人,会篮球的有15人,会羽毛球的有17人,会网球的有12人,既会篮球又会羽毛球的有11人,既会羽毛球又会网球的有7人,既会篮球又会网球的有9人,那么三种运动都会的有多少人? (2) 《西游记》《三国演义》《红楼梦》三大名著,至少读过其中一本的有20人,读过《西游记》的有10人,读过《三国演义》的有12人,读过《红楼梦》的有15人,读过《西游记》《三国演义》的有8人,读过《三国演义》《红楼梦》的有9人,读过《西游记》《红楼梦》的有7人。问三本书全都读过的有多少人?
152 0
算法提高:组合数学| 容斥原理常见应用
|
算法
算法提高:计算几何基础 | 详解凸包问题
点集Q的凸包(convex hull)是指一个最小凸多边形,满足Q中的点或者在多边形边上,或者在其内
177 0
算法提高:计算几何基础 | 详解凸包问题
|
6月前
|
Shell
【高数定积分求解旋转体体积】 —— (上)高等数学|定积分|柱壳法|学习技巧
【高数定积分求解旋转体体积】 —— (上)高等数学|定积分|柱壳法|学习技巧
113 0
|
11月前
差分方程模型:兔子繁殖问题(斐波拉契数列)
差分方程模型:兔子繁殖问题(斐波拉契数列)
132 0
|
机器学习/深度学习 算法
算法提高:组合数学| 卡特兰数的实现
卡特兰数列是组合数学中在各种计数问题中常出现的数列,其前几项为1,1,2,5,14,42,132,429,1430,4862,16796,58786,208012…… 卡特兰数首先是由欧拉在计算对凸n边形的不同的对角三角形剖分的个数问题时得到的,即在一个凸n边形中,通过不相交于n边形内部的对角线,把n边形拆分成若干三角形,不同的拆分数用Hn表示,Hn即卡特兰数。
145 0
算法提高:组合数学| 卡特兰数的实现
|
算法 Java C++
蓝桥杯 算法训练 小生物的逃逸(球坐标公式+暴力求解)
蓝桥杯 算法训练 小生物的逃逸(球坐标公式+暴力求解)
|
机器学习/深度学习 算法
模拟退火-n皇后问题
模拟退火-n皇后问题
|
存储 算法
数学知识:高斯消元(一)
复习acwing算法基础课的内容,本篇为讲解数学知识:高斯消元,关于时间复杂度:目前博主不太会计算,先鸽了,日后一定补上。
163 0
数学知识:高斯消元(一)
数学知识:高斯消元(二)
AcWing 884. 高斯消元解异或线性方程组
115 0
数学知识:高斯消元(二)