描述
线段相交有两种情形:一种是“规范相交”,另一种是“非规范相交”。规范相交是指两条线段恰有唯一一个不是端点的公共点。即如果一条线段的端点在另一条线段上则不视为相交。如果两条线段有部分重合,也不视为相交。而非规范相交则把以上两种情况都视为相交。如下图所示:
规范相交认为a,b两种情况都是不相交的,而非规范相交认为a,b两种情况都是相交的。
本题要求判断两条线段是否相交。如果是规范相交则输出YES,并输出交点坐标,如果是非规范相交则只需输出YES,如果不相交则输出NO。
输入
输入有多组数据,T表示输入数据的组数。每组测试数据有两行第一行输入一条线段的两个端点的坐标,第二行输入另一个线段的两个端点的坐标。
输出
对于每组测试数据,输出一行。如果是规范相交则输出YES,并输出交点坐标(小数点后面保留3位),如果是非规范相交则只需输出YES,如果不相交则输出NO。
样例输入
4
0 0 1 1
0 1 1 0
0 0 2 2
2 2 3 3
0 0 2 2
1.5 1.5 3 3
0 0 1 1
2 2 3 3
样例输出
YES (0.500,0.500)
YES
YES
NO
题解
使用快速排斥实验进行初步判断是否可能相交,若是则利用跨立实验来进一步确认是规范相交还是非规范相交
快速排斥实验:两线段的投影若不相交则两线段不可能相交
跨立实验:设两向量a,b叉积为x,则有x>0,a到b方向为逆时针,x==0,la,b向量共线,x<0,a到b为顺时针方向,所以线段两端点若在另一线段两侧则叉积为一正一负
线段交点求法:设p1,p2为一条线段两端点,p3,p4为另一条线段两端点,两线段相交且交点为x,我们可以先求出p1p2和p1p3所在平行四边形的高与p2p4和p1p2所在平行四边形的高的比值h,此时,通过p3x=p1p3*h+ox就可以求出交点x,其中o为原点
代码
#include <bits/stdc++.h>
using namespace std;
struct point
{
double x , y;
}a,b,c,d;
double cross(point p1,point p2,point p3)
{
return (p2.x-p1.x)*(p3.y-p1.y)-(p2.y-p1.y)*(p3.x-p1.x);
}
void solve()
{
if(min(a.x,b.x)>max(c.x,d.x)||min(a.y,b.y)>max(c.y,d.y)||min(c.x,d.x)>max(a.x,b.x)||min(c.y,d.y)>max(a.y,c.y))
{
cout<<"NO"<<endl;
return;
}
double e = cross(a,b,c);
double f = cross(a,b,d);
double g = cross(c,d,a);
double k = cross(c,d,b);
if(e*f>0||g*k>0)
{
cout<<"NO"<<endl;
return;
}
if(e*f==0||g*k==0)
{
cout<<"YES"<<endl;
return;
}
double d1 = fabs(e);
double d2 = fabs(f);
double h = d1/(d1+d2);
double x = (d.x-c.x)*h+c.x;
double y = (d.y-c.y)*h+c.y;
printf("YES (%.3f,%.3f)\n",x,y);
}
int main()
{
int t;
cin>>t;
while(t--)
{
cin>>a.x>>a.y>>b.x>>b.y>>c.x>>c.y>>d.x>>d.y;
solve();
}
}