题目编号:1004
题目大意:一个人射箭,在(0,0)点开始射箭,往(x,y)处射,有个g的重力加速度,给你个初速度v求最小的射中角度。
解题思路:额,初看题目,怎么说呢,不像是编程题,更像个物理数学计算的题目,看给出的输入例子,可以看到小数点后面有好多。这是个信号,要你保留到小数点几位几位, 而且很容易让人联想到用二分查找的方法来求解。至于求解过程就很简单了,大体的框架就是运用物理学的公式计算保留位数,然后再一个循环判断,二分查找,直到 达到要求,如果找不到直接返回NO,所以关键是物理公式要写对就行了。用公式,根据正交分解坐标系,得出方程的通式。想 x^2*g/(2*v^2)*tan^2(ß) - x*tan(ß) +y + x^2*g/(2*v^2) = 0; 即:a = g*pow(x,2)/(2*pow(v,2)); b = -x; c = y + g*pow(x,2)/(2*pow(v,2));然后根据求根的公式求出来根就行了,但要注意后面的讨论,有无解情况。
解题感想:感想呢,怎么说呢,挺惭愧的,物理公式忘得七七八八了,还是看老师的公式才隐约记得知道这么回事,编程的一些高级算法果然和数学和物理搭边啊!这个是老师的 代码,我只是看懂了,偷偷模仿写的提交的,一开始确实写不出来。
#include <iostream>
#include <stdio.h>
#include <math.h>
using namespace std;
int main()
{
int t;
double a,b,c,angle,z;
double x,y,v,g = 9.8,T,ans1,ans2;
scanf("%d",&t);
while(t--)
{
scanf("%lf%lf%lf",&x,&y,&v);
if(x==0&&y==0)
printf("0\n");
else if(x==0&&y>0)
printf("90\n");
else
{
a = g*pow(x,2)/(2*pow(v,2));
b = -x;
c = y+a;
T = pow(b,2) - 4*a*c;
angle = 0;
if(T<0)
printf("-1\n");
else
{
ans1 = ((-b)+pow(T,1.0/2))/(2*a);
ans2 = ((-b)-pow(T,1.0/2))/(2*a);
if(ans1>=0) angle = atan(ans1);
if(ans2>=0)
{
z = atan(ans2);
if(z<angle) angle = z;
printf("%.6f\n",angle);
}
if(ans1<0&&ans2<0)
printf("-1\n");
}
}
}
return 0;
}