叉积法判断点在凸多边形内:遍历逆时针边,计算点对每条边的叉积,若全≥0则在内部或边界,全>0为严格内部;需校验顶点顺序并谨慎处理浮点误差。

用叉积判断点与每条边的相对位置
核心思路是:对凸多边形按顺序(顺时针或逆时针)遍历每条有向边,计算点相对于该边的叉积符号。若所有叉积同号(全 ≥ 0 或全 ≤ 0),点就在内部或边界上;否则在外部。
假设多边形顶点存为 std::vector<:pair double>></:pair>,按逆时针排列,点为 (px, py),边从 v[i] 指向 v[(i+1)%n],则叉积为:
double cross = (v[i+1].first - v[i].first) * (py - v[i].second)
- (v[i+1].second - v[i].second) * (px - v[i].first);
注意:cross == 0 表示点在该边上;严格内部要求所有 cross > 0(逆时针时)。
- 必须保证顶点顺序一致,混用顺/逆时针会导致符号混乱
- 浮点比较慎用
== 0,建议用std::abs(cross) 判共线 - 若多边形是顺时针存储,应统一检查
cross ,而非翻转逻辑
使用 std::all_of 封装成一行判断
现代 C++ 可用 STL 算法快速实现,避免手写循环出错:
bool pointInConvexPolygon(const std::vector<:pair>>& poly,
double px, double py) {
int n = poly.size();
if (n double {
auto& a = poly[i];
auto& b = poly[(i + 1) % n];
return (b.first - a.first) * (py - a.second)
- (b.second - a.second) * (px - a.first);
};
double sgn = cross(0);
return std::all_of(
std::begin(poly), std::end(poly),
[&](const auto& p) { return cross(&p - &poly[0]) * sgn >= 0; }
);
}
</:pair>
这里用首个叉积 sgn 定义“期望符号”,后续每个叉积与其同号或为零即通过。注意:&p - &poly[0] 是安全获取索引的方式,但更推荐用索引循环避免指针算术。
- 不要直接捕获
i进 lambda——循环变量会变化,导致未定义行为 -
sgn本身可能是 0(点恰在首边上),此时要求所有叉积都为 0,即点在整条边上——这仅当多边形退化时发生,通常可接受 - 若需区分内部/边界,把
>= 0拆成两路判断:先查是否所有cross == 0(在边上),再查是否所有cross > 0(严格内部)
警惕浮点误差导致的误判
叉积接近零时,因浮点舍入可能由正变负,使本应在边上的点被判为外部。典型现象是:点明明在顶点或边上,函数返回 false。
组合式C++代码评审方案,融合静态分析、AI推理、多轮迭代评审和C++专项检查,适用于PR审查、增量代码审查、全项目评审和代码质量评分,触发词包括review cpp、cpp代码评审、C++review、代码审查。
解决方式不是盲目加大容差,而是分层处理:
- 先用宽松容差(如
1e-9)判断是否共线:std::abs(cross) - 对共线情况,额外做线段包含判断:检查
px是否在v[i].first和v[i+1].first的范围内(需考虑 x/y 都可能恒定) - 主逻辑中,统一用带容差的符号比较:
cross >= -1e-9替代cross >= 0
不推荐全局用 std::nextafter 调整边界——它增加复杂度且不解决根本问题。
为什么不用射线法或面积法?
射线法(如水平射线交点计数)虽通用,但对凸多边形是杀鸡用牛刀:它要处理边平行、顶点重合、浮点交点精度等一堆边界,时间复杂度 O(n),常数更大;面积法(比较子三角形面积和)涉及多次开方或乘法,易累积误差,且无法自然区分边界。
叉积法专为凸性设计,O(n) 但常数极小,无开方、无除法,只用加减乘,适合高频调用(如碰撞检测、光栅化预判)。
真正容易被忽略的是:**顶点顺序的校验不能省略**。传入一个你认为“逆时针”的多边形,实际可能是顺时针(比如从图像坐标系读入),此时所有符号翻转,结果全错。上线前至少加一句断言:assert(cross(0) * cross(1) > 0 || std::abs(cross(0))
C++免费学习笔记(深入):立即使用
在学习笔记中,你将探索 C++ 的入门与实战技巧!










