zoj 1104

1. /**
2. * 由于抛物线经过点（0，0）， （D， 0),设抛物线方程为 y = -Ax(x-D), 其中（D， 0）为落地点。
3. * 再由参数方程 x = v0*t*cosO
4.                y = -0.5*g*t^2 + v0*t*sin0
6.   消去参数得   y = -[g/(2*v0^2*cos0^2)]*x^2 + tan0 *x
8.   与前面方程比较得 AD = tan0, v0 = sqrt( g*cos0*cos0/(2*A) ).
10.   所以， 只要求得 A就可以了。
12.   对每个建筑物的两个最高点， 检验通过该点的形如 y = -Ax(x-D)的抛物线是否会碰到其他建筑物， 
13.   在不会碰到其他建筑物的抛物线中， 取height最小的。
15. *
16. *
17. **/
19. #include <iostream>
20. #include <cmath>
21. #include <cstdlib>
22. using namespace std;
24. const double EP = 1E-10;
25. const double INF = 1E200;
26. const double PI = 3.1415926;
28. int main()
29. {
30.     int n;
31.     double d[100], h[100], x[200], y[200], D = 0, A;
32.     double a, yy;
33.     bool fine = true;
35.     while (cin>>n, cin.get() != EOF)
36.     {
37.         D = 0;
38.         for (int i=0; i<n; i++)
39.         {
40.             cin>>h[i]>>d[i];
41.             D += d[i];
42.         }
43.         x[0] = x[1] = 0;
44.         n *= 2;
46.         //计算出2n个最高点
47.         for (int i=2; i<n; i+=2)
48.         {
49.             x[i] = x[i-2] + d[i/2-1];
50.             x[i+1] = x[i] + d[i/2];
51.             y[i] = y[i+1] = h[i/2];
52.         }
53.     
54.         a=INF;
55.         for (int i=2; i<n-2; i++)
56.         {
57.             fine = true;
58.             A = y[i]/((D - x[i])*x[i]);
59.             for (int j=2; j<n-2; j++)
60.             {
61.                 //yy为在当前抛物线在x[i]处的高度， 若它比
62.                 //该处建筑物矮， 则舍弃当前抛物线
63.                 yy = A*x[j]*(D-x[j]);
64.                 if ( yy < y[j] )
65.                 {
66.                     fine = false;
67.                     break;
68.                 }
69.             }
70.             if (fine && a > A) {
71.                 a = A;
72.             }
73.         }
74.         double theta = atan(a*D);
75.         double angle = 180*theta/PI;
76.         double v0 = sqrt(4.9*(1+a*D*a*D)/a);
77.     
78.         printf("%.2f %.2f/n",angle, v0);
79.     }
80. }