几何模板判断线段与直线相交模板

判断线段与直线是否相交

  对于线段相交的问题网上模板一大堆，但线段与直线相交的却很少，做几何计算的时候偶尔需要用到，如果没有相关的模板往往难以下手。这里就关于线段与直线相交的问题整理出一份模板，希望有所帮助。若有不足之处，请路过的大牛门不惜赐教。

  对于线段相交问题有两种版本：

       ①  包括端点相交和部分重合

       ②  不包括端点相交和部分重合

  那么对于直线与线段问题同样可以分为两种版本：

       ①  包括端点（线段端点）相交和重合

        ②  不包括端点和重合

    原理：判断线段两个端点是否在直线异侧

  预备函数：

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);    }    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;    }} ;struct line//线{    point a,b;};int sgn(double x){    if(fabs(x)<eps) return 0;    if(x<0) return -1;    return 1;}

  叉积

double mult(point p0,point p1,point p2)//叉积{    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);}

  版本一：  判断两点(p1,p2)是否在直线异侧，不包括端点和重合

bool opposite_side(point p1,point p2,point l1,point l2){    return mult(l1,p1,l2)*mult(l1,p2,l2)<-eps;}

  版本二：  判断线段l1与直线l2是否相交，包括端点相交

bool seg_inter_line(line l1,line l2){    return sgn((l1.a-l2.b)^(l2.a-l2.b))*sgn((l1.b-l2.b)^(l2.a-l2.b))<=0;}

  test:

#include<bits/stdc++.h>using namespace std;const double eps=1e-8;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);    }    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;    }} ;struct line//线{    point a,b;};int sgn(double x){    if(fabs(x)<eps) return 0;    if(x<0) return -1;    return 1;}double mult(point p0,point p1,point p2)//叉积{    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);}bool opposite_side(point p1,point p2,point l1,point l2){    return mult(l1,p1,l2)*mult(l1,p2,l2)<-eps;}bool seg_inter_line(line l1,line l2){    return sgn((l1.a-l2.b)^(l2.a-l2.b))*sgn((l1.b-l2.b)^(l2.a-l2.b))<=0;}bool inter(point s1,point e1,point s2,point e2)//判断线段相交{//    return opposite_side(s1,e1,s2,e2) && opposite_side(s2,e2,s1,e1);    return (max(s1.x,e1.x)>=min(s2.x,e2.x)) &&           (max(s2.x,e2.x)>=min(s1.x,e1.x)) &&           (max(s1.y,e1.y)>=min(s2.y,e2.y)) &&           (max(s2.y,e2.y)>=min(s1.y,e1.y)) &&           (mult(s1,s2,e1)*mult(s1,e1,e2)>0)&&           (mult(s2,s1,e2)*mult(s2,e2,e1)>0);}int main(){    line Seg,Line;    Seg.a.x=0,Seg.a.y=0;    Seg.b.x=1,Seg.b.y=1;    Line.a.x=0,Line.a.y=2;    Line.b.x=2,Line.b.y=0;    printf("test1: ");    if(opposite_side(Seg.a,Seg.b,Line.a,Line.b)) puts("opposite_side");    else puts("in_side!");    printf("test2: ");    if(seg_inter_line(Seg,Line)) puts("opposite_side");    else puts("in_side!");}