描述
线段相交有两种情形:一种是“规范相交”,另一种是“非规范相交”。规范相交是指两条线段恰有唯一一个不是端点的公共点。即如果一条线段的端点在另一条线段上则不视为相交。如果两条线段有部分重合,也不视为相交。而非规范相交则把以上两种情况都视为相交。如下图所示:
规范相交认为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
题目分析
两个线段规范相交指两个线段的交点既不是两个线段的端点也不在两个线段的延长线上。判断两个线段是否规范相交可以按照下列步骤进行:
1. 计算两个线段的斜率。如果其中一个线段的斜率为无穷大,则代表这个线段是竖直的,而另一个线段不是竖直的,它们一定相交,因为无论垂直于哪一条线段,都可以看做是两条斜率不同的线段之间的相交。
2. 判断两条线段是否平行。如果它们平行,则它们不可能相交。
3. 判断两条线段是否在同一直线上。如果它们在同一直线上,则它们可以相交,但需要进一步判断它们是否重叠或有一个线段在另一个线段的延长线上。
4. 如果两个线段不平行且不在同一直线上,则它们相交,计算它们的交点。如果交点既不在两个线段的端点也不在两个线段的延长线上,则它们是规范相交。否则,它们不是规范相交。
#include <bits/stdc++.h>
using namespace std;
typedef double db;
struct Node{
db x,y;
};
int check(Node t, Node a1, Node a2, Node b1, Node b2){
if (t.x > min(b1.x, b2.x) && t.x < max(b1.x, b2.x) && t.y > min(b1.y, b2.y) && t.y < max(b1.y, b2.y) && t.x > min(a1.x, a2.x) && t.x < max(a1.x, a2.x) && t.y > min(a1.y, a2.y) && t.y < max(a1.y, a2.y)) return 2;
else if (t.x >= min(b1.x, b2.x) && t.x <= max(b1.x, b2.x) && t.y >= min(b1.y, b2.y) && t.y <= max(b1.y, b2.y) && t.x >= min(a1.x, a2.x) && t.x <= max(a1.x, a2.x) && t.y >= min(a1.y, a2.y) && t.y <= max(a1.y, a2.y)) return 1;
else return 0;
}
void solve()
{
Node a1, a2, b1, b2;
cin >> a1.x >> a1.y >> a2.x >> a2.y >> b1.x >> b1.y >> b2.x >> b2.y;
if ((a1.y - a2.y) * (b1.x - b2.x) == (a1.x - a2.x) * (b1.y - b2.y)){
if (a1.x - a2.x == 0) {
if (a1.x == b1.x) {
if (check(a1, a1, a2, b1, b2) || check(a2, a1, a2, b1, b2)) cout<<"YES\n";
else cout<<"NO\n";
}
else cout<<"NO\n";
}
else {
db num1 = a1.y - (a1.y - a2.y) / (a1.x - a2.x) * a1.x,num2 = b1.y - (b1.y - b2.y) / (b1.x - b2.x) * b1.x;
if (num1 == num2) {
if (check(a1, a1, a2, b1, b2) || check(a2, a1, a2, b1, b2)) cout << "YES" << endl;
else cout<<"NO\n";
}
else cout<<"NO\n";
}
}
else {
if (a1.x == a2.x){
db num = b1.y - (b1.y - b2.y) / (b1.x - b2.x) * b1.x;
Node t = {a1.x, (b1.y - b2.y) / (b1.x - b2.x) * a1.x + num};
if (check(t, a1, a2, b1, b2) == 2) printf("YES (%.3f,%.3f)\n", t.x, t.y);
else if (check(t, a1, a2, b1, b2) == 1) cout<<"YES\n";
else cout<<"NO\n";
}
else if (b1.x == b2.x){
db num = a1.y - (a1.y - a2.y) / (a1.x - a2.x) * a1.x;
Node t = {b1.x, (a1.y - a2.y) / (a1.x - a2.x) * b1.x + num};
if (check(t, a1, a2, b1, b2) == 2) printf("YES (%.3f,%.3f)\n", t.x, t.y);
else if (check(t, a1, a2, b1, b2) == 1) cout<<"YES\n";
else cout<<"NO\n";
}
else{
db num1 = a1.y - (a1.y - a2.y) / (a1.x - a2.x) * a1.x,num2 = b1.y - (b1.y - b2.y) / (b1.x - b2.x) * b1.x;
db k1 = (a1.y - a2.y) / (a1.x - a2.x);
db k2 = (b1.y - b2.y) / (b1.x - b2.x);
Node t;
t.x = (num2 - num1) / (k1 - k2);
t.y = k1 * t.x + num1;
if (check(t, a1, a2, b1, b2) == 2) printf("YES (%.3f,%.3f)\n", t.x, t.y);
else if (check(t, a1, a2, b1, b2) == 1) cout<<"YES\n";
else cout<<"NO\n";
}
}
}
signed main()
{
int _=1;
cin>>_;
while(_--) solve();
}