计算几何之线段性质（二）：求线段交点

求线段的交点即求线段交点的坐标。若已确定两线段相交，可以通过解析几何中的直线方程来求解交点，这里介绍另外的一种方法。叉积multi（p1,p2,p0）可以看做是由点P1，P2，P0和P1+P2,这四个点围成的平行四边形的面积（如下图所示），及三角形P0P1P2的面积S=1/2 multi（P1，P2，P0）

在下图中，线段AB和CD相交于点P。

我们分别从D点和C点做AB的垂线DD1和CC1，由于△DD1P∽△CC1P，因此我们有|DD1| / |CC1|=|DP| / |CP|。又因为S△ABD=DD1*AB/2，S△ABC=CC1*AB/2，所以：

|DP| / |CP|=S△ABD/S△ABC=|AD × AB| / |AC × AB|=|multi（D,B,A）| / |multi（C,B,A）|

又因为|DP| / |CP|=（XD-XP）/（XP-XC）=（YD-YP）/（YP-YC）

所以有：

XP=（S△ABD * XC + S△ABC * XD）/（S△ABD + S△ABC）

    =（multi（D,B,A）* XC - multi（C,B,A）* XD）/（multi（D,B,A）- multi（C,B,A））

YP=（S△ABD * YC + S△ABC * YD）/（S△ABD + S△ABC）

    =（multi（D,B,A）* YC - multi（C,B,A）* YD）/（multi（D,B,A）- multi（C,B,A））

由此得出计算线段u和v交点的算法，参考代码如下：

#include <cstdio>
using namespace std;
typedef struct node
{
    double x,y;
}point;
typedef struct line
{
    point start,end;
}segment;
double multi(point p1,point p2,point p0)
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
point intersection(segment u,segment v)
{
    point p;
    p.x=(multi(v.end,u.end,u.start)*v.start.x-multi(v.start,u.end,u.start)*v.end.x)/(multi(v.end,u.end,u.start)-multi(v.start,u.end,u.start));
    p.y=(multi(v.end,u.end,u.start)*v.start.y-multi(v.start,u.end,u.start)*v.end.y)/(multi(v.end,u.end,u.start)-multi(v.start,u.end,u.start));
    return p;
}
int main()
{
    segment u,v;
    scanf("%lf%lf%lf%lf",&u.start.x,&u.start.y,&u.end.x,&u.end.y);
    scanf("%lf%lf%lf%lf",&v.start.x,&v.start.y,&v.end.x,&v.end.y);
    point p=intersection(u,v);
    printf("%lf  %lf\n",p.x,p.y);
    return 0;
}