TZOJ 3348: 线段相交Ⅲ(叉乘+快速排斥实验+跨立实验)

2023-10-11 19:50

本文主要是介绍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

题目来源

TZOJ

思路:判断线段相交用叉乘+快速排斥

第一步:快速排斥实验

如果两线段在x,y的投影都不重合,是不可能会相交的。

第二步:跨立实验 

两个坐标A(x1,y1),B(x2,y2),那么AxB的向量积就是x1y2-y1x2。
我们假定一个向量积R,R=x1y2-y1x2。
R<0 说明A在B的逆时针方向
R=0 说明A与B共线,可能正向也可能反向
R>0 说明A在B的顺时针方向

判断要两组,以A为顶点的三条向量 + 以C为顶点的三条向量

若此时A线段向量在L1,L2的中间或L1,L2的边上,就能说明B跨立A
即L1,L2在A的不同的顺逆时针方向,就可以分别求出两个L的向量积,再将他们相乘,如果结果<=0,即向量积异号或有0。

判断规范相交

 判断不规范相交

最后就是用二次函数y=kx+b或者先前有的向量叉乘直接求点坐标

 

 下午做得人都要傻掉了也没做出来,痛、太痛了

AC的代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=30;
const double eps=1e-6;
double ax,ay,cx,cy,dx,dy,bx,by;
int check(){
    if(min(ax,bx)>max(cx,dx)||min(ay,by)>max(cy,dy))
        return 0;
    if(min(cx,dx)>max(ax,bx)||min(cy,dy)>max(ay,by))
        return 0;
    return 1;
}
double ex(double x0,double y0,double x1,double y1,double x2,double y2){
    return (x1-x0)*(y2-y0)-(x2-x0)*(y1-y0);
}
int judge(double x){
    if(fabs(x)<eps)
        return 0;
    return x>0? 1:-1;
}
void find(){
    double a1=dy-cy,a2=cy-ay,a3=by-ay;
    double b1=dx-cx,b2=cx-ax,b3=bx-ax;
    double t=(a3*b2-a2*b3)/(a1*b3-a3*b1);
    double y=a1*t+cy;
    double x=b1*t-cx;
    printf(" (%.3f,%.3f)",x,y);
    return;
}
void solve(){
    //(A-C) x (D-C) * (B-C) x (D-C) < 0
    //C D 在向量ab的两侧 
    //(C-A) x (B-A) * (D-A) x (B-A) < 0
    double s1=ex(ax,ay,bx,by,cx,cy);
    double s2=ex(ax,ay,bx,by,dx,dy);
    double s3=ex(cx,cy,dx,dy,ax,ay);
    double s4=ex(cx,cy,dx,dy,bx,by);
    int d1=judge(s1),d2=judge(s2),d3=judge(s3),d4=judge(s4);
    //判断规范相交
    if((d1*d2<0)&&(d3*d4<0)){
        //find(); //交点坐标 
        double x=(cx*s2-dx*s1)/(s2-s1);
        double y=(cy*s2-dy*s1)/(s2-s1);
        printf("YES (%.3f,%.3f)",x,y); 
        return;
    }
    //判断非规范相交
    //d1==0, 则证明a,b,c三点共线; 
    s1=ex(cx,cy,ax,ay,bx,by);
    s2=ex(dx,dy,ax,ay,bx,by);
    s3=ex(ax,ay,cx,cy,dx,dy);
    s4=ex(bx,by,cx,cy,dx,dy);
    int k1=judge(s1),k2=judge(s2),k3=judge(s3),k4=judge(s4);
    if(d1==0&&k1<=0||d2==0&&k2<=0||d3==0&&k3<=0||d4==0&&k4<=0){
        cout<<"YES";
        return;
    }
    cout<<"NO";
    return;
}
int main()
{
    int t;
    cin>>t;
    while(t--){
        cin>>ax>>ay>>bx>>by;
        cin>>cx>>cy>>dx>>dy;
        if(!check()){
            cout<<"NO";
            if(t) cout<<endl;
            continue;
        }
        solve();
        if(t) cout<<endl;
    }
    return 0;
}

