15.计算几何：点、直线、线段的位置关系及其它关系

直线的表示

直线有多种表示方法。编程时，可以灵活使用这些方法：

struct Line
{
    Point p1,p2;    //线上的两个点
    Line(){}
    Line(Point p1,Point p2):p1(p1),p2(p2){}
    
    
   //根据一个点和倾斜角angle——>确定直线
   //其中：0<=angle<pi （就是0到180度）
    Line(Point p,double angle)
    {
        p1 = p;
        if(sgn(angle – pi/2) == 0)
        {
            p2 = (p1 + Point(0,1));
        }
        else
        {
            p2 = (p1 + Point(1,tan(angle)));
        }
    }
    

   //ax+by+c=0
    Line(double a,double b,double c)
    {
        
        if(sgn(a) == 0)
        {
            //y=-c/b
            //p1是把x=0带入y=-c/b，得-c/b
            //p2是把x=1带入y=-c/b，得-c/b
            p1 = Point(0,-c/b);
            p2 = Point(1,-c/b);
        }
        else if(sgn(b) == 0)
        {
            //x=-c/a
            p1 = Point(-c/a,0);
            p2 = Point(-c/a,1);
        }
        else
        {
            //p2-p1刚好是直线：y=-a/b*x-c/b的斜率
            //p1是把x=0带入y=-a/b*x-c/b，得-c/b
            //p2是把x=1带入y=-a/b*x-c/b，得(-c-a)/b
            p1 = Point(0,-c/b);
            p2 = Point(1,(-c-a)/b);
        }
    }
};

---

线段的表示

可以用两个点来表示线段，起点是p1，终点是p2。
直接用直线的数据结构来定义线段：typedef Line Segment；

---

位置关系

1. 点和直线的位置关系：点在直线左侧、在右侧、在直线上。

2. 点和线段的位置关系：点 p 是否在线段 v 上

---

点和直线的其它关系

1. 点到直线的距离

---

2. 点关于直线的对称点

补充：【常规做法】

简单做法：
先求点p在直线上的投影q，再用中点公式求对称点p’

点在直线上的投影
Point Point_line_proj(Point p, Line v)
{
    double k=Dot(v.p2-v.p1,p-v.p1)/Len2(v.p2-v.p1);
    return v.p1+(v.p2-v.p1)*k;
}

Point Point_line_symmetry(Point p, Line v)
{
    Point q = Point_line_proj(p,v);
    return Point(2*q.x-p.x,2*q.y-p.y);
}

代码解释

---

3. 点到线段的距离

我们设现在有一个点p以及线段AB，那么点p到线段AB的距离只要从3个距离中取最小值即可。

---

4. 两条直线的位置关系

两条直线的位置关系可能有：
重合
平行
相交

int Line_relation(Line v1, Line v2)
{
	//判断直线是否相交 为0：不相交
    if(sgn(Cross(v1.p2-v1.p1,v2.p2-v2.p1)) == 0)
    {
           if(Point_line_relation(v1.p1,v2)==0) return 1;  //1 重合
          else return 0;                                   //0 平行
    }
    return 2;                                              //2 相交
} 

---

5. 判断两个线段是否相交

关系两线段是否相交，我们可以利用叉积有正负的特点。
如果一条线段的两端在另一条线段的两侧，那么两个端点与另一线段的产生的两个叉积正负相反，也就是说两个叉积相乘为负。如果两条线段互相满足这一点，那么就是相交的。

//Line1:ab,  Line2:cd
bool Cross_segment(Point a,Point b,Point c,Point d)
{
    double c1=Cross(b-a,c-a),c2=Cross(b-a,d-a);
    double d1=Cross(d-c,a-c),d2=Cross(d-c,b-c);
    return sgn(c1)*sgn(c2)<=0 && sgn(d1)*sgn(d2)<=0;  //1相交；0不相交
}

---

6. 求两条线段的交点

这个简单：先判断两条线段是否相交。若相交，问题转化成两条直线求交点；若不相交，则无交点