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.判断两条直线是否相交

快速排斥实验

 cd中y的最大值也要比ab中y的最小值要小，所以不可能相交；

ef中x的最小值也要比ab中x的最大值要大，所以不可能相交；

反之同理； 

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))
    {
        cout<<"NO\n";
        return ;
    }

 跨立实验

我们利用叉积来判断线段的方向

double cross(ss a,ss b,ss c)
{
	return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}

 如果CA*CB<0,说明A点与B点在直线CD的两侧；

 如果CA*CB==0，说明A点与B点有一个点在CD的直线上；

 如果 CA*CB>0，说明A点与B点在CD的同侧，也就不可能相交；

如果没有做快速排斥实验，有可能出现下面的这种情况，明明A点与B点不在C的同侧，仍旧不相交。快速排斥实验是必不可少的。

 2.求交点

不得不说，求交点是一个看起来非常简单的东西（然鹅我遇到了发现我不会）；

在已知四个点的前提下，求交点最方便的就是利用定比分点直接求；

定比分点

设坐标轴上一有向线段的起点和终点的坐标分别为x1和x2，分点M分此有向线段的比为λ，那么，分点M的坐标x=(x1+λx2)/(1+λ);

 

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
struct ss
{
    double x,y;
};
double cross(ss a,ss b,ss c)
{
	return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}//叉积
void solve()
{
	ss a,b,c,d;
    scanf("%lf %lf %lf %lf %lf %lf %lf %lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y,&d.x,&d.y);
	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))
    {
        cout<<"NO\n";
        return ;
    }//快速排斥
	double x=cross(a,b,c);
	double g=cross(a,b,d);
	double y=cross(c,d,a);
	double k=cross(c,d,b);
	if(x*g>0||y*k>0)
	{
		printf("NO\n");
        return ;
	}//跨立
	if(x*g==0||y*k==0)
	{
		printf("YES\n");
		return ;
	}
		double s1=fabs(x)/fabs(g);
		printf("YES (%.3lf,%.3lf)\n",(c.x+s1*d.x)/(1+s1),(c.y+s1*d.y)/(1+s1));//计算交点

}
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        solve();
    }
}