hihoCoder#1142(三分求极值)

时间限制: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

#include <cstdio>
#include <cmath>
using namespace std;
const double EPS=1.0e-6;
int a,b,c,x,y;
double dist(double x0)
{
    return sqrt((x0-x)*(x0-x)+(a*x0*x0+b*x0+c-y)*(a*x0*x0+b*x0+c-y));
}
int main()
{
    scanf("%d%d%d%d%d",&a,&b,&c,&x,&y);
    if(a==0&&b==0)
    {
        double res=fabs(y-c);
        printf("%.3f\n",res);
    }
    else
    {
        double l=-200.0;
        double r=200.0;
        while(fabs(r-l)>EPS)
        {
            double div3=(r-l)/3;
            double lm=l+div3;
            double rm=r-div3;
            if(dist(lm)>dist(rm))
            {
                l=lm;
            }
            else
            {
                r=rm;
            }
        }
        printf("%.3f\n",dist(r));
    }
    return 0;
}