hihocoder 1142 三分·三分求极值

#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 <bits/stdc++.h>
using namespace std;
#define eqs 0.000001
int a,b,c,x,y;
double pansum(double x1)
{
    double y1=a*x1*x1+b*x1+c;
    return sqrt((x-x1)*(x-x1)+(y-y1)*(y-y1));
}
double pan(double l,double r)
{
    while(r-l>eqs)
    {
        double ml=l+(r-l)/3.0;
        double mr=r-(r-l)/3.0;
        if(pansum(ml)<pansum(mr))
        {
            r=mr;
        }
        else
        {
            l=ml;
        }
    }
    return pansum(l);
}
int main()
{
    while(~scanf("%d%d%d%d%d",&a,&b,&c,&x,&y))
    {
        double k=-b*1.0/(2*a);
        if(k<x)
        {
            printf("%.3f\n",pan(k,205));
        }
        else
        {
            printf("%.3f\n",pan(-205,k));
        }
    }
}