现场赛的时候,读了三遍都读不懂的一个题,当时还是队友跟我解释的题意,虽然那时有想法,不过是另一个队友搞掉的。
题目:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=5074
大意是有两个同心圆,圆心都是原点,另外有个硬币,硬币也看成一个圆,初始时跟两个同心圆没有相交。
然后给个硬币运动的方向(向量表示),问硬币在运动过程中,跟大圆有任何一点接触的时间是多长。如果碰到里面的小圆,会发生碰撞,但无能量损耗。
小圆半径Rm,大圆半径R,硬币半径r,圆心(x,y),速度方向(vx,vy)
分三种情况考虑:
1、根本不会接触到大圆,直接输出0.000
2、接触到大圆,但是没有跟小圆碰撞,这时答案是跟大圆接触的时间
3、跟小圆发生碰撞,因为这是个对称的过程,所以可以计算硬币刚进入大圆的时间以及碰到小圆的时间,时间差的两倍就是答案了。
先将圆相交碰撞转换成圆心的距离,所以当硬币跟大圆相交时,圆心距离d1<=R+r,跟小圆碰撞时d2=Rm+r
所以在某时刻t,圆心的位置就是(x+vx*t, y+vy*t)
对于大圆相交,我们可以列出一个关于t的一元二次不等式,
d1 = (x+vx*t)^2+(y+vy*t)^2 = (vx^2+vy^2)*t^2 + 2*(x*vx+y*vy)*t + x^2+y^2<=(R+r)^2
计算出△=b*b-4ac,判断是否有实根,没有就证明不接触,如果有实根,两个根对应的时刻正是开始跟大圆接触还有刚好离开大圆的两个时刻,因为硬币一开始不跟大圆接触,所以两个时刻的正负是一致的,如果是负数说明是反向移动接触到的,实际还是没有接触;
在前面的时刻是正的基础上,用同样的方式再判断是否可以跟小圆碰撞,分别对应上述的第2、3情况。
#include<cstdio>
#include<cmath>
double rm, r0, r, x, y, vx, vy, d, a, b, c, t1, t2, t3, t4;
int main(){
while(~scanf("%lf %lf %lf %lf %lf %lf %lf", &rm, &r0, &r, &x, &y, &vx, &vy)){
rm+=r;
r0+=r;
a = vx*vx+vy*vy;
b = 2*(x*vx+y*vy);
c = x*x+y*y-r0*r0;
d = b*b-4*a*c;
if(d<=0){
puts("0.000");
}
else{
t1 = (b*(-1)-sqrt(d))/(2*a);
if(t1<=0){//不跟大圆接触
puts("0.000");
}
else{
t2 = (b*(-1)+sqrt(d))/(2*a);
c = x*x+y*y-rm*rm;
d = b*b-4*a*c;
if(d<=0){//不和小圆碰撞
printf("%.3lf\n", t2-t1);
}
else{
t3 = (b*(-1)-sqrt(d))/(2*a);
printf("%.3lf\n", 2*(t3-t1));
}
}
}
}
return 0;
}