牛牛战队的比赛地(三分求极值)

题目

链接:https://ac.nowcoder.com/acm/contest/3006/B
来源:牛客网

由于牛牛战队经常要外出比赛,因此在全国各地建立了很多训练基地,每一个基地都有一个坐标(x,y)。
这周末,牛牛队又要出去比赛了,各个比赛的赛点都在x轴上。牛牛战队为了方便比赛,想找一个到达训练基地最大距离最小的地方作为比赛地。
这个问题对于牛牛战队太简单了,它就交给了你,你来帮他算一下~

输入描述

输入数据第一行包含一个整数N(1≤N≤100000),表示牛牛战队训练基地的数量。

接下来N行,每行包括2个整数x,y(−10000≤x,y≤10000),表示每一个训练基地的坐标。

输出描述

输出一个小数,表示选择的比赛地距离各训练基地最大距离的最小值。

如果你的答案是a,标准答案是b,当|a-b|/max(1,|b|)≤0.0001) 时,你的答案将被判定为正确。

思路

因为dis 2 ^2 2=x 2 ^2 2+y 2 ^2 2,y值不变 所以所求值只和x值有关。这是一个凹函数即有极小值。我们三分答案,l,r 当用l求的的答案大于等于r的时候,我们就可知答案在l和ri之间,反之则在le和r之间。

代码

#include<bits/stdc++.h>
using namespace std;
#define x first
#define y second
typedef long long ll;
typedef pair<int,int> pii;
typedef unsigned long long ull;
const int N=1e5+7;
const int inf=0x3f3f3f3f;
const double eps=1e-4;
int n;
struct node
{
	double x,y;
}point[N];

double dis(node a, node b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

double getmax(double x)
{
	double maxv=eps;
	node q={x,0};
	for(int i=0;i<n;i++)
	{
		double p=dis(q,point[i]);
		maxv=max(maxv,p);
	}
	return maxv;
}

int main()
{
	//freopen("test.in","r",stdin);//设置 cin scanf 这些输入流都从 test.in中读取
    //freopen("test.out","w",stdout);//设置 cout printf 这些输出流都输出到 test.out里面去
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	cin>>n;
	double le=inf,ri=-inf;
	for(int i=0;i<n;i++)
	{
		cin>>point[i].x>>point[i].y;
		le=min(le,point[i].x);
		ri=max(ri,point[i].x);
	}
	while((ri-le)>eps)
	{
		double mid=(ri-le)/3;
		double l,r;
		l=le+mid,r=le+2*mid;
		if(getmax(l)>=getmax(r)) le=l;
		else ri=r;
	}
	cout<<getmax(le)<<endl;
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值