计算几何模板（2）：线与线

比点与向量部分要难一点，会有一定的技巧与证明。

---

1.首先是直线。

小学生都知道两点确定一条直线。

所以我们不用两点式，用点斜式。

一个直线上的点+一个向量=一条直线。

struct Line
{
    Point x;
    Vector v;
    Line(){}
    Line(Point x,Vector v):x(x),v(v){}
};

直线

2.两直线求交点

我们知道两直线是这样的：

求点$C$。

现在来做几条辅助线：

现在新产生了三个三角形，两大一小。

两个大的面积相等，因为同底等高。

所以面积可求，$AC$可求。

Point L_L(Line a,Line b)
{
    double t = (b.v^(b.x-a.x))/(b.v^a.v);
    return a.x+a.v*t;
}

两直线交点

 3.点到直线距离

有标准式方程就可以套公式了，但是我们只有（假的）点斜式。

其实更简单，搞个三角形叉积求面积，然后求高就好了。

就像这样：

（叉积求的是有向面积，结果要加绝对值）

double Dis_P_L(Point p,Line l)
{
    return fabs((l.v^(p-l.x))/lth(l.v));
}

点->直线

4.点到线段

（线段直接用端点给出）

端点处有钝角=垂足不在线段上。

钝角直接用$sin$值判。

double Dis_P_S(Point p,Point a,Point b)
{
    if(a==b)return lth(p-a);
    Vector v1 = b-a,v2 = p-a,v3 = p-b;
    if(v1*v2<=0)return lth(v2);
    else if(v1*v3>=0)return lth(v3);
    else return (v1^v2)/lth(v1);
}

点->线段

5.过定点求直线垂足（非人语

依然是面积法，与上面那个差不多（？）就不证了。

Point P_L(Point p,Line l)
{
    return p+l.v*(((p-l.x)*l.v)/(l.v*l.v));
}

垂足

6.线段是否相交（交点不在端点上）

相交和不相交：

然后连起来：

所以判对于任意一条线段，另一条线段的两端点是否在其两侧即可。

bool S_S(Point a1,Point a2,Point b1,Point b2)
{
    double d1 = (a2-a1)^(b1-a1),d2 = (a2-a1)^(b2-a1);
    double d3 = (b2-b1)^(a1-b1),d4 = (b2-b1)^(a2-b1);
    return dcmp(d1)*dcmp(d2)<0&&dcmp(d3)*dcmp(d4)<0;
}

线段相交

7.点在线段上（不能在端点上）

叉积判点在直线上，点积判端点在其两侧。

bool P_S(Point p,Point a,Point b)
{
    Vector A = a-p,B = b-p;
    return dcmp(A*B)==0&&dcmp(A^B)<0;
}

点与线段

8.求多边形有向面积

选一个点，然后把多边形分成$n-2$个三角形然后叉积。

typedef vector<Point> Pol;
double S_(Pol p)
{
    double ret = 0;
    for(int i=1,lim=(int)p.size();i<lim-1;i++)
        ret+=((p[i]-p[0])^(p[i+1]-p[0]));
    return ret/2;
}

多边形面积