[HiHoCoder]三分·三分求极值[hiho一下 第四十周]

时间限制:10000ms

单点时限:1000ms

内存限制:256MB

描述

在直角坐标系中有一条抛物线y=ax^2+bx+c和一个点P(x,y),求点P到抛物线的最短距离d。

提示:三分法

输入

第1行:5个整数a,b,c,x,y。前三个数构成抛物线的参数,后两个数x,y表示P点坐标。-200≤a,b,c,x,y≤200

输出

第1行:1个实数d,保留3位小数(四舍五入)

样例输入
2 8 2 -2 6
样例输出
2.437
import java.math.BigDecimal;
import java.util.Scanner;

public class Main {

	private static String minLen(int a,int b,int c,int x,int y){
		double zhou = -((double)b)/(2*a);
		double left = -1;
		double right = -1;
		if(b*b-4*a*(c-y)<0){
			left = zhou<x?zhou:x;
			right = zhou<x?x:zhou;
		}else{
			double as = x<zhou?(-b-Math.sqrt(b*b-4*a*(c-y)))/(2*a):(-b+Math.sqrt(b*b-4*a*(c-y)))/(2*a);
			left = as<x?as:x;
			right = as<x?x:as;
		}
	/*	left = x>zhou?zhou:left;
		right = x<=zhou?zhou:right;*/
		BigDecimal res = null;
		while(true){
			double lm = left+(right-left)/3;
			double rm = left+(right-left)/3*2;
			double lmval = len(a,b,c,x,y,lm);
			double rmval = len(a,b,c,x,y,rm);
			BigDecimal lvStr = new BigDecimal(lmval).setScale(3,BigDecimal.ROUND_HALF_UP);
			BigDecimal rvStr = new BigDecimal(rmval).setScale(3,BigDecimal.ROUND_HALF_UP);
			if(lvStr.equals(rvStr)&&right-left<0.0001){
				res = lvStr;
				break;
			}
			if(lmval<rmval){
				right = rm;
			}else{
				left = lm;
			}
		}
		return res.setScale(3,BigDecimal.ROUND_HALF_UP)+"";
	}

	private static double len(int a,int b,int c,int x,int y,double xx){
		return Math.sqrt((xx-x)*(xx-x)+Math.pow(a*xx*xx+b*xx+c-y, 2));
	}

	public static void main(String args[]){
		Scanner sc = new Scanner(System.in);
		int a = sc.nextInt();
		int b = sc.nextInt();
		int c = sc.nextInt();
		int x = sc.nextInt();
		int y = sc.nextInt();
		System.out.println(minLen(a, b, c, x, y));
		sc.close();
	}
}
时间: 2024-10-12 08:19:31

[HiHoCoder]三分·三分求极值[hiho一下 第四十周]的相关文章

hiho一下 第四十周 题目1 : 三分&#183;三分求极值

时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 这一次我们就简单一点了,题目在此: 在直角坐标系中有一条抛物线y=ax^2+bx+c和一个点P(x,y),求点P到抛物线的最短距离d. 提示:三分法 × 提示:三分法 在之前的几周中我们了解到二分法作为分治中最常见的方法,适用于单调函数,逼近求解某点的值. 但当函数是凸形函数时,二分法就无法适用,这时就需要用到三分法. 从三分法的名字中我们可以猜到,三分法是对于需要逼近的区间做三等分: 我们发现lm这个点比rm要低,那

HihoCoder 1142-三分求极值(三分模板)

描述 这一次我们就简单一点了,题目在此: 在直角坐标系中有一条抛物线y=ax^2+bx+c和一个点P(x,y),求点P到抛物线的最短距离d. 输入 第1行:5个整数a,b,c,x,y.前三个数构成抛物线的参数,后两个数x,y表示P点坐标.-200≤a,b,c,x,y≤200 输出 第1行:1个实数d,保留3位小数(四舍五入) 样例输入 2 8 2 -2 6 样例输出 2.437 思路:三分裸题.注意一下变量类型. 1 #include<cstdio> 2 #include<cmath&g

