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

题目大意

判断线段是否相交，是则判断是否规范相交。不相交输出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;
    }
}