ZOJ 3728 Collision
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=5074
题意:两个圆,小圆为实体,具有碰撞性。其中一个内含于另外一个,另有一枚硬币在大圆外,呈射线发射,求该硬币在大圆内的时间。
分析:
原先思路:圆心和直线的距离dist和R进行比较,R<dist则硬币和圆不相交。
若射线往圆的反方向射,则不会相交,但是用这种方法会判断出相交。
第二个错误的思想在于,虽然将直线转换为射线,没有求出交点,来求出时间t,但却没有判断t>0,若t<0的话,说明射线往反方向走
正确思路,若与大圆么没有两个交点,则时间为0,否则判断和小圆的交点,若没有两个交点,则距离为大圆两个交点距离,否则由于小圆反射
就是大圆和小圆的距离差
交点的数学原理:
圆:圆心为o,其半径为r,则||p-o||=r
射线:起点为p0,其速度方向为u,则p=p0+ut
若射线与圆有交点,则存在某个点pt,(p0+ut-o)^2=r^2
u^2t+2u(p0-o)t+(p0-o)^2-r^2=0;
这样就转换为判断delta来确定交点数量
求出来的t=(-b+-sqrt(b^2-4ac))就为时间
其交点就是p1=p0+ut
#include<stdio.h> #include<string.h> #include<math.h> #define eps 1e-9 struct point { double x,y; }; struct line { point a,b; }; struct circle { point p; double r; }; double distance(point p0,point p1) { return sqrt((p0.x-p1.x)*(p0.x-p1.x)+(p0.y-p1.y)*(p0.y-p1.y)); } double xmult(point p1,point p2,point p0) { return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y); } double disptoline(circle p,line l) { // printf("%lf",xmult(p.p,l.a,l.b)); //printf("%lf",distance(l.a,l.b)); return fabs(xmult(p.p,l.a,l.b))/distance(l.a,l.b); } int GetLineInsertion(line l,circle cir,point v,double &t1,double &t2) { double a=v.x,c=v.y; double b=l.a.x-cir.p.x,d=l.a.y-cir.p.y;//p0-o double f=2*(b*a+d*c);//2*t*u*(po-o) double e=a*a+c*c;//射线方向的平方 double g=b*b+d*d-cir.r*cir.r;//(po-o)^2-r^2 double delta=f*f-4*e*g; if(delta<0) return 0; if(delta==0) { t1=t2=-f/(2*e); } t1=(-f-sqrt(delta))/(2*e); t2=(-f+sqrt(delta))/(2*e); if(t1<0 || t2<0) return 0;//当时间为负数的时候,射线反方向执行 //虽然这样可以符合delta>0 return 2; } int main() { double R,RM,r; double k1,k2,k3,k4; point p1,v,p2; circle cir_R,cir_RM; cir_R.p.x=cir_R.p.y=cir_RM.p.x=cir_RM.p.y=0.0; line l; point t1,t2,t3,t4; while(scanf("%lf%lf%lf%lf%lf%lf%lf",&RM,&R,&r,&p1.x,&p1.y,&v.x,&v.y)!=EOF) { cir_R.r=R+r; cir_RM.r=RM+r; p2.x=p1.x+v.x*10000; p2.y=p1.y+v.y*10000; l.a.x=p1.x,l.a.y=p1.y; l.b.x=p2.x,l.b.y=p2.y; double dist=disptoline(cir_R,l);//圆心到直线的距离 int flag=GetLineInsertion(l,cir_R,v,k1,k2); if(flag==0) { printf("0.000\n"); } else { t1.x=p1.x+v.x*k1; t1.y=p1.y+v.y*k1; t2.x=p1.x+v.x*k2; t2.y=p1.y+v.y*k2; flag=GetLineInsertion(l,cir_RM,v,k3,k4); if(flag==0)//没有和小圆相交 { double dist_time; dist_time=distance(t1,t2); dist_time/=sqrt(v.x*v.x+v.y*v.y); printf("%.3lf\n",dist_time); } else { t3.x=p1.x+v.x*k3; t3.y=p1.y+v.y*k3; t4.x=p1.x+v.x*k4; t4.y=p1.y+v.y*k4; double dist_time; dist_time=distance(t1,t2)-distance(t3,t4); dist_time/=sqrt(v.x*v.x+v.y*v.y); printf("%.3lf\n",dist_time); } } } return 0; }