hihoCoder hiho一下 第四十八周 题目1 : 拓扑排序&#183;二

题意:给定一个拓扑图,其中部分结点含有1个病毒,每个结点只要收到病毒就会立即往出边所能到达的点传播,病毒数可叠加,求所有结点的病毒数总和. 思路:根据拓扑的特点,每个入度为0的点肯定不会再被传播病毒,而且会将自己的所有病毒向与其相连的结点传播.那么可以从入度味为0的点着手,逐个删除入度为0的结点,在删除的过程中,更新与其相连的结点的病毒数(即将病毒数累加到该结点),到最后所有结点都没了,各个结点所累积的病毒数的和就是答案. 1 #include<bits/stdc++.h> 2 using n

hihoCoder hiho一下 第四十六周 博弈游戏&#183;Nim游戏&#183;三( sg函数 )

题意:给出几堆石子数量,每次可以取走一堆中任意数量的石头,也可以将一堆分成两堆,而不取.最后取走者胜. 思路:先规矩地计算出sg值,再对每个数量查SG值就可以了.最后求异或和. 1 #include <bits/stdc++.h> 2 using namespace std; 3 const int N=105, limit=20004; 4 int a[N],n,sg[limit]={0,1,2}; 5 bool B[limit]; 6 int main() 7 { 8 //freopen(

hiho一下 第四十二周

题目1 : 骨牌覆盖问题·二 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 上一周我们研究了2xN的骨牌问题,这一周我们不妨加大一下难度,研究一下3xN的骨牌问题? 所以我们的题目是:对于3xN的棋盘,使用1x2的骨牌去覆盖一共有多少种不同的覆盖方法呢? 首先我们可以肯定,奇数长度一定是没有办法覆盖的:对于偶数长度,比如2,4,我们有下面几种覆盖方式: 提示:3xN骨牌覆盖 输入 第1行:1个整数N.表示棋盘长度.1≤N≤100,000,000 输出 第1行:1个

hiho一下 第五十周 题目1 : 欧拉路&#183;二

时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 在上一回中小Hi和小Ho控制着主角收集了分散在各个木桥上的道具,这些道具其实是一块一块骨牌. 主角继续往前走,面前出现了一座石桥,石桥的尽头有一道火焰墙,似乎无法通过. 小Hi注意到在桥头有一张小纸片,于是控制主角捡起了这张纸片,只见上面写着: 将M块骨牌首尾相连放置于石桥的凹糟中,即可关闭火焰墙.切记骨牌需要数字相同才能连接. ——By 无名的冒险者 小Hi和小Ho打开了主角的道具栏,发现主角恰好拥有M快骨牌. 小H

hiho一下 第四十九周(欧拉路的判定)49

/////////////////////////////////////////////////////////////////////////////////////////////////////// 作者:tt267 声明:本文遵循以下协议自由转载-非商用-非衍生-保持署名|Creative Commons BY-NC-ND 3.0 转载请注明:http://blog.csdn.net/tt2767/article/details/45420067 ///////////////////

hiho一下 第五十周 欧拉路&#183;二

题目链接:click here~~ [题目大意] 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 在上一回中小Hi和小Ho控制着主角收集了分散在各个木桥上的道具,这些道具其实是一块一块骨牌. 主角继续往前走,面前出现了一座石桥,石桥的尽头有一道火焰墙,似乎无法通过. 小Hi注意到在桥头有一张小纸片,于是控制主角捡起了这张纸片,只见上面写着: 将M块骨牌首尾相连放置于石桥的凹糟中,即可关闭火焰墙.切记骨牌需要数字相同才能连接. --By 无名的冒险者 小Hi和小Ho打

hiho一下 第五十周(欧拉路径)50

/////////////////////////////////////////////////////////////////////////////////////////////////////// 作者:tt267 声明:本文遵循以下协议自由转载-非商用-非衍生-保持署名|Creative Commons BY-NC-ND 3.0 查看本文更新与讨论请点击:http://blog.csdn.net/tt2767/article/details/46493073 ////////////