线段ab与cd相交当且仅当a、b在直线cd两侧且c、d在直线ab两侧;通过叉积符号判断点在线段的左右侧,cross(p,a,b)>0表示p在ab左侧。

用叉积判断线段相交的核心逻辑
两条二维线段是否相交,不能只靠求直线交点再判断是否在线段范围内——那样容易因浮点误差或共线情况出错。真正稳定的做法是用向量叉积做方向性判断,即「跨立实验」:线段 AB 与 CD 相交 ⇔ A、B 在直线 CD 两侧,且 C、D 在直线 AB 两侧。
关键点在于:叉积符号反映点相对于有向线段的左右方位,cross(p, a, b) 表示向量 a→p 和 a→b 的二维叉积(即 (p.x - a.x) * (b.y - a.y) - (p.y - a.y) * (b.x - a.x))。结果 > 0 表示 p 在 ab 左侧,
实现 segments_intersect 的完整步骤
写一个健壮的判断函数,必须同时处理「规范相交」(内部交叉)和「端点接触」(如一个端点落在另一线段上)两种常见需求。多数工程场景默认只要求规范相交,但需明确约定。
- 先封装
cross函数,输入三个Point(假设为struct { double x, y; }),返回double - 定义
on_segment(p, a, b)判断点 p 是否在线段 ab 上:先用cross(p, a, b) == 0确保共线,再用点积dot(p-a, p-b) 或坐标范围检查确保不在线段延长线上 - 主逻辑:计算
c1 = cross(c, a, b),c2 = cross(d, a, b),c3 = cross(a, c, d),c4 = cross(b, c, d) - 规范相交条件:
c1 * c2 (严格异号,排除共线) - 若需包含端点接触,则额外检查:
(c1 == 0 && on_segment(c, a, b)) || (c2 == 0 && on_segment(d, a, b)) || (c3 == 0 && on_segment(a, c, d)) || (c4 == 0 && on_segment(b, c, d))
浮点数比较和零值陷阱
直接用 == 0 判断叉积是否为零,在浮点下极不可靠。必须引入容差(epsilon),但 epsilon 选多大很讲究。
- 不要用固定小值如
1e-9:当坐标值很大(比如 1e6 级别)时,叉积本身可能高达 1e12,此时1e-9比较毫无意义 - 推荐用相对误差:定义
const double eps = 1e-10;,然后用fabs(x) 做归一化,或更简单地用 <code>fabs(x) - 所有
cross结果、所有on_segment中的点积和坐标比较,都要走同一套 epsilon 判定逻辑,否则行为不一致 - 如果输入坐标是整数,可改用
long long做精确叉积(避免 double 舍入),此时== 0安全
性能与边界场景提醒
这个判断本身是 O(1),但实际使用中容易忽略几个隐性成本和失效点:
- 重复调用时,别在函数内反复构造
Point临时对象;传引用或直接传double x1, y1, x2, y2, x3, y3, x4, y4更快 - 共线但不重叠的线段(如 AB 和 CD 同向但分离)会被正确判为不相交,但若没写好
on_segment,可能误判为相交 - 两线段垂直或接近垂直时,叉积数值稳定;但接近平行时,叉积趋近于 0,此时 epsilon 容差的影响被放大,务必用前述动态容差
- 图形学中常需批量判断(如碰撞检测),建议预计算每条线段的 bounding box 并先做 AABB 快速剔除,避免大量调用叉积
最易被忽略的是:你写的 on_segment 是否真的能区分「点在线段上」和「点在直线延长线上」——很多初版实现只检查共线+坐标范围,却忘了线段可能是竖直或水平的,导致 x 或 y 范围判断失效。务必用点积或参数 t ∈ [0,1] 的方式验证。
C++免费学习笔记(深入):立即使用
在学习笔记中,你将探索 C++ 的入门与实战技巧!











