(p1-p0) x (p3-p2)结果的意义:

   正: <p0p1>在<p2,p3>顺时针(0,pi)内

   负: <p0p1>在<p2,p3>逆时针(0,pi)内

   0 : <p0p1>, <p2p3>共线,夹角为0或pi

double xmult(point p0,point p1,point p2){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

   可用于判p2与<p0 p1>的关系

 

double xmult(point p0,point p1,point p2,point p3){
    return (p1.x-p0.x)*(p3.y-p2.y)-(p3.x-p2.x)*(p1.y-p0.y);
}

   可用于判<p0p1>, <p2p3>位置关系

double xmult(line u,line v){
    return (u.b.x-u.a.x)*(v.b.y-v.a.y)-(v.b.x-v.a.x)*(u.b.y-u.a.y);
}

//计算两直线交点,注意事先判断直线是否平行!
point intersection(line u,line v){
    point ret=u.a;
    double t=((u.a.x-v.a.x)*(v.a.y-v.b.y)-(u.a.y-v.a.y)*(v.a.x-v.b.x))
            /((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));
    ret.x+=(u.b.x-u.a.x)*t;
    ret.y+=(u.b.y-u.a.y)*t;
    return ret;
}

int relat_and_intersection(line u,line v,point &ret)
{
    /**//*
    LINE 重合
    NONE 平行
    INTERSECT 相交并返回交点
    */
    double a1,b1,c1,a2,b2,c2;
    //系数
    a1=u.a.y-u.b.y;
    b1=-(u.a.x-u.b.x);
    c1=u.a.x*u.b.y-u.a.y*u.b.x;

    a2=v.a.y-v.b.y;
    b2=-(v.a.x-v.b.x);
    c2=v.a.x*v.b.y-v.a.y*v.b.x;

    if(fabs(a1*b2-a2*b1)<eps)
    {
        if( fabs(a1*c2-a2*c1)<eps && fabs(b1*c2-b2*c1)<eps )
            return LINE;
        else return NONE;
    }
    else 
    {
        point ret;
        ret.x=(b1*c2-b2*c1)/(a1*b2-a2*b1);
        ret.y=(a2*c1-a1*c2)/(a1*b2-a2*b1);
        return INTERSECT;
    }
}