计算几何之线段相交（51nod1264）

设两线段端点P（p1, p2)， Q（q1, q2），分两步确定两线段是否相交

1. 快速排斥实验
   设以线段 P1P2 为对角线的矩形为 R，设以线段 Q1Q2 为对角线的矩形为 T，若 R、T 不相交，则两线段不可能相交

   假设 P1 = (x1, y1), P2 = (x2, y2), Q1 = (x3, y3), Q2 = (x4, y4)，设矩形 R 的 x 坐标的最小边界为 minRX = min(x1, x2)，以此类推，将矩形表示为 R = (minRX, minRY, maxRX, maxRY) 的形式，

   若两矩形相交，则相交的部分构成了一个新的矩形 F，如图所示，我们可以知道 F 的 minFX = max(minRX, minTX), minFY = max(minRY, minTY), maxFX = min(maxRX, maxTX), maxFY = min(maxRY, maxTX)，

   得到 F 的各个值之后，只要判断矩形 F 是否成立就知道 R 和 T 到底有没有相交了，若 minFX > maxFX 或 minFY > maxFy 则 F 无法构成，RT不相交，否则 RT相交。具体如图所示：
   
   2.跨立实验：
   如果两线段相交，则两线段必然相互跨立对方。

   若P1P2跨立Q1Q2 ，则矢量 ( P1 - Q1 ) 和( P2 - Q1 )位于矢量( Q2 - Q1 ) 的两侧，即( P1 - Q1 ) × ( Q2 - Q1 ) * ( P2 - Q1 ) × ( Q2 - Q1 ) < 0。上式可改写成( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) > 0。当 ( P1 - Q1 ) × ( Q2 - Q1 ) = 0 时，说明 ( P1 - Q1 ) 和 ( Q2 - Q1 )共线，但是因为已经通过快速排斥试验，所以 P1 一定在线段 Q1Q2上；同理，( Q2 - Q1 ) ×(P2 - Q1 ) = 0 说明 P2 一定在线段 Q1Q2上。

   所以判断P1P2跨立Q1Q2的依据是：( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) >= 0。同理判断Q1Q2跨立P1P2的依据是：( Q1 - P1 ) × ( P2 - P1 ) * ( P2 - P1 ) × ( Q2 - P1 ) >= 0。具体情况如下图所示
   

   3. 满足快速排斥实验再满足跨立实验即可证明相交
      

代码：例（51nod1264）

#include <bits/stdc++.h>
#pragma GCC diagnostic error "-std=c++11"
#define d(x) cout << (x) << endl
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int mod = 1e9 + 7;
const int N = 14;

struct point	//端点
{
    double x, y;
};
struct linseg	//线段
{
    point a, b;
    linseg(){
        scanf("%lf%lf%lf%lf", &a.x, &a.y, &b.x, &b.y);
    }
};
double xmulti(point p0, point p1, point p2)     //p0p1 与 p0p2 叉乘
{
    return ((p1.x - p0.x) * (p2.y - p0.y) - (p2.x - p0.x) * (p1.y - p0.y));
}
bool solve(linseg p, linseg q)
{
    if(max(p.a.x, p.b.x) >= min(q.a.x, q.b.x)&&     //快速排斥实验
        max(q.a.x, q.b.x) >= min(p.a.x, p.b.x)&&
        max(p.a.y, p.b.y) >= min(q.a.y, q.b.y)&&
        max(q.a.y, q.b.y) >= min(p.a.y, p.b.y)&&
        xmulti(q.a, p.a, q.b) * xmulti(q.a, p.b, q.b) <= 0 &&   //跨立实验
        xmulti(p.a, q.a, p.b) * xmulti(p.a, q.b, p.b) <= 0 
    ){  return true;
    }else{
        return false;
    }
}
int main()
{
    int t;
    for (cin >> t; t; t--){
        linseg p = linseg();
        linseg q = linseg();
        if(solve(p, q)){
            d("Yes");
        }else{
            d("No");
        }
    }
    return 0;
}