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

本体考察线段是否相交，关于线段相交，我们非常容易想到快速排斥；

那么最开始的门槛过了，接下来便是求出两条线段是怎么一个情况

接下来便是求出两条线段的方向问题；

经过同学指点，最后得知这个可以运用叉积来判断；

叉积的原理在我之前的博客中有详细描述过了，这里就简略跳过；

两条直线，分别以点a与点c为顶点

以qa，qb，qc，qd来确定两两直线的关系（bc与cd，ac与cd，ab与ad，ab与ac）

如果 qa与qb（或者qc与qd）乘积大于0，代表关于c（或关于a）两条线的叉积呈散开状，也就没有交点；

如果等于0 则有一点在线上，如果小于0，则相交；

于是乎相交的问题便全部解决了；

那么如何求出两直线的交点呢？

 

以上内容摘自《算法艺术与信息学竞赛》。

思路利用叉积求得点P分线段DC的比，然后利定比分点坐标公式求得分点P的坐标。

既然全部解决了，便是写代码问题了

#include<bits/stdc++.h>
using namespace std;
struct dian
{
    double x;
	double y;
};
double cross(dian a,dian b,dian c)
{
	return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}
int kspc(dian a,dian b,dian c,dian d)
{
	if(max(a.x,b.x)<min(c.x,d.x))return 1;
	if(min(a.x,b.x)>max(c.x,d.x))return 1;
	if(max(a.y,b.y)<min(c.y,d.y))return 1;
	if(min(a.y,b.y)>max(c.y,d.y))return 1;
	return 0;
}
void solve()
{
		dian 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(kspc(a,b,c,d)==1)
	    {
	        cout<<"NO"<<endl;
	    }
	    else
		{
			double qa=cross(c,d,b);
			double qb=cross(c,d,a);
			double qc=cross(a,b,d);
			double qd=cross(a,b,c);
			if(qd*qc>0||qb*qa>0)
			{
				cout<<"NO"<<endl; 
			}
			else if(qd*qc==0||qb*qa==0)
			{
				cout<<"YES"<<endl;
			}
			else
			{
				double s1=fabs(qd)/fabs(qc);
				printf("YES (%.3f,%.3f)\n",(c.x+s1*d.x)/(1+s1),(c.y+s1*d.y)/(1+s1));
			}
    	}
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        solve();
	} 
		
	
}