zoj 2340 Little Jumper（三分 物理好题）

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1395
大意：青蛙能从一个点跳到第三个点，如图，需要跳两次。问整个过程的最大起跳速度。

分析：对于每一次跳跃，设置水平的运动距离总长度是L，到达ti，bi时已经运动的水平距离是S，那么由物体的竖直抛体运动可以得到相应的等式。

$$L=v_xt \\ v_y=\frac{gt}{2} \\ t=\frac{L}{v_x}=\frac{2v_y}{g} \\ \longrightarrow \\ v_xv_y=\frac{gL}{2} \\ v=\sqrt{v_x^2+v_y^2}\geq \sqrt{2v_xv_y}=\sqrt{gL}$$
即当 $v_x=v_y$时，v最小。速度夹角成45度。
（1）这是理想情况：
$$h=v_yt-\frac{1}{2}gt^2 =v_y\frac{S}{v_x}-\frac{1}{2}g(\frac{S}{v_x})^2 =S-\frac{S^2}{l} \\ B \leq h\leq T$$
当非理想情况时，那么计算出来的h不在相应的区间内。想要成功的跳跃，只有将速度改变（大小，方向）。
（2） $h>T:$
设速度夹角delta:

想要不撞上墙，需要将速度夹角变小。尽量接近理想状态。
$$h=v_yt-\frac{1}{2}gt^2=T \\v_yt=T+\frac{1}{2}gt^2 =T+\frac{1}{2}g(\frac{S}{v_x})^2 \\v_y\frac{S}{v_x}=T+\frac{1}{2}g\frac{S^2}{v_x^2} \\v_xv_yS=Tv_x^2+\frac{1}{2}gS^2 \\ \frac{gLS}{2}=Tv_x^2+\frac{1}{2}gS^2 \\ v_x^2=\frac{gS(L-S)}{2T} \\ v_y^2=\frac{1}{t^2}(T+\frac{1}{2}g\frac{S^2}{v_x^2})^2 =\frac{v_x^2}{S^2}(T+\frac{1}{2}g\frac{S^2}{v_x^2})^2 \\ v=\sqrt{v_x^2+v_y^2}$$
（3）当h < B时，我们相同过程推出:
$$v_x^2=\frac{gS(L-S)}{2B} \\ v_y^2=\frac{v_x^2}{S^2}(B+\frac{1}{2}g\frac{S^2}{v_x^2})^2 \\ v=\sqrt{v_x^2+v_y^2}$$

solution:
上诉过程的基本假设是已知水平运动的距离，由此我们直接枚举那个距离，因为存在最小的速度，所以有技巧的枚举：三分。
注意此题需要人为的设定迭代次数，不然会超时。

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const double eps=1e-15;
double b1,t1,b2,t2,l,ds,df,g;
double calc(double L,double S,double T,double B){
     double h=S-S*S/L;
     if(h-T>eps){
        double xx=g*S*(L-S)/2/T;
        double yy=xx/S/S*(T+0.5*g*S*S/xx)*(T+0.5*g*S*S/xx);
        return sqrt(xx+yy);
     }
     if(h-B<-eps){
        double xx=g*S*(L-S)/2/B;
        double yy=xx/S/S*(B+0.5*g*S*S/xx)*(B+0.5*g*S*S/xx);
        return sqrt(xx+yy);
     }
     return sqrt(g*L);
}
double gv(double x){
     double v1=calc(x+ds,ds,t1,b1),v2=calc(l-x+df,df,t2,b2);  //L:总的长度 S:越点  T  B
     return v1-v2>eps?v1:v2;
}
int main()
{
    //freopen("cin.txt","r",stdin);
    while(~scanf("%lf%lf",&b1,&t1)){
        scanf("%lf %lf %lf %lf %lf %lf",&b2,&t2,&l,&ds,&df,&g);
        double lt=0,rt=l,mid1,mid2;
        double y1,y2;
        int c=0;
        while(rt-lt>eps && c++<100){
            mid1=lt+(rt-lt)/3;
            mid2=rt-(rt-lt)/3;
            y1=gv(mid1);  y2=gv(mid2);
            if(y1<y2) rt=mid2;
            else lt=mid1;
        }
        printf("%.6lf\n",gv(lt));
    }
    return 0;
}