给定两条线段(p1, q1)和(p2, q2),判断给定的线段是否相交。
在讨论解决方案之前,让我们先定义方向的概念。平面中有序点三元组的方向可以是
–逆时针
–顺时针
–共线
下图显示了(a,b,c) 的不同可能方向
方向在这里有什么用处?
两条线段(p1,q1)和(p2,q2)相交,当且仅当以下两个条件之一得到验证时
1. 一般情况:
– ( p1 , q1 , p2 ) 和 ( p1 , q1 , q2 ) 具有不同的方向,并且
– ( p2 , q2 , p1 ) 和 ( p2 , q2 , q1 ) 具有不同的方向。
例子:
2.特殊情况
– ( p1 , q1 , p2 )、( p1 , q1 , q2 )、( p2 , q2 , p1 ) 和 ( p2 , q2 , q1 ) 均为共线,且 – ( p1 , q1 ) 和 ( p2 , q2 )
的 x 投影相交 – ( p1 , q1 ) 和 ( p2 , q2 )的 y 投影相交
例子:
以下是基于上述想法的实现:
// C# program to check if two given line segments intersect
using System;
using System.Collections.Generic;
class GFG
{
public class Point
{
public int x;
public int y;
public Point(int x, int y)
{
this.x = x;
this.y = y;
}
};
// Given three collinear points p, q, r, the function checks if
// point q lies on line segment 'pr'
static Boolean onSegment(Point p, Point q, Point r)
{
if (q.x <= Math.Max(p.x, r.x) && q.x >= Math.Min(p.x, r.x) &&
q.y <= Math.Max(p.y, r.y) && q.y >= Math.Min(p.y, r.y))
return true;
return false;
}
// To find orientation of ordered triplet (p, q, r).
// The function returns following values
// 0 --> p, q and r are collinear
// 1 --> Clockwise
// 2 --> Counterclockwise
static int orientation(Point p, Point q, Point r)
{
// See https://www.geeksforgeeks.org/orientation-3-ordered-points/
// for details of below formula.
int val = (q.y - p.y) * (r.x - q.x) -
(q.x - p.x) * (r.y - q.y);
if (val == 0) return 0; // collinear
return (val > 0)? 1: 2; // clock or counterclock wise
}
// The main function that returns true if line segment 'p1q1'
// and 'p2q2' intersect.
static Boolean doIntersect(Point p1, Point q1, Point p2, Point q2)
{
// Find the four orientations needed for general and
// special cases
int o1 = orientation(p1, q1, p2);
int o2 = orientation(p1, q1, q2);
int o3 = orientation(p2, q2, p1);
int o4 = orientation(p2, q2, q1);
// General case
if (o1 != o2 && o3 != o4)
return true;
// Special Cases
// p1, q1 and p2 are collinear and p2 lies on segment p1q1
if (o1 == 0 && onSegment(p1, p2, q1)) return true;
// p1, q1 and q2 are collinear and q2 lies on segment p1q1
if (o2 == 0 && onSegment(p1, q2, q1)) return true;
// p2, q2 and p1 are collinear and p1 lies on segment p2q2
if (o3 == 0 && onSegment(p2, p1, q2)) return true;
// p2, q2 and q1 are collinear and q1 lies on segment p2q2
if (o4 == 0 && onSegment(p2, q1, q2)) return true;
return false; // Doesn't fall in any of the above cases
}
// Driver code
public static void Main(String[] args)
{
Point p1 = new Point(1, 1);
Point q1 = new Point(10, 1);
Point p2 = new Point(1, 2);
Point q2 = new Point(10, 2);
if(doIntersect(p1, q1, p2, q2))
Console.WriteLine("Yes");
else
Console.WriteLine("No");
p1 = new Point(10, 1); q1 = new Point(0, 10);
p2 = new Point(0, 0); q2 = new Point(10, 10);
if(doIntersect(p1, q1, p2, q2))
Console.WriteLine("Yes");
else
Console.WriteLine("No");
p1 = new Point(-5, -5); q1 = new Point(0, 0);
p2 = new Point(1, 1); q2 = new Point(10, 10);;
if(doIntersect(p1, q1, p2, q2))
Console.WriteLine("Yes");
else
Console.WriteLine("No");
}
}
/* This code contributed by PrinciRaj1992 */
输出:
No
Yes
No
时间复杂度: O(1)
空间复杂度:O(1)
资料来源:
http://www.dcs.gla.ac.uk/~pat/52233/slides/Geometry1x1.pdf
《算法导论》第三版,作者:Clifford Stein、Thomas H. Cormen、Charles E. Leiserson、Ronald L. Rivest