判断线段相交模板

 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;
12 }
13 
14 double Max(double a,double b){
15     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;//两线段不相交返回false
53 }