码迷,mamicode.com
首页 > 其他好文 > 详细

【TOJ 1288】计算几何练习题――线段相交(模板题)

时间:2018-05-16 00:46:51      阅读:316      评论:0      收藏:0      [点我收藏+]

标签:条件   测试数据   double   http   name   space   min   style   cin   

描述

线段相交测试在计算几何中是经常用到的,给定线段P1P2(P1和P2是线段的两端点,且不重合)、P3P4(P3和P4是线段的两端点,且不重合),判断P1P2和P3P4是否相交。P1P2和P3P4不重合,即指只存在一个点P,它既落在P1P2上又落在P3P4上(含线段的端点)。

输入

输入数据有多组,第一行为测试数据的组数N,下面包括2N行,每组测试数据含2行,第一行为P1P2的坐标值,第二行为P3P4的坐标值,比如下面的数据
1
0 0 1 1
2 2 3 3
表示P1、P2、P3、P4的坐标分别为:P1(0,0),P2(1,1),P3(2,2),P4(3,3)

输出

判断每组数据中的线段P1P2和P3P4是否相交,如果相交输出YES,否则输出NO。每组数据输出占一行。

样例输入

1
0 0 1 1
2 2 3 3

样例输出

NO

题解

判断两线段是否相交:

① 快速排斥

就是初步的判断一下,两条线段是不是相交,以两条线段为对角线的矩形,如果不重合的话,那么两条线段一定不可能相交。

技术分享图片                

1.线段ab的低点低于cd的最高点(可能重合)

2.cd的最左端小于ab的最右端(可能重合)

3.cd的最低点低于ab的最高点(加上条件1,两线段在竖直方向上重合)

4.ab的最左端小于cd的最右端(加上条件2,两直线在水平方向上重合)

 

综上4个条件,两条线段组成的矩形是重合的。

② 跨立实验

如果两条线段相交,那么必须跨立,就是以一条线段为标准,另一条线段的两端点一定在这条线段的两段

即A、B两点在线段CD的两端,C、D两点在线段AB的两端

这里就用到了向量X乘的知识点,有向量X乘的物理意义知:AB x CD=-CD x AB

看下图:

 技术分享图片

 

(CA x CD)·(CB x CD)<=0 则说明向量CA、CB对于向量CD的方向不同,则A、B在线段CD的两侧,由此可以判断其他点

#include<bits/stdc++.h>
using namespace std;
struct node{
    int x,y;
}a,b,c,d;
bool judge(node a,node b,node c,node d)
{
    //快速排斥 
    if(!(min(a.x,b.x)<=max(c.x,d.x)&&min(c.x,d.x)<=max(a.x,b.x)&&min(a.y,b.y)<=max(c.y,d.y)&&min(c.y,d.y)<=max(a.y,b.y))) 
        return false;

    //跨立实验 
    double u,v,w,z;
    u=(c.x-a.x)*(b.y-a.y)-(b.x-a.x)*(c.y-a.y);  //AC×AB
    v=(d.x-a.x)*(b.y-a.y)-(b.x-a.x)*(d.y-a.y);  //AD×AB
    w=(a.x-c.x)*(d.y-c.y)-(d.x-c.x)*(a.y-c.y);  //CA×AB
    z=(b.x-c.x)*(d.y-c.y)-(d.x-c.x)*(b.y-c.y);  //CB×AB
    
    if(u*v<=1e-9&&w*z<=1e-9)
        return true;
    return false; 
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        cin>>a.x>>a.y>>b.x>>b.y>>c.x>>c.y>>d.x>>d.y;
        if(judge(a,b,c,d))
            printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}

 

【TOJ 1288】计算几何练习题――线段相交(模板题)

标签:条件   测试数据   double   http   name   space   min   style   cin   

原文地址:https://www.cnblogs.com/kannyi/p/9043424.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!