Time Limit:3000MS Memory Limit:0KB
其实就是暴力所有情况,对每个弹跳次数搜索所有障碍高。
对于弹跳次数b,那么每次抛物线向前移动d/(b+1)次,设为x0。那么对于每个障碍将他们的水平距离 x,x%x0即可将他们都移到过原点的抛物线来计算。因为x0不一定是整数,需要自己写一个去模函数。
设a为障碍水平距离,b为障碍高,则vy2=g*x0^2*b / 2a(x0-a) vx2=g^2*x0^2 / 4vy2
对于弹跳次数b,算出最大的vy2+vx2作为其v2,然后在所有的v2里面找到最小的就是结果。
需要注意的是如果vy2小于vx2,那么出射角小于45度,这种情况下最小的v是在45度时取到的,即vx2=vy2=g*x0 / 2
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
#define exp 0.00000001
#define inf 1e17
struct ob
{
double x,h;
};
double mod(double a,double b)
{
if (a<b) return a;
int t=floor(a/b);
return (a-b*t);
}
int main()
{
int n,b,d;
double ans;
while(cin>>d>>n>>b)
{
ob p[11];
for(int i=0;i<n;i++)
cin>>p[i].x>>p[i].h;
double maxx,minn=inf,vx2,vy2;
for(int i=0;i<=b;i++)
{
maxx=0;
double x0=(d*1.0)/(i+1.0);
bool yep=1;
for(int j=0;j<n;j++)
{
double a=mod(p[j].x,x0),b=p[j].h;
if (fabs(a-0)<exp) {yep=0;break;}
maxx=max(maxx,x0*x0*b/(2.0*a*(x0-a)));
}
if (!yep) continue;
vy2=maxx;
vx2=x0*x0/(4.0*vy2);
if (vx2>vy2) vx2=vy2=x0/2.0;
minn=min(vy2+vx2,minn);
}
ans=sqrt(minn);
printf("%.5lf\n",ans);
}
}