练习二 1004 Toxophily

<span style="font-size:14px;">题目：</span>

<span style="font-size:14px;"></span><div class="panel_content">The recreation center of WHU ACM Team has indoor billiards, Ping Pang, chess and bridge, toxophily, deluxe ballrooms KTV rooms, fishing, climbing, and so on.<br>We all like toxophily.<br><br>Bob is hooked on toxophily recently. Assume that Bob is at point (0,0) and he wants to shoot the fruits on a nearby tree. He can adjust the angle to fix the trajectory. Unfortunately, he always fails at that. Can you help him?<br><br>Now given the object's coordinates, please calculate the angle between the arrow and x-axis at Bob's point. Assume that g=9.8N/m. <br></div><div class="panel_bottom"> </div>
<div class="panel_title" align="left">Input</div><div class="panel_content">The input consists of several test cases. The first line of input consists of an integer T, indicating the number of test cases. Each test case is on a separated line, and it consists three floating point numbers: x, y, v. x and y indicate the coordinate of the fruit. v is the arrow's exit speed.<br>Technical Specification<br><br>1. T ≤ 100.<br>2. 0 ≤ x, y, v ≤ 10000. <br></div><div class="panel_bottom"> </div>
<div class="panel_title" align="left">Output</div><div class="panel_content">For each test case, output the smallest answer rounded to six fractional digits on a separated line.<br>Output "-1", if there's no possible answer.<br><br>Please use radian as unit. <br></div><div class="panel_bottom"> </div>
<div class="panel_title" align="left">Sample Input</div><div class="panel_content"><pre><div style="FONT-FAMILY: Courier New,Courier,monospace">3<br>0.222018 23.901887 121.909183<br>39.096669 110.210922 20.270030<br>138.355025 2028.716904 25.079551<br></div>

Sample Output

  
  

   
   1.561582<br>-1<br>-1<br>
  
  

<span style="font-size:14px;">题意：bob在（0，0）点想射到在点（x,y）的水果，初速度是v，且g=9.8,求能射到苹果的最小角度。</span>

<span style="font-size:14px;">解题思路：由题意可写出物理关系式y=tan(t)*x+g*x*x/(2*v*v*cos(t)*cos(t)),t就是速度方向与x的夹角，由关系式可知在（0，π/2）内先增后减</span>

<span style="font-size:14px;">可以先用三分求出t的最大值，然后再用二分在（0，t0）内求出符合要求的值。</span>

<span style="font-size:14px;">感想：二分加三分，三分用的还是不好。</span>

<span style="font-size:14px;">代码：</span>

<span style="font-size:14px;">#include<iostream>  
#include<cstdio> 
#include<iomanip> 
#include<cmath>   
using namespace std;   
const double PI=acos(-1.0);  
const double g=9.8;  
double x,y,v;  
double f(double t)  
{  
    return x*tan(t)-g*x*x/2/(v*v*cos(t)*cos(t));  
}  
double ef(double left,double right)
{
	double mid;
	while(right-left>=1e-9)
	{
		mid=(left+right)/2;
		if(f(mid)>=y)
			right=mid;
		else left=mid;
	}
	return left;
}
double sf(double left,double right)  
{  
    double mid,midmid;
	while(right-left>=1e-9) 
    {  
        mid=left+(right-left)/3;  
        midmid=right-(right-left)/3;  
        if(f(mid)<=f(midmid))  
            left=mid;  
        else  
            right=midmid;  
    }  
    return (mid+midmid)/2;  
}  
  
int main()  
{  
    int t;  
    double t0;  
    cin>>t;  
    while(t--)  
    {  
        cin>>x>>y>>v;  
        t0=sf(0,PI/2);  
        if(f(t0)<y)  
            cout<<"-1"<<endl;
		else cout<<fixed<<setprecision(6)<<ef(0,t0)<<endl; 
    }  
    return 0;  
}  </span>