UVALive 6835 - Space Golf（简单解析几何）

题目链接：
UVALive 6835 - Space Golf
题意：
一颗子弹最多可以跳跃b次击中d，有n个木板，在跳跃过程中如果有木板挡住了，就不能往后跳跃了。
求最小初速度。
分析；
因为d已知，可以枚举跳跃次数，那么每次跳跃的长度就知道了。也就是抛物线的开口长度。然后枚举每个木板，得到要跳跃这个木板的最低抛物线高度，然后就知道了需要跳过所有木板的最低高度。
通过物理学公式和数学知识，可以推出抛物线的最低高度和初速度Vx和Vy的关系。整理下可以得到:
V^2=m*m/8/y_max+2*y_max.左侧是最低速度的平方，m是开口宽，y_max是最低抛物线高度。
注意到没有上面是个对勾函数！(双飞燕函数)，所以不是简单的取y_max就是速度最小值，还要和函数极点(m/4)比较！不信的话，试试第二个样例就知道了！

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <climits>
#include <iostream>
using namespace std;
const int MAX_N=10010;
const double eps=1e-6;

double d;
int n,b;
int x[50],y[50];

int main()
{
    while(~scanf("%lf%d%d",&d,&n,&b)){
        for(int i=0;i<n;i++){
            scanf("%d%d",&x[i],&y[i]);
        }
        double ans=1000000000.0;
        for(int i=1;i<=b+1;i++){
            double m=d/i;
            int ok=1;
            double y_max=-1;
            for(int j=0;j<n;j++){
                double tx=x[j],ty=y[j];
                while(tx>=m){
                    tx-=m;
                }
                if(fabs(tx)<=eps){ //判断是不是有木板恰好出于过程中的落地点
                    ok=0;
                    break;
                }
                double tmp_y_max=m*m*ty/(m-tx)/tx/4.0;
                y_max=max(y_max,tmp_y_max);
            }
            if(ok==0) continue; //落地点有木板肯定不行
            if(y_max<m/4+eps) ans=min(m,ans);
            else ans=min(m*m/8/y_max+2*y_max,ans);
        }
        printf("%.6f\n",sqrt(ans));
    }
    return 0;
}