6
最近我正在处理一些计算几何,并且试图找到一种检查两个线段是否相交的方法。我认为我可以使用逆时针方向(简称CCW)来确定。这是我到目前为止的代码:用于确定两个线段是否相交的C++程序
struct point { double x, y };
double CCW(point a, point b, point c)
{ return (b.x-a.x)*(c.y-a.y) - (b.y-a.y)*(c.x-a.x); }
int intersect(point a, point b, point c, point d)
{ return (CCW(a,b,c)*CCW(a,b,d)<0 && CCW(c,d,b)*CCW(c,d,a)<0); }
上面的代码工作了我进入测试案例,它是相当可读的,很容易实现。但是在网上搜索之后,我找到了另一种解决细分交叉问题的方法。该代码与我的代码类似,但它还有一些我的实现忽略的语句。下面是代码:
struct line { point s, e; };
int middle(int a, int b, int c) {
int t;
if (a > b) {
t = a;
a = b;
b = t;
}
if (a <= c && c <= b) return 1;
return 0;
}
int intersect(line a, line b) {
if ((CCW(a.s, a.e, b.s) * CCW(a.s, a.e, b.e) < 0) &&
(CCW(b.s, b.e, a.s) * CCW(b.s, b.e, a.e) < 0)) return 1;
if (CCW(a.s, a.e, b.s) == 0 && middle(a.s.x, a.e.x, b.s.x) && middle(a.s.y, a.e.y, b.s.y)) return 1;
if (CCW(a.s, a.e, b.e) == 0 && middle(a.s.x, a.e.x, b.e.x) && middle(a.s.y, a.e.y, b.e.y)) return 1;
if (CCW(b.s, b.e, a.s) == 0 && middle(b.s.x, b.e.x, a.s.x) && middle(b.s.y, b.e.y, a.s.y)) return 1;
if (CCW(b.s, b.e, a.e) == 0 && middle(b.s.x, b.e.x, a.e.x) && middle(b.s.y, b.e.y, a.e.y)) return 1;
return 0;
}
有人能解释一下这两个实现之间的差异,这是使用更安全?提前致谢。