探索三维几何：计算空间两线段的公共垂直线与分数类算法优化

作者：叶子美容美体养生馆os | 来源：互联网 | 2024-10-23 17:53

分析：首先判断线段俩直线是否平行（或重合），如果是的话直接求。考虑4个端点到另外一条线段的距离，取最小值即可

分析&＃xff1a;

首先判断线段俩直线是否平行&＃xff08;或重合&＃xff09;&＃xff0c;如果是的话直接求。考虑4个端点到另外一条线段的距离&＃xff0c;取最小值即可。

如果不平行或重合&＃xff0c;说明俩条直线是异面直线&＃xff0c;这时最短距离既可能是某端点到另外一条线段的距离&＃xff0c;也可能是异面直线的最短距离。

如何求异面直线的最短距离&＃xff1f;假设俩条直线分别为l1&＃61;(p1,v1&＃xff09;和l2&＃61;(p2,v2),那么最短距离会在某个q1&＃61;p1&＃43;sv1和q2&＃61;p2&＃43;tv2上取到&＃xff0c;其中q1和q2分别在l1和l2上&＃xff0c;且q1q2是这俩条异面直线的公垂线。

向量q1q1&＃61;q2-q1&＃61;p2-p1&＃43;tv2-sv1垂直于V1,因此Dot(p2-p1&＃43;tv2-sv1&＃xff0c;v1)&＃61;0,根据分配率&＃xff0c;有Dot(p2-p1,v1)&＃43;t*Dot(v2,v1)-s*Dot(v1,v1)&＃61;0.注意这里的三个点积都是可以直接算出来的&＃xff0c;因此实际上得到的是一个关于t和s的一次方程。根据q1q2垂直于v2还可以得到一个一次方程&＃xff0c;联立求解即可。下面的代码直接使用了经过复杂化简以后的结果。

注意本题比较特殊&＃xff0c;要求以分数方式输出。所以可以考虑定义一个Rat类&＃xff08;代表Rational&＃xff09;来保存和计算有理数&＃xff0c;并且重载加法、减法、和乘法,然后把本题用到的俩个关键函数&＃xff1a;点到线段距离Distace2Tosegment和异面直线的最小距离LineDistace3D改写成返回有理数的版本。

代码&＃xff1a;

#include #include #include #include #include #include #include #include #include #include #include #include #define LL long long using namespace std; const double eps&＃61;1e-6; int dcmp(double x) {return fabs(x)0?1:-1); } struct Point3 {LL x,y,z;Point3(LL x,LL y,LL z):x(x),y(y),z(z){}Point3() {}void in(){cin>>x>>y>>z;} }; typedef Point3 Vector3; Vector3 operator &＃43;(Vector3 A,Vector3 B) {return Vector3(A.x&＃43;B.x,A.y&＃43;B.y,A.z&＃43;B.z); } Vector3 operator -(Vector3 A,Vector3 B) {return Vector3(A.x-B.x,A.y-B.y,A.z-B.z); } Vector3 operator * (Vector3 A,int p) {return Vector3(A.x*p,A.y*p,A.z*p); } Vector3 operator / (Vector3 A,int p) {return Vector3(A.x/p,A.y/p,A.z/p); } bool operator &＃61;&＃61; (Vector3 A,Vector3 B) {return A.x&＃61;&＃61;B.x&&A.y&＃61;&＃61;B.y&&A.z&＃61;&＃61;B.z; } int Dot(Vector3 A,Vector3 B) {return A.x*B.x&＃43;A.y*B.y&＃43;A.z*B.z; } double Length(Vector3 A) {return sqrt(Dot(A,A)); } int Length2(Vector3 A) {return Dot(A,A); } Vector3 Cross(Vector3 A,Vector3 B) {return Vector3(A.y*B.z-A.z*B.y,A.z*B.x-A.x*B.z,A.x*B.y-A.y*B.x); } LL gcd(LL a,LL b){return b?gcd(b,a%b):a;} LL lcm(LL a,LL b){return a/gcd(a,b)*b;} struct Rat {LL a,b;Rat(LL a&＃61;0):a(a),b(1) {}Rat(LL x,LL y):a(x),b(y){if(b<0) a&＃61;-a,b&＃61;-b;LL d&＃61;gcd(a,b);if(d<0) d&＃61;-d;a/&＃61;d;b/&＃61;d;} }; Rat operator &＃43; (const Rat &A,const Rat &B) {LL x&＃61;lcm(A.b,B.b);return Rat(A.a*(x/A.b)&＃43;B.a*(x/B.b),x); } Rat operator - (const Rat &A,const Rat &B) {return A&＃43;Rat(-B.a,B.b); } Rat operator *(const Rat &A,const Rat &B) {return Rat(A.a*B.a,A.b*B.b); }void updatemin(Rat &A,const Rat &B) {if(A.a*B.b>B.a*A.b) A.a&＃61;B.a,A.b&＃61;B.b; } Rat Rat_Distace2ToSegment(const Point3 &P,const Point3 &A,const Point3 &B) {if(A&＃61;&＃61;B) return Length2(P-A);Vector3 v1&＃61;B-A,v2&＃61;P-A,v3&＃61;P-B;if(Dot(v1,v2)<0) return Length2(v2);else if(Dot(v1,v3)>0) return Length2(v3);else return Rat(Length2(Cross(v1,v2)),Length2(v1)); } bool Rat_LineDistace3D(const Point3 &p1,const Point3 &u,const Point3 &p2,const Vector3 &v, Rat &s) {LL b &＃61;(LL)Dot(u,u)*Dot(v,v)-(LL)Dot(u,v)*Dot(u,v);if(b&＃61;&＃61;0) return false;LL a&＃61;(LL)Dot(u,v)*Dot(v,p1-p2)-(LL)Dot(v,v)*Dot(u,p1-p2);s&＃61;Rat(a,b);return true; } void Rat_GetPointOnLine(const Point3 &A,const Point3 &B,const Rat &t,Rat &x,Rat &y,Rat &z) {x&＃61;Rat(A.x)&＃43;Rat(B.x-A.x)*t;y&＃61;Rat(A.y)&＃43;Rat(B.y-A.y)*t;z&＃61;Rat(A.z)&＃43;Rat(B.z-A.z)*t; } Rat Rat_Distance2(const Rat &x1,const Rat &y1,const Rat &z1,const Rat &x2,const Rat &y2,const Rat &z2) {return (x1-x2)*(x1-x2)&＃43;(y1-y2)*(y1-y2)&＃43;(z1-z2)*(z1-z2); } int main() {int T;cin>>T;while(T--){Point3 A,B,C,D;A.in();B.in();C.in();D.in();Rat s,t;bool ok&＃61;false;Rat ans&＃61;Rat(1000000000);if(Rat_LineDistace3D(A,B-A,C,D-C,s))if(s.a>0&&s.a0&&t.a}

推荐阅读

less
贪心与优先队列：最小化加法代价问题

本题要求在一组数中反复取出两个数相加，并将结果放回数组中，最终求出最小的总加法代价。这是一个经典的哈夫曼编码问题，利用贪心算法可以有效地解决。 ... [详细]

蜡笔小新 2024-12-20 23:20:38
function
深入理解Redis的数据结构与对象系统

本文详细探讨了Redis中的数据结构和对象系统的实现，包括字符串、列表、集合、哈希表和有序集合等五种核心对象类型，以及它们所使用的底层数据结构。通过分析源码和相关文献，帮助读者更好地理解Redis的设计原理。 ... [详细]

蜡笔小新 2024-12-25 04:11:22
go
图论问题解析：POJ2762 从u到v或从v到u的可达性判断（强连通分量缩点与单向连通性检测）

本文深入探讨了POJ2762问题，旨在通过强连通分量缩点和单向连通性的判断方法，解决有向图中任意两点之间的可达性问题。文章详细介绍了算法原理、实现步骤，并附带完整的代码示例。 ... [详细]

蜡笔小新 2024-12-24 10:44:24
char
如何高效清空C++标准输入缓冲区

本文探讨了在C++中如何有效地清空输入缓冲区，确保程序只处理最近的输入并丢弃多余的输入。我们将介绍一种不阻塞的方法，并提供一个具体的实现方案。 ... [详细]

蜡笔小新 2024-12-23 15:44:01
go
紫荆花之恋：动态树上的小精灵友谊问题

本题来自WC2014，题目编号为BZOJ3435、洛谷P3920和UOJ55。该问题描述了一棵不断生长的带权树及其节点上小精灵之间的友谊关系，要求实时计算每次新增节点后树上所有可能的朋友对数。 ... [详细]

蜡笔小新 2024-12-22 14:36:54
web
UC浏览器分享图片显示问题及解决方案

本文探讨了在UC浏览器中调用分享面板后，图片无法正常显示的问题，并提供了详细的解决方法和代码示例。 ... [详细]

蜡笔小新 2024-12-24 17:24:36
web
C#中SQL Server与Access批量数据插入性能对比

本文探讨了使用C#在SQL Server和Access数据库中批量插入多条数据的性能差异。通过具体代码示例，详细分析了两种数据库的执行效率，并提供了优化建议。 ... [详细]

蜡笔小新 2024-12-23 13:03:32
go
丽江客栈选择问题

本文介绍了一道经典的算法题，题目涉及在丽江河边的n家特色客栈中选择住宿方案。两位游客希望住在色调相同的两家客栈，并在晚上选择一家最低消费不超过p元的咖啡店小聚。我们将详细探讨如何计算满足条件的住宿方案总数。 ... [详细]

蜡笔小新 2024-12-22 20:15:25
go
JSOI2010 蔬菜庆典：树结构中的无限大权值问题

本文探讨了 JSOI2010 的蔬菜庆典问题，主要关注如何处理非根非叶子节点的无限大权值情况。通过分析根节点及其子树的特性，提出了有效的解决方案，并详细解释了算法的实现过程。 ... [详细]

蜡笔小新 2024-12-22 18:57:48
function
使用Fetch进行HTTP请求的基本示例

本文介绍了如何使用JavaScript的Fetch API与Express服务器进行交互，涵盖了GET、POST、PUT和DELETE请求的实现，并展示了如何处理JSON响应。 ... [详细]

蜡笔小新 2024-12-22 12:55:37
c语言
使用栈实现字符串中括号匹配的C++程序

本文介绍如何利用栈数据结构在C++中判断字符串中的括号是否匹配。通过顺序栈和链栈两种方式实现，并详细解释了算法的核心思想和具体实现步骤。 ... [详细]

蜡笔小新 2024-12-22 11:44:35
go
2-SAT问题学习笔记

本文介绍了一种解决二元可满足性（2-SAT）问题的方法。通过具体实例，详细解释了如何构建模型、应用算法，并提供了编程实现的细节和优化建议。 ... [详细]

蜡笔小新 2024-12-24 21:48:43
char
Tetris 排名系统 (拓扑排序与并查集的应用)

本题旨在通过给定的评级信息，利用拓扑排序和并查集算法来确定全球 Tetris 高手排行榜。题目要求判断是否可以根据提供的信息生成一个明确的排名表，或者是否存在冲突或信息不足的情况。 ... [详细]

蜡笔小新 2024-12-24 21:03:51
char
Linux 文件输入输出操作

本文介绍了Linux系统中的文件IO操作，包括文件描述符、基本文件操作函数以及目录操作。详细解释了各个函数的参数和返回值，并提供了代码示例。 ... [详细]

蜡笔小新 2024-12-24 02:33:04
char
深入解析Java多线程与并发库的应用：空中网实习生面试题详解

本文详细探讨了Java多线程与并发库的高级应用，结合空中网在挑选实习生时的面试题目，深入分析了相关技术要点和实现细节。文章通过具体的代码示例展示了如何使用Semaphore和SynchronousQueue来管理线程同步和任务调度。 ... [详细]

蜡笔小新 2024-12-21 09:43:07

叶子美容美体养生馆os

这个家伙很懒，什么也没留下！

Tags | 热门标签

RankList | 热门文章