TZOJ 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

判断线段是否相交有两个步骤：

1、快速排斥实验

先判断两条线段在x方向和y方向上是否都有投影重合，只有这样两条线段才有可能相交

2、跨立实验

判断一条线段的两端点在另一条线段的两边，用到了叉乘

 

 需要判断两次，否则会出现这种情况

关于非规范相交：特判叉乘等于零的情况

#include<bits/stdc++.h>
using namespace std;
#define x first
#define y second
typedef pair<double,double>PII;
PII p[7];
double cross(PII a,PII b,PII c)
{
    PII a1,b1;
    a1={a.x-b.x,a.y-b.y};
    b1={a.x-c.x,a.y-c.y};
    return a1.x*b1.y-a1.y*b1.x;
}//叉乘
bool reject(PII a,PII b,PII c,PII d)
{
    if(max(a.x,b.x)>=min(c.x,d.x)&&min(a.x,b.x)<=max(c.x,d.x)&&
    max(a.y,b.y)>=min(c.y,d.y)&&min(a.y,b.y)<=max(c.y,d.y))
    return true;
    else return false;
}//快速排斥实验
void solve()
{
    if(reject(p[1],p[2],p[3],p[4]))
    {
        double t1,t2,t3,t4;
        t1=cross(p[1],p[2],p[3]);
        t2=cross(p[1],p[2],p[4]);
        t3=cross(p[3],p[4],p[1]);
        t4=cross(p[3],p[4],p[2]);
        
        if(t1*t2<0&&t3*t4<0)//叉乘结果异号，说明此点在另一条线段两点之间
        {
            double t=fabs(t1)/(fabs(t1)+fabs(t2));//a X b==|a||b|sin,平行四边形，面积之比等于两高之比
            double a,b;
            a=t*(p[4].x-p[3].x)+p[3].x;//a-p[3].x=h1/(h1+h2)*(p[4].x-p[3].x) 
            b=t*(p[4].y-p[3].y)+p[3].y;//这里的式子和t有关右边是p[4]-p[3]
            printf("YES (%.3lf,%.3lf)\n",a,b);
        }
        else if(t1*t2==0||t3*t4==0)
        {
            cout<<"YES"<<endl;
        }
        else 
        {
            cout<<"NO"<<endl;
        }
    }
    else
    {
        cout<<"NO"<<endl;
    }
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        for(int i=1;i<=4;i++)cin>>p[i].x>>p[i].y;
        solve();
    }

    return 0;
}