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

判断线段与直线是否相交

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

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

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

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

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

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

        ②  不包括端点和重合

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

  预备函数：

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!");
}