演算法 三分法

給定若干個點，求包含所有點且與x軸相切的圓的最小半徑

#include<bits/stdc++.h>
using namespace std;
const int MAX=1e5+5;
const double eps=1e-8;
int n;
pair<double,double>a[MAX];
double l,r,midl,midr,lr,rr;
double dis(double x1,double x2,double y2)
{
    return (x1-x2)*(x1-x2)+y2*y2;
}
double getr(double x)
{
    double R=0.0;
    for(int i=0;i<n;++i)
    {
        double d=dis(x,a[i].first,a[i].second);
        R=max(R,fabs(d/a[i].second/2.0));
    }
    return R;
}
int main()
{
    while(~scanf("%d",&n))
    {
        bool yl=0,yr=0;int yo=0;
        l=rr=1e7,r=lr=-1e7;
        for(int i=0;i<n;++i)
        {
            scanf("%lf%lf",&a[i].first,&a[i].second);
            l=min(l,a[i].first);
            r=max(r,a[i].first);
            if(a[i].second>0) yl=1;
            else if(a[i].second<0) yr=1;
            else ++yo;
        }
        if((yl&&yr)||yo>1) {printf("-1\n");continue;}
        if(fabs(l-r)<eps)
        {
            double R=0.0;
            for(int i=0;i<n;++i) R=max(r,fabs(a[i].second));
            printf("%.8f\n",R/2.0);
            continue;
        }
        for(int i=0;i<100;++i)//100次逼近1e-6
        {
            midl=l+(r-l)/3.0;midr=r-(r-l)/3.0;//三分
            lr=getr(midl),rr=getr(midr);
            if(lr-rr>eps) l=midl;
            else r=midr;
        }
        printf("%.8f\n",(lr+rr)/2.0);
    }
    return 0;
}