TZOJ3348--线段相交Ⅲ

描述

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

规范相交认为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叉积为x,则有x>0,a到b方向为逆时针,x==0,la,b向量共线,x<0，a到b为顺时针方向，所以线段两端点若在另一线段两侧则叉积为一正一负

线段交点求法：设p1,p2为一条线段两端点，p3,p4为另一条线段两端点，两线段相交且交点为x，我们可以先求出p1p2和p1p3所在平行四边形的高与p2p4和p1p2所在平行四边形的高的比值h，此时，通过p3x=p1p3*h+ox就可以求出交点x,其中o为原点

代码

  

#include <bits/stdc++.h>
using namespace std;
struct point
{
	double x , y;
}a,b,c,d;
double cross(point p1,point p2,point p3)
{
	return (p2.x-p1.x)*(p3.y-p1.y)-(p2.y-p1.y)*(p3.x-p1.x);
}
void solve()
{
    if(min(a.x,b.x)>max(c.x,d.x)||min(a.y,b.y)>max(c.y,d.y)||min(c.x,d.x)>max(a.x,b.x)||min(c.y,d.y)>max(a.y,c.y))
    {
    	cout<<"NO"<<endl;
    	return;
	}
    double e = cross(a,b,c);
    double f = cross(a,b,d);
    double g = cross(c,d,a);
    double k = cross(c,d,b);
    if(e*f>0||g*k>0)
    {
    	cout<<"NO"<<endl;
    	return;
	}
	if(e*f==0||g*k==0)
	{
		cout<<"YES"<<endl;
		return;
	}
	double d1 = fabs(e);
	double d2 = fabs(f);
	double h = d1/(d1+d2);
	double x = (d.x-c.x)*h+c.x;
	double y = (d.y-c.y)*h+c.y;
	printf("YES (%.3f,%.3f)\n",x,y);
}
int main()
{
	int t;
	cin>>t;
	while(t--)
	{
		cin>>a.x>>a.y>>b.x>>b.y>>c.x>>c.y>>d.x>>d.y;
		solve();
	}
}