这篇关于TZOJ 3348: 线段相交Ⅲ(叉乘+快速排斥实验+跨立实验)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/190436

相关文章

使用Python实现快速搭建本地HTTP服务器

《使用Python实现快速搭建本地HTTP服务器》:本文主要介绍如何使用Python快速搭建本地HTTP服务器,轻松实现一键HTTP文件共享,同时结合二维码技术,让访问更简单,感兴趣的小伙伴可以了... 目录1. 概述2. 快速搭建 HTTP 文件共享服务2.1 核心思路2.2 代码实现2.3 代码解读3.

springboot security快速使用示例详解

《springbootsecurity快速使用示例详解》:本文主要介绍springbootsecurity快速使用示例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录创www.chinasem.cn建spring boot项目生成脚手架配置依赖接口示例代码项目结构启用s

C++快速排序超详细讲解

《C++快速排序超详细讲解》快速排序是一种高效的排序算法,通过分治法将数组划分为两部分,递归排序,直到整个数组有序,通过代码解析和示例,详细解释了快速排序的工作原理和实现过程,需要的朋友可以参考下... 目录一、快速排序原理二、快速排序标准代码三、代码解析四、使用while循环的快速排序1.代码代码1.由快

Win32下C++实现快速获取硬盘分区信息

《Win32下C++实现快速获取硬盘分区信息》这篇文章主要为大家详细介绍了Win32下C++如何实现快速获取硬盘分区信息,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 实现代码CDiskDriveUtils.h#pragma once #include <wtypesbase

Spring AI与DeepSeek实战一之快速打造智能对话应用

《SpringAI与DeepSeek实战一之快速打造智能对话应用》本文详细介绍了如何通过SpringAI框架集成DeepSeek大模型,实现普通对话和流式对话功能,步骤包括申请API-KEY、项目搭... 目录一、概述二、申请DeepSeek的API-KEY三、项目搭建3.1. 开发环境要求3.2. mav

Python如何快速下载依赖

《Python如何快速下载依赖》本文介绍了四种在Python中快速下载依赖的方法,包括使用国内镜像源、开启pip并发下载功能、使用pipreqs批量下载项目依赖以及使用conda管理依赖,通过这些方法... 目录python快速下载依赖1. 使用国内镜像源临时使用镜像源永久配置镜像源2. 使用 pip 的并

SpringBoot快速接入OpenAI大模型的方法(JDK8)

《SpringBoot快速接入OpenAI大模型的方法(JDK8)》本文介绍了如何使用AI4J快速接入OpenAI大模型,并展示了如何实现流式与非流式的输出,以及对函数调用的使用,AI4J支持JDK8... 目录使用AI4J快速接入OpenAI大模型介绍AI4J-github快速使用创建SpringBoot

使用Python快速实现链接转word文档

《使用Python快速实现链接转word文档》这篇文章主要为大家详细介绍了如何使用Python快速实现链接转word文档功能,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 演示代码展示from newspaper import Articlefrom docx import

shell脚本快速检查192.168.1网段ip是否在用的方法

《shell脚本快速检查192.168.1网段ip是否在用的方法》该Shell脚本通过并发ping命令检查192.168.1网段中哪些IP地址正在使用,脚本定义了网络段、超时时间和并行扫描数量,并使用... 目录脚本:检查 192.168.1 网段 IP 是否在用脚本说明使用方法示例输出优化建议总结检查 1

Rust中的Option枚举快速入门教程

《Rust中的Option枚举快速入门教程》Rust中的Option枚举用于表示可能不存在的值,提供了多种方法来处理这些值,避免了空指针异常,文章介绍了Option的定义、常见方法、使用场景以及注意事... 目录引言Option介绍Option的常见方法Option使用场景场景一:函数返回可能不存在的值场景