题目
描述
线段相交有两种情形:一种是“规范相交”,另一种是“非规范相交”。规范相交是指两条线段恰有唯一一个不是端点的公共点。即如果一条线段的端点在另一条线段上则不视为相交。如果两条线段有部分重合,也不视为相交。而非规范相交则把以上两种情况都视为相交。如下图所示:
规范相交认为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
题目大意
判断线段是否相交,是则判断是否规范相交。不相交输出NO,相交输出YES,如果是规范相交再输出交点坐标。
思路与算法
首先进行快速排斥实验,符合在进行跨立实验,也符合则线段相交。再利用叉乘判断是否存在一条线段的端点在另一条线段上,有则为非规范相交。
交点坐标
(=(xc-xa)*(yd-yc)-(xd-xc)*(yc-ya))/((xb-xa)*(yd-yc)-(xd-xc)*(yb-ya))
代码
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
double xa,ya,xb,yb,xc,yc,xd,yd;
int t;
double cross(double x1,double y1,double x2,double y2,double x0,double y0){
return (x1-x0)*(y2-y0)-(x2-x0)*(y1-y0);
}
int fastjudge(double xa,double ya,double xb,double yb,double xc,double yc,double xd,double yd){
if(max(xc,xd)<min(xa,xb) || max(yc,yd)<min(ya,yb) ||max(xa,xb)<min(xc,xd) || max(ya,yb)<min(yc,yd))
return 0;
return 1;
}
int judge(double xa,double ya,double xb,double yb,double xc,double yc,double xd,double yd){
if(fastjudge(xa,ya,xb,yb,xc,yc,xd,yd)&&cross(xb,yb,xc,yc,xa,ya)* cross(xb,yb,xd,yd,xa,ya)<=0&& cross(xc,yc,xa,ya,xd,yd)*cross(xc,yc,xb,yb,xd,yd)<=0)
return 1;
return 0;
}
int judgegui(double x,double y,double x1,double y1,double x2,double y2){
if(cross(x,y,x1,y1,x2,y2)==0)
return 1;
return 0;
}
int main(){
cin>>t;
while(t--){
cin>>xa>>ya>>xb>>yb>>xc>>yc>>xd>>yd;
if(judge(xa,ya,xb,yb,xc,yc,xd,yd)){
if(judgegui(xa,ya,xc,yc,xd,yd)||judgegui(xb,yb,xc,yc,xd,yd)||judgegui(xc,yc,xa,ya,xb,yb)||judgegui(xd,yd,xa,ya,xb,yb))
cout<<"YES"<<endl;
else{
double r=((xc-xa)*(yd-yc)-(xd-xc)*(yc-ya))/((xb-xa)*(yd-yc)-(xd-xc)*(yb-ya));
double x=xa+r*(xb-xa),y=ya+r*(yb-ya);
cout<<fixed<<setprecision(3)<<"YES"<<" ("<<x<<","<<y<<")"<<endl;
}
}
else
cout<<"NO"<<endl;
}
}