1 /* 2 部分重合，端点重合，点边重合都会返回true 3 */ 4 struct PT{ 5     double x,y;//点的坐标 6     PT(double x = 0, double y = 0) : x(x), y(y){} 7     //bool operator 8 };//线段的2个端点，MAX表示线段条数的上限 9 10 double Min(double a,double b){11     return a
     
      b?a:b;16 }17 18 double Dir(PT p1,PT p2,PT p0){
   //判断线段p0p2相对于线段p0p1的位置19     //向量p0p1叉乘向量p0p2，转化为行列式进行计算20     //结果若大于0，则p0p2沿顺时钟方向旋转可得到与p0p1同方向的向量21     //结果若小于0，则p0p2沿逆时钟方向旋转可得到与p0p1同方向的向量22     //结果若等于0，则p0p2与p0p1共线23     return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);24 }25 26 int OnSgm(PT p1,PT p2,PT p0){27     double minx,maxx,miny,maxy;28     minx=Min(p1.x,p2.x);29     maxx=Max(p1.x,p2.x);30     miny=Min(p1.y,p2.y);31     maxy=Max(p1.y,p2.y);32     return minx<=p0.x&&p0.x<=maxx&&miny<=p0.y&&p0.y<=maxy;33 }34 35 bool Intersect(PT p1,PT p2,PT p3,PT p4){
   //两线段相交有两种情况36     double d1,d2,d3,d4;37     d1=Dir(p4,p1,p3);38     d2=Dir(p4,p2,p3);39     d3=Dir(p2,p3,p1);40     d4=Dir(p2,p4,p1);41     if(d1*d2<0&&d3*d4<0)//两线段互相跨过42         return 1;43     //一条直线有一个端点在另外一条直线上44     if(d1==0&&OnSgm(p3,p4,p1))45         return 1;46     if(d2==0&&OnSgm(p3,p4,p2))47         return 1;48     if(d3==0&&OnSgm(p1,p2,p3))49         return 1;50     if(d4==0&&OnSgm(p1,p2,p4))51         return 1;52     return 0;//两线段不相交返回false53 }