3348: 线段相交Ⅲ

描述

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

规范相交认为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

题解：

这题主要是学习向量与向量，向量与点之间的关系，涉及叉积等知识

#include<iostream> 
#include<algorithm> 
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<iomanip>
#include<string.h>
#include<set>
#include<string>
#include<map>
#include<bitset>
#pragma warning(disable : 4996)
using namespace std;
//定义点
class point
{
public:
    double x;
    double y;
//两个构造
    point()
    {
    }
    point(double x, double y)
    {
        this->x = x;
        this->y = y;
    }
};
//点之间的加减乘除，重构运算符
point operator-(point a, point b)
{
    return point(a.x - b.x, a.y - b.y);
}
point operator+(point a, point b)
{
    return point(a.x + b.x, a.y + b.y);
}
point operator/(point a, int x)
{
    return point(a.x / x, a.y / x);
}
point operator*(point a, double b)
{
    return point(a.x * b, a.y * b);
}
//这个特殊一点，实际上这是个叉积式子，及两个向量的叉积，x1*y2-y1*x2
double operator*(point a, point b)
{
    return a.x * b.y - a.y * b.x;
}
//也是叉积式子,不过三个变量都是点，上面的叉积式两个变量是向量
double cross(point a, point b, point c)
{
    return (b - a) * (c - a);
}
//求两点之间的长度
double len(point a, point b)
{
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}
//取得交点
point getNode(point a, point u, point b, point v)
{
    
    double t = (a - b) * v / (v * u);
    return a + u * t;
}
int main()
{
    //加快cin和cout的速度
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int k;
    cin >> k;
    while (k--)
    {
        //四个点
        point p[4];
        for (int i = 0; i < 4; i++)
        {
            cin >> p[i].x >> p[i].y;
        }
        //求叉积*叉积，设一个向量a,如果向量b在其逆时针方向，叉积为正，反之为负，同线为0，可借此判断两点相对于直线的位置
        double k1 = cross(p[0], p[1], p[2]) * cross(p[0], p[1], p[3]);
        double k2 = cross(p[2], p[3], p[0]) * cross(p[2], p[3], p[1]);
        //只要一个大于0，就不相交
        if (k1 > 0 || k2 > 0)
        {
            cout << "NO" << endl;
        }
        //如果两者小于零，规范相交
        else if (k1 < 0 && k2 < 0)
        {
            //运用比例求交点，算初中数学？这里的比例可以用叉积来表示，叉积的几何意义是两向量组合成的平行四边形面积，可以找到公共边算比例
            point k = getNode(p[0], p[1] - p[0], p[2], p[3] - p[2]);
            cout << fixed << setprecision(3) << "YES " << "(" << k.x << "," << k.y << ")" << endl;
        }
        //如果都等于0，则在同一直线上
        else if (k1 == 0 && k2 == 0)
        {
            //这里借用两线段的半径和两线段中心点之间的距离判断是否相交
            point p1 = (p[0] + p[1]) / 2;
            point p2 = (p[2] + p[3]) / 2;
            double L = len(p1, p2), R1 = len(p1, p[0]), R2 = len(p2, p[2]);
            if (R1 + R2 < L)
            {
                cout << "NO" << endl;
            }
            else
            {
                cout << "YES" << endl;
            }
        }
        //如果一个等于0，一个小于0，则是T字形相交
        else
        {
            cout << "YES" << endl;
        }
    }
        
    
    return 0;
}