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
");
        }
        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
",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
",dist_time);
            }
        }
    }
    return 0;
}