CFgym 100803D -Space Golf

CFgym 100803D -Space Golf

题意：一子弹从水平地面上原点发出，每次的水平、竖直速度都相同即每次做相同的抛物线运动最后到达终点，中间有n个不计宽度的障碍物，给定终点的位置和每个障碍物的位置和其高度，求整个过程中的可能最小初始速度，子弹落地的次数不能超过b次。

分析：对于起点终点固定的抛物线运动，可以证明当速度为45度时会达到最小，如果不能通过的话，那么也可以证明当运动曲线刚好经过障碍物最高点时速度达到最小。

落地点不超过15次，暴力即可，对于每一个落地次数，知道了每段抛物线运动的水平距离，首先判断45度时能不能通过所有的障碍物，不能的话，有n种可能的速度，在求得刚好经过每一个障碍物最高点时的速度之后判断是不是能够通过所有的障碍物，然后取最小值即可。

整个过程主要是物理公式的推导，这个应该不难，还有就是浮点数比较的精度问题，以前都是懒得eps，直接大于小于等于，这次果断WA了，浮点数的比较还是使用eps避免精度误差！！！

#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#define INF 0x3f3f3f3f
#define EPS 1e-7
using namespace std;

int n, b;
double d, dis, vx, vy;
struct obstacles{
    double h, x;
}ob[20];
double x[20];

void cal(double x, double h)
{
    double a = h / ((x - dis) * x);
    double b = -a * dis;
    double k = dis / 2;
    vy = sqrt(2*(a*k*k + b*k));
    vx = k / vy;
}

int judge()
{
    for (int i = 0; i < n; i++){
        double h = vy*x[i] / vx - 0.5*x[i]*x[i] / (vx*vx);
        if (ob[i].h-h >= EPS) return 0;
    }
    return 1;
}

double solve(int cnt)
{
    dis = d/cnt;
    for (int i = 0; i < n; i++) x[i] = fmod(ob[i].x, dis);
    vy = vx = sqrt(dis/2);
    if (judge()) return vx * sqrt(2.0);
    double ma = INF;
    for (int i = 0; i < n; i++){
        cal(x[i], ob[i].h);
        if (judge()) ma = min(ma, sqrt(vx*vx + vy*vy));
    }
    return ma;
}

int main()
{
    scanf("%lf %d %d", &d, &n, &b);
    for (int i = 0; i < n; i++) scanf ("%lf %lf", &ob[i].x, &ob[i].h);
    double ans = INF;
    for (int i = 1; i <= b+1; i++) ans = min(ans, solve(i));
    printf ("%.8f\n", ans);
    return 0;
}