线段相交Ⅲ

线段相交有两种情形:一种是“规范相交”,另一种是“非规范相交”。规范相交是指两条线段恰有唯一一个不是端点的公共点。即如果一条线段的端点在另一条线段上则不视为相交。如果两条线段有部分重合,也不视为相交。而非规范相交则把以上两种情况都视为相交。如下图所示:

规范相交认为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,C,D,线段分别为AB,CD。可以通过计算ABxCD是否为0判断两线段(向量)是否平行。

线段平行的情况判断是否为非规范相交证明平行直线AB,CD非规范相交可以通过证明AB,CD共线且AB,CD在x轴上的投影是否相交。

线段不平行的情况先通过叉乘找交点再判断交点在不在线段上。

#include <bits/stdc++.h>
using namespace std;
class tt
{
public:
    double x, y;
};
tt d[10];
tt ff(tt a, tt b)
{
    return {b.x - a.x, b.y - a.y};
}
double check(tt a, tt b)
{
    return a.x * b.y - a.y * b.x;
}
void solve()
{
    for (int i = 1; i <= 4; i++)
        cin >> d[i].x >> d[i].y;
    if (fabs(check(ff(d[1], d[2]), ff(d[3], d[4]))) < 1e-6)
    {
        if (fabs(check(ff(d[1], d[2]), ff(d[1], d[4]))) < 1e-6 && max(d[1].x, d[3].x) <= min(d[2].x, d[4].x) + 1e-6)
            cout << "YES" << endl;
        else
            cout << "NO" << endl;
    }
    else
    {
        double k1 = check(ff(d[1], d[2]), ff(d[1], d[3])),
               k2 = check(ff(d[1], d[2]), ff(d[1], d[4])),
               k3 = check(ff(d[3], d[1]), ff(d[3], d[4])),
               k4 = check(ff(d[3], d[2]), ff(d[3], d[4]));
        if (fabs(k1 * k2) < 1e-6 || fabs(k3 * k4) < 1e-6) // 点在线段上
        {
            cout << "YES" << endl;
        }
        else if (k1 * k2 < 0 || k3 * k4 < 0) // 两点在线段两侧
        {
            double s1 = check(ff(d[1], d[4]), ff(d[3], d[4])), s2 = fabs(check(ff(d[1], d[2]), ff(d[3], d[4])));
            double k = fabs(s1 / s2);
            tt t = {k * ff(d[1], d[2]).x, k * ff(d[1], d[2]).y};
            printf("YES (%.3f,%.3f)\n", d[1].x + t.x, d[1].y + t.y);
        }
        else
        {
            cout << "NO" << endl;
        }
    }
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int T = 1;
    cin >> T;
    while (T--)
    {
        solve();
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值