ACM 几何基础（2）

判断两条线段是否相交：

矢量

      如果一条线段的端点是有次序之分的话，那么这种线段就称为 有向线段，如果有向线段p1p2的起点p1在坐标的原点，则可以把它称为矢量p2

矢量的加减

      设二维矢量 P = (x1, y1), Q = (x2, y2)，则P + Q = (x1 + x2, y1 + y2), P - Q = (x1 - x2, y1 - y2)，且有P + Q = Q + P, P - Q = -(Q - P)

矢量叉积

      设矢量 P = (x1, y1), Q = (x2, y2)，则P * Q = x1 * y2 - x2 * y1;其结果是一个由(0, 0), P, Q, P + Q所组成的平行四边形的 带符号的面积，P * Q = -(Q * P), P * (- Q) = -(P * Q)

      叉积的一个非常重要的性质是可以通过它的符号来判断两矢量相互之间的顺逆时针关系：

            若 P * Q > 0，则 P 在 Q的顺时针方向

            若 P * Q < 0, 则 P 在 Q的逆时针方向

            若 P * Q = 0，则 P 与 Q共线，但不确定P, Q的方向是否相同

折线段的拐向判断

      如图，假设有折线段p0p1p2，要确定p1p2相对于p0p1是向左拐还是向右拐，可以通过计算(p2 - p0)*(p1 - p0)的符号来确定折线的拐向（点p2 - p0实际上就是向量p2，但这里要注意线段和矢量的区别）

判断点是否在线段上

      设点 Q = (x, y), P1 = (x1, y1), P2 = (x2, y2)，若(Q - P1)*(P2 - P1) = 0且min(x1, x2) <= x <= max(x1, x2) && min(y1, y2) <= y <= max(y1, y2)，则点Q在线段P1P2上

判断两线段是否相交

1）快速排斥试验

      设以线段 P1P2 为对角线的矩形为R，设以线段Q1Q2为对角线的矩形为T，若R、T不相交，则两线段不可能相交

      假设 P1 = (x1, y1), P2 = (x2, y2), Q1 = (x3, y3), Q2 = (x4, y4)，设矩形R的x 坐标的最小边界为minRX = min(x1, x2)，以此类推，将矩形表示为R = (minRX, minRY, maxRX, maxRY)的形式，若两矩形相交，则相交的部分构成了一个新的矩形F，如图所示，我们可以知道F的minFX = max(minRX, minTX), minFY = max(minRY, minTY), maxFX = min(maxRX, maxTX), maxFY = min(maxRY, maxTX)，得到F的各个值之后，只要判断矩形F是否成立就知道R和T 到底有没有相交了，若minFX > maxFX或minFY > maxFy则F无法构成，RT不相交，否则RT相交

 

 

 

2）跨立试验

      若 P1P2 跨立 Q1Q2，则 P1,P2 分别在 Q1Q2 所在直线的两端，则有(P1 - Q1)*(Q2 - Q1) * (Q2 - Q1)*(P2 - Q1) > 0，当(P1 - Q1)*(Q2 - Q1) = 0 时，说明(P1 - Q1)与(Q2 - Q1)共线，但由于已经经过快速排斥试验，所以Q1必为P1P2与 Q1Q2的交点，依然满足线段相交的条件，则跨立试验可改为：

      当 (P1 - Q1)*(Q2 - Q1) * (Q2 - Q1)*(P2 - Q1) >= 0，则P1P2跨立Q1Q2，当Q1Q2也跨立P1P2的时候，则P1P2相交

      （注意式子中被隔开的*代表两个叉积的值的相乘，而另外的两个*则代表计算矢量叉积）

 

 

 

 

判断线段与直线是否相交：

首先在直线上随便取两点做为线段l1;

再与线段L2进行判断，因为是直线所以现在只需判断L2是否跨立l1就行了。

 

代码如下：

bool inter_line(line l1,line l2)//判断直线与线段是否相交

{

    return sgn((l2.s-l1.e)^(l1.s-l1.e))*sgn((l2.e-l1.e)^(l1.s-l1.e))<=0;

}

 

 

代码实现：

再这之前可以先去学下，C++的运算符重载部分

运算符重载详解  http://blog.chinaunix.net/uid-21411227-id-1826759.html

#include <iostream>

#include <cstdio>

#include <cmath>

using namespace std;

const double eps=1e-8;

const double PI=acos(-1.0);

int sgn(double x)

{

    if(fabs(x)<eps) return 0;

    if(x < 0) return -1;

    else return 1;

}

 

struct point

{

    double x,y;

    point(){}

    point (double _x,double _y)

    {

        x=_x,y=_y;

    }

    point operator -(const point &b)const

    {

        return point(x-b.x,y-b.y);

    }

    point operator +(const point &b)const

    {

        return point (x+b.x,y+b.y);

    }

    double operator ^(const point &b)const //叉积

    {

        return x*b.y-y*b.x;

    }

    double operator *(const point &b)const //点积

    {

        return x*b.x+y*b.y;

    }

    void transxy(double B)  //绕原点旋转角度B（弧度值），后的x，y

    {

        double tx=x,ty=y;

        x= tx*cos(B) - ty*sin(B);

        y= tx*sin(B) + ty*cos(B);

    }

};

 

struct line//线

{

    point s,e;

    line(){}

    line(point _s,point _e)

    {

        s=_s,e=_e;

    }

};

 

double dist(point a,point b)//两点间距离

{

    return sqrt((a-b)*(a-b));

}

 

bool inter(line l1,line l2)//判断线段相交

{

    return max(l1.s.x,l1.e.x)>=min(l2.s.x,l2.e.x)&&

    max(l2.s.x,l2.e.x)>=min(l1.s.x,l1.e.x)&&

    max(l1.s.y,l1.e.y)>=min(l2.s.y,l2.e.y)&&

    max(l2.s.y,l2.e.y)>=min(l1.s.y,l1.e.y)&&

    sgn((l2.s-l1.e)^(l1.s-l1.e))*sgn((l2.e-l1.e)^(l1.s-l1.e))<=0&&

    sgn((l1.s-l2.e)^(l2.s-l2.e))*sgn((l1.e-l2.e)^(l2.s-l2.e))<=0;

}

 

bool inter_line(line l1,line l2)//判断直线与线段是否相交

{

    return sgn((l2.s-l1.e)^(l1.s-l1.e))*sgn((l2.e-l1.e)^(l1.s-l1.e))<=0;

}

 

int main()

{

    double ax,ay,bx,by,cx,cy,dx,dy;

    while(scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&ax,&ay,&bx,&by,&cx,&cy,&dx,&dy)!=EOF)//分别输入两条线段的首尾端点

    {

        point k1(ax,ay);

        point k2(bx,by);

        point k3(cx,cy);

        point k4(dx,dy);

        line l1(k1,k2);

        line l2(k3,k4);

        if(inter(l1,l2))

            printf("Yes\n");

        else

            printf("No\n");

}

return 0;

}