hdu2529Shot+三分

Description
“超人”霍华德在NBA扣篮大赛上要求把篮筐升高表演扣篮，但是却没有得到批准，现在我们的hhb也想要给大家表演一个，但是这次他想要表演的是投篮，篮筐放得越高，hhb的表演当然就越精彩，所以现在请你帮助hhb计算出篮筐离地最高能放多高。

假设把球投到篮框的高度就算球进。忽略球和框的大小。

Input
多组测试数据
每组测试数据包括3个浮点数h, l, v (1<=h<=2,1<=l<=100,0 < v<=100) , 分别表示hhb的出手点高度，出手点离篮筐的水平距离和hhb投出的球的速度，h=l=v=0表示输入结束（球在运动过程中只受重力的作用，不受其它任何力的作用，重力加速度取g=9.8，题目中所有物理量均是国际单位制）

Output
对于每组数据，输出一行，包含一个数，在球可以投进篮筐的情况下，篮筐可以离地的最高高度（输入保证这个值一定大于0）。输出保留到2位小数。

Sample Input

1.5 5.0 7.0
0 0 0

Sample Output

1.50

直接就是凸性函数。。
vcosO*t=l
y=vsinO*t-1/2*g*t^2

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pi acos(-1.0)
#define eps 1e-6
double h,l,v;
double geth(double c){
    return l*sin(c)/cos(c)-l*l*9.8/2/v/v/cos(c)/cos(c);
}

int main(){
    while(cin>>h>>l>>v){
        if(h==0&&l==0&&v==0) break;
        double left=0,right=pi/2;
        while(abs(right-left)>eps){
            double m1=(2*left+right)/3;
            double m2=(left+2*right)/3;
            if(geth(m1)>geth(m2)) right=m2;
            else left=m1;
        }
        cout<<setprecision(2)<<std::fixed<<geth(left)+h<<endl;
    }
    return 0;
}