洛谷 P1663 山

一看就是有意思的题目；
然后去考fop_zz；
他说他做过了；

https://www.luogu.org/problem/show?pid=1663
我们怎么搞呢；
显然就是求很多半平面的交；
我们考虑二分这个答案；
那么怎么去验证呢？
我们对于每一条直线解析式
y=ax+b
如果a<0,那么x的显然要在这条线的右边；
如果a>0，那么x要在其左边；
如果a==0，那么y就要在这条线的上面；
所以我们通过这3个条件去做就好了；
对于二分的每一个y;
找到x可取范围左边x1;
右边x2;
如果x2>=x1就是可以的；

#include<bits/stdc++.h>
using namespace std;
const int N=5e3+5;
int x[N],y[N];
double a1[N],b1[N],a2[N],b2[N],l,r,mid;
int n,m,t1,t2;
void make(double x1,double y1,double x2,double y2){
    if(fabs(y1-y2)<1e-8){l=max(l,y1);return;}
    double A=(y1-y2)/(x1-x2);
    double B=y1-A*x1;
    if(A<0){t1++;a1[t1]=A;b1[t1]=B;}
    else {t2++;a2[t2]=A;b2[t2]=B;}
}
bool check(double y){
    double x1=0,x2=x[n],X;
    for(int i=1;i<=t1;i++)x1=max(x1,(y-b1[i])/a1[i]);
    for(int i=1;i<=t2;i++)x2=min(x2,(y-b2[i])/a2[i]);
    return (x2>x1)||(fabs(x2-x1)<1e-8);
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%d%d",&x[i],&y[i]);
    for(int i=1;i<n;i++)make(x[i],y[i],x[i+1],y[i+1]);
    r=1e6;
    while((r-l)>1e-4){
        mid=(l+r)/2.;
        if(check(mid))r=mid;else l=mid;
    }
    printf("%.2lf",mid);
}