HDOJ 题目2603 Wiskey's Power（数学）

Wiskey's Power

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 482    Accepted Submission(s): 210

Problem Description

Come back school from the 33 ^rdACM / ICPC Asia ChenDu, everyone is exhausted, in particular the captain Wiskey. As saying that night, Wiskey drink a lot of wine, just as he blurred, fall to sleep. All of a sudden, Wiskey felt a slight discomfort in the chest, and then vomiting, vomitted all over. The next day, Wiskey is extremely sluggish, vomit still on the train.

Your task is to calculate the coordinates of vomit.
We assume that the quality of vomit is m, and its size can be ignored.
As the figure below:

The vomit start from the blue point A, whose speed is V, and its angle with X-axis is a. If the vomit hit the ceiling, then its value of the speed don't changed and if before the collision the angle of the speed with X-axis is b, then after the collision the angle of the speed with X-axis is b , too.
Ignore air resistance, acceleration due to gravity g = 9.87m/s ², calculate and output Q.
(you can assume that the vomit will not fall to the higher berth)

 

Input

Each line will contain three numbers V , m and a (0 <= a < 90)described above.
Process to end of file.

 

Output

For each case, output Q in one line, accurate up to 3 decimal places.

 

Sample Input

  
  

   
   100 100 45
  
  

 

Sample Output

  
  

   
   3.992
  
  

 

Author

WhereIsHeroFrom

 

Source

HDU 1st “Vegetable-Birds Cup” Programming Open Contest

 

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2603

ac代码

#include<stdio.h>
#include<math.h>
int main()
{
	double v,m,a;
	double pi=acos(-1.0);
	while(scanf("%lf%lf%lf",&v,&m,&a)!=EOF)
	{
		double vx,vy,t,sm,g=9.87,vm,t1,t2,vm2,s;
		a=a/180*pi;
		vy=v*sin(a);
		vx=v*cos(a);
		sm=vy*vy/(2*g);
		
		if(sm>0.5)
		{
			double t1=2*(vy-sqrt(vy*vy-2*0.5*g))/g;  
            double t2=(sqrt(vy*vy+2*g*3)-vy)/g;  
			t=t1+t2;
		}
		else
		{
			t=2*vy/g+(sqrt(vy*vy+2*g*3)-vy)/g;
		}
		s=vx*t;
		printf("%.3lf\n",s);
	}
}