线段交 【判断两线段是否有交点】

题目链接

题目描述

给定N个线段。求有交点的线段对数。
保证没有两条线段共线

输入

一行一个整数N，表示线段的个数
第2~N+1行，每行四个实数，x1,y1,x2,y2，表示线段的两个端点(x1,y1)和(x2,y2)

输出

一行一个整数，表示有交点的线段对数。

样例输入

复制样例数据

3
0.00 0.00 1.00 1.00
0.00 1.00 1.00 0.00
0.00 0.00 1.00 0.00

样例输出

3

提示

(0,0)(1,1)和(0,1)(1,0)有交点
(0,0)(1,1)和(0,0)(1,0)有交点
(0,1)(1,0)和(0,0)(1,0)有交点

对于100%的数据，N≤100
点的坐标范围(−10000,10000)

判断两个线段ab和cd是否有交点

一、快速排斥实验

分别以ab 、cd为对角线做矩阵，若这两个矩阵不相交，则这两个线段一定不相交；即使两个矩阵相交，这两个线段也不一定相交。

二、跨立实验

若ab跨立cd，那么a、b一定分别在cd的两边，则有 (d-a)×(d-c)*(d-b)×(d-c)<0;若 (d-a)×(d-c)*(d-b)×(d-c)=0，那么a、b中有一个点在线段cd上；若cd跨立ab，同理。

AC代码

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <set>
using namespace std;
#define io ios::sync_with_stdio(0),cin.tie(0)
#define ms(arr) memset(arr,0,sizeof(arr))
#define mc(a,b) memcpy(a,b,sizeof(b))
#define inf 0x3f3f3f
#define fin freopen("in.txt", "r", stdin)
#define fout freopen("out.txt", "w", stdout)
typedef long long ll;
typedef unsigned long long ULL;
const int mod=1e9+7;
const int N=1e5+7;
struct node
{
    double x;
    double y;
};
int n;
node ps[N],pe[N];
bool is_xiangjiao(node a,node b,node c,node d)
{
     if(max(c.x,d.x)<min(a.x,b.x)||max(a.x,b.x)<min(c.x,d.x)||max(c.y,d.y)<min(a.y,b.y)||max(a.y,b.y)<min(c.y,d.y))
        return false;
     if(((d.x-a.x)*(d.y-c.y)-(d.y-a.y)*(d.x-c.x))*((d.x-b.x)*(d.y-c.y)-(d.y-b.y)*(d.x-c.x))>0.0000000001)
        return false;
     if(((c.x-a.x)*(b.y-a.y)-(c.y-a.y)*(b.x-a.x))*((d.x-a.x)*(b.y-a.y)-(d.y-a.y)*(b.x-a.x))>0.0000000001)
        return false;
     return true;
}
int main()
{
//    fin;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%lf%lf%lf%lf",&ps[i].x,&ps[i].y,&pe[i].x,&pe[i].y);
    }
    int ans=0;
    for(int i=1;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            if(is_xiangjiao(ps[i],pe[i],ps[j],pe[j])) ans++;
        }
    }
    printf("%d\n",ans);
    return 0;
}