Codeforces#514D(三分,简单二维几何)

#include<bits/stdc++.h>
using namespace std;
const double eps=1e-8;
int n;
struct node{
    double x,y;
}a[1000007];
long double find_r(long double x){
    long double r=0;
    for(int i=1;i<=n;i++){
        r=max(r,(x-a[i].x)*(x-a[i].x)/a[i].y/2+a[i].y/2);//寻找最大的半径
    }
    return r;
}
int main(){
    int flag=0,flag1=0,flag2=0;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%lf%lf",&a[i].x,&a[i].y);
        if(a[i].y>0&&flag1==0)
            flag1=1,flag++;
        if(a[i].y<0&&flag2==0)
            flag2=1,flag++;
        if(a[i].y<0)
            a[i].y-=2*a[i].y;
    }
    int can=0;
    if(flag==2){
        can=1;
    }
    if(can==1){
        puts("-1");
    }
    else{
        long double l=-1e7,r=1e7,lx,rx,dx;
        while(r-l>eps){//寻找最适中的圆心横坐标
            dx=(r-l)/3.0;
            lx=l+dx;
            rx=r-dx;
            if(find_r(lx)-find_r(rx)<0)
                r=rx;
            else
                l=lx;
        }
        long double x=find_r(r);
        printf("%.8lf\n",x);//Lf需用C++14及以上提交,C++11会WA
    }
    return 0;
}

转载于:https://www.cnblogs.com/ldudxy/p/9769132.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值