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<cstdio>
#include<cmath>
#define eps 1e-4
using namespace std;
double a,b,c,x,y;
double l,r,p,mid1,mid2,k1,k2,ans1,ans2;
double f(double fx)
{
     double fy=a*fx*fx+b*fx+c;
     return sqrt((x-fx)*(x-fx)+(y-fy)*(y-fy));
}
int main()
{
    scanf("%lf%lf%lf%lf%lf",&a,&b,&c,&x,&y);
    l=-200;r=200;
    while(r-l>eps)
    {
        //mid1=(2*l+r)/3;mid2=(l+2*r)/3;
//这种三分姿势不能过
       //改成这样可以    mid1=l+(r-l)/3;mid2=r-(r-l)/3; 
　　　　mid1=(l+r)/2;mid2=(mid1+r)/2;
        k1=f(mid1);k2=f(mid2);
        if(k1<=k2) r=mid2;
        else l=mid1;
    }
    printf("%.3lf",f(l));
    return 0;
}

问：为什么可以三分横坐标

因为横坐标卡到4位小数，算出的答案不止4位了