ZOJ 1104 Leaps Tall Buildings

It's a bird! It's a plane! It's coming right at us!

Although it sometimes seems like it, Superman can't fly (without a plane). Instead, he makes super-human leaps, especially over tall buildings. Since he never knows when he will need to catch a criminal, he can't register flight paths. To avoid hitting planes, he tries to keep his jumps as low to the ground as he can. Given a city-scape as input, find the angle and velocity of Superman's jump that minimizes his maximum altitude.

Recall that gravity provides an acceleration of 9.8 m/s² downwards and the formula for Superman's vertical distance from his starting location is d(t)=v t + 0.5 a t²where v is his initial velocity, a is his acceleration and t is time in seconds since the start of the leap.

Input:

Input consists of a sequence of city-scapes, each of the form

n
0 d₁
h₂ d₂
:
h_(n-1) d_(n-1)
0 d_n

Superman starts at ground level and leaps d₁+...+d_n metres, landing at ground level and clearing all of the buildings at heights h₂ to h_(n-1), each with the given widths. n will be at most 100.

Output:

Output is the angle and initial velocity that minimizes the height that Superman attains, both appearing on the same line. The values should be given to two decimal places and be accurate within 0.01 degrees or m/s, as appropriate.

Sample Input:

3
0 5
10 5
0 5
5
0 10.5
20 11.5
25 10
10 15
0 7

Diagram for Second City-scape

 
(Not to scale.)

Sample Output:

71.57 15.65
67.07 27.16

二分最高点，判断即可。
#include<map>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<math.h>
using namespace std;
const int N=1e5+10;
const double eps = 1e-8;
double x[N],y[N];
int n,sz;
double v,c;

bool check(double h)
{
    double v1=sqrt(2*9.8*h);
    double v2=x[sz-1]/(v1/4.9);

    v=sqrt(v1*v1+v2*v2);
    c=180*acos(v2/v)/acos(-1.0);

    for (int i = 0;i<sz;i++)
    {
        double xx = x[i] > (x[sz-1]/2)?x[sz-1]-x[i]:x[i];
        double t = xx / v2;
        double hh = v1*t-4.9*t*t;
        if (hh<y[i]) return 0;
    }
    return 1;
}

int main()
{

    while(~scanf("%d",&n))
    {
        double xx = sz = 0,a,b;
        for (int i=1;i<=n;i++)
        {
            scanf("%lf%lf",&a,&b);
            x[sz] = xx; y[sz++] = a;
            //printf("%lf %lf\n",x[sz-1],y[sz-1]);
            x[sz] = xx+=b; y[sz++] =a;
            //printf("%lf %lf\n",x[sz-1],y[sz-1]);
        }
        double l=0,r=1LL<<60;
        while (l+eps<r)
        {
            double mid=(l+r)/2;
            if (check(mid)) r=mid;
            else l=mid;
        }
        //printf("%lf\n",r);
        check(r);
        printf("%lf %lf\n",c,v);
    }
    return 0;
}