题目大意:有一个与地面垂直的管子,管口与地面相距H,管子里面有很多弹性球,从t=0时,第一个球从管口求开始下落,然后每1s就会又有球从球当前位置开始下落,球碰到地面原速返回,球与球之间相碰会发生完全弹性碰撞(各自方向改变,速率变为相碰时另一个球的速率)问最后所有球的位置?
这一题又是一道弹性碰撞模型的题目,和Liner World有点像,但是这一题不要求输出球的名字,而是要你求得各个球的高度
这道题我们只用明白一个道理就很方便了:
首先我们来看一个球的情况:
一个球从H的高度下落,碰到地面原速返回,那么他在t时刻的位置由下列方程决定:
t=sqrt(2*H/g);
k=(int)(T/t);
Hi=H-(T-k*t)*(T-k*t)*g/2; (当K是偶数)
Hi=H-(k*t+t-T)*(k*t+t-T)*g/2; (当K是奇数)
如果有两个球:
如果第一个球没有碰地,那么两个球的高度将由Hi=H-(T-k*t)*(T-k*t)*g/2唯一确定
如果第一个球碰地了,那么在接下来的某一时刻两个球一定会相碰且交换速度,可以看成是擦身而过,但问题还没完,如果相碰后不再相碰,要确定高度,那么究竟怎么确定呢?
我们知道因为各个球可以看成独立的,根据题意,第2个球的高度应该要加上两个球的半径并除以100.0(单位是cm看清楚!),然后交换速度过后我们虽然2和1交换了状态,但是状态不变,我们还是可以把他们当成一个球看,但是同时注意,这两个球是有半径的!(一开始我没有想到半径打死都想出来为什么,模拟了一下才明白),最后答案为从低到高排列的小球逐次加上0和2r的高度即可(如果存在交换一定要排列,否则两个球就会重叠了就不对了)。
把上面这两种情况推广到多个球,我们就可以得到解决方法了
#include <iostream>
#include <algorithm>
#include <functional>
#include <math.h> using namespace std; double cal(const double, const double, const double);
const double g = 10.0;
static double h[]; int main(void)//很巧妙的一道物理模型题
{
int case_sum, N;
double H, R, T, seg_t;
scanf("%d", &case_sum); while (case_sum--)
{
scanf("%d%lf%lf%lf", &N, &H, &R, &T);
seg_t = sqrt( * H / g); for (int i = ; i < N; i++)
h[i] = cal(seg_t, T - i, H);//每一个间隔减少1s
sort(h, h + N);
for (int i = ; i < N; i++)
printf("%.2f ", h[i] + * R*i / 100.0);
printf("\n");
}
return EXIT_SUCCESS;
} double cal(const double t, const double T, const double H)
{
if (T < )
return H;//相当于还没开始下落就直接
int k =(int)(T / t);//看究竟是上升阶段还是下降阶段。
if (k % == )//下降阶段
return H - g*(T - k*t)*(T - k*t) / ;
else
return H - g*(k*t + t - T)*(k*t + t - T) / ; }