tzoj 3348: 线段相交Ⅲ

题目：

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

规范相交认为a，b两种情况都是不相交的，而非规范相交认为a，b两种情况都是相交的。

本题要求判断两条线段是否相交。如果是规范相交则输出YES，并输出交点坐标，如果是非规范相交则只需输出YES，如果不相交则输出NO。

输入：

输入有多组数据，T表示输入数据的组数。每组测试数据有两行第一行输入一条线段的两个端点的坐标，第二行输入另一个线段的两个端点的坐标。

输出：

输入有多组数据，T表示输入数据的组数。每组测试数据有两行第一行输入一条线段的两个端点的坐标，第二行输入另一个线段的两个端点的坐标。

样例输入：

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

题意：

有非规范相交、规范相交、以及不相交三种情况，非规范相交就是图片中给的两种情况，其他的相交   都是规范相交，然后让我们进行判断

思路：

 首先我们可以把不相交筛出来，通过使用快速排斥实验+跨立实验  两种方法可以两条直线是否不相交。

快速排斥实验：

我们先假设第一条线段两端为p1,p2   第二条线段为t1,t2，以第一条线段为对角线可以做出一个矩形P，以第二条为对角线可以做出另一个矩形T，如果两个矩形不想交，那么两条直线必定不相交。

bool fir_reject(point p1,point p2,point t1,point t2)
{
	if(max(p1.x,p2.x)<min(t1.x,t2.x)||
	   max(p1.y,p2.y)<min(t1.y,t2.y)||
	   max(t1.x,t2.x)<min(p1.x,p2.x)||
	   max(t1.y,t2.y)<min(p1.y,p2.y) 
	)
	return true;//说明不相交；
	return false; 
}

在第一个实验中无法筛去，某种不相交的情况，所以要再深入一步判断。

 跨立实验：

 

先看第一张图，当两直线相交时，p1,p2两点一定会位于t1t2直线的两端，我们可以通过叉积来判断p1,p2是否位于该直线两端，只有当  和的结果相同，说明是在两侧的，但看这一个条件是不够的。

我们可以看一下图二，这种情况是不是也成立，所以我们还要再判，t1,t2两点是否位于p1p2直线两端，判断方法和上述一致。

	//跨立实验
	double a=chaji(t2,p1,t2,t1);
	double b=chaji(t2,t1,t2,p2);
	double c=chaji(p1,t1,p1,p2);
	double d=chaji(p1,p2,p1,t2);
	if(a*b<0||c*d<0)
	{
		cout<<"NO"<<endl;
		return;
	} 
	if(a*b==0||c*d==0)
	{
		cout<<"YES"<<endl;
		return;
		//说明不规范相交有共线; 
	}

不规范相交看图可知大概就是有共线的情况，所以就是当我们在跨立实验时，有叉积为0的情况就是不规范相交。

最后如何求出规范相交的x,y点呢？

暂且只看懂了这一种办法 哎

 因为叉积    ，所以可以直接把三角形面积求出来。

通过代码：

#include<bits/stdc++.h>
#define IOS ios::sync_with_stdio(false);cin.tie(NULL);cout.tie(NULL)
#define x first
#define y second
using namespace std;
typedef long long ll;
const int N=10;
struct point
{
	double x,y;
}; 
int n,m;
bool fir_reject(point p1,point p2,point t1,point t2)
{
	if(max(p1.x,p2.x)<min(t1.x,t2.x)||
	   max(p1.y,p2.y)<min(t1.y,t2.y)||
	   max(t1.x,t2.x)<min(p1.x,p2.x)||
	   max(t1.y,t2.y)<min(p1.y,p2.y) 
	)
	return true;//说明不相交；
	return false; 
}
double chaji(point a,point b,point c,point d)
{
	point aa,bb;
	aa.x=a.x-b.x,aa.y=a.y-b.y; 
	bb.x=c.x-d.x,bb.y=c.y-d.y;
	return aa.x*bb.y-aa.y*bb.x; 
}
void solve()
{
	point p1,p2,t1,t2;
	cin>>p1.x>>p1.y>>p2.x>>p2.y>>t1.x>>t1.y>>t2.x>>t2.y;
	if(fir_reject(p1,p2,t1,t2)) 
	{
		cout<<"NO"<<endl;
		return; 
	}
	//跨立实验
	double a=chaji(t2,p1,t2,t1);
	double b=chaji(t2,t1,t2,p2);
	double c=chaji(p1,t1,p1,p2);
	double d=chaji(p1,p2,p1,t2);
	if(a*b<0||c*d<0)
	{
		cout<<"NO"<<endl;
		return;
	} 
	if(a*b==0||c*d==0)
	{
		cout<<"YES"<<endl;
		return;
		//说明不规范相交有共线; 
	}
	//否则输出相交坐标；
	double t=fabs(chaji(t1,p1,t1,p2))/(fabs(chaji(t2,p2,t2,p1))+fabs(chaji(t1,p1,t1,p2))); 
	double x=t*(t2.x-t1.x)+t1.x;
	double y=t*(t2.y-t1.y)+t1.y;
	printf("YES (%.3f,%.3f)\n",x,y);
}
int main()
{
	//IOS;
	int T;
	T=1;
	cin>>T;
	while(T--)
	{
		solve();
	}
}

睡觉了，困~