描述:
线段相交有两种情形:一种是“规范相交”,另一种是“非规范相交”。规范相交是指两条线段恰有唯一一个不是端点的公共点。即如果一条线段的端点在另一条线段上则不视为相交。如果两条线段有部分重合,也不视为相交。而非规范相交则把以上两种情况都视为相交。如下图所示:
规范相交认为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与点c为顶点
以qa,qb,qc,qd来确定两两直线的关系(bc与cd,ac与cd,ab与ad,ab与ac)
如果 qa与qb(或者qc与qd)乘积大于0,代表关于c(或关于a)两条线的叉积呈散开状,也就没有交点;
如果等于0 则有一点在线上,如果小于0,则相交;
于是乎相交的问题便全部解决了;
那么如何求出两直线的交点呢?
以上内容摘自《算法艺术与信息学竞赛》。
思路利用叉积求得点P分线段DC的比,然后利定比分点坐标公式求得分点P的坐标。
既然全部解决了,便是写代码问题了
#include<bits/stdc++.h>
using namespace std;
struct dian
{
double x;
double y;
};
double cross(dian a,dian b,dian c)
{
return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}
int kspc(dian a,dian b,dian c,dian d)
{
if(max(a.x,b.x)<min(c.x,d.x))return 1;
if(min(a.x,b.x)>max(c.x,d.x))return 1;
if(max(a.y,b.y)<min(c.y,d.y))return 1;
if(min(a.y,b.y)>max(c.y,d.y))return 1;
return 0;
}
void solve()
{
dian a,b,c,d;
scanf("%lf %lf %lf %lf %lf %lf %lf %lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y,&d.x,&d.y);
if(kspc(a,b,c,d)==1)
{
cout<<"NO"<<endl;
}
else
{
double qa=cross(c,d,b);
double qb=cross(c,d,a);
double qc=cross(a,b,d);
double qd=cross(a,b,c);
if(qd*qc>0||qb*qa>0)
{
cout<<"NO"<<endl;
}
else if(qd*qc==0||qb*qa==0)
{
cout<<"YES"<<endl;
}
else
{
double s1=fabs(qd)/fabs(qc);
printf("YES (%.3f,%.3f)\n",(c.x+s1*d.x)/(1+s1),(c.y+s1*d.y)/(1+s1));
}
}
}
int main()
{
int t;
cin>>t;
while(t--)
{
solve();
}
}