不要打错啊,通过抛物线性质检查代码
l=0.0;
l1=calc(x[i],y[i],l);
r=(double)z[i];
r1=calc(x[i],y[i],r);
while (r-l>0.000001)
{
m1=l+(r-l)/3;
m2=r-(r-l)/3;
n1=calc(x[i],y[i],m1);
n2=calc(x[i],y[i],m2);
if (n1<n2)
{
r=m2;
r1=n2;
}else
{
l=m1;
l1=n1;
}
}
ans=min(ans,l1);