经典算法:最短点对

                                                                          软件架构师何志丹

说明

旧文新发,改了错别字,死链等。尽量保持“原汁原味”。

难点

如何测试。我的解决方式是:a,三种解法,看结果是否一致。b,小数据(100个点),人工排查。第一种方法,暴力法适合小数据。第二种方法:我的改进型。第三种方法:经典方法(分治法)。实验证明1000万数据时,我的算法有优势。

改进


暴力算法,O(n2)。我的改进型要点:先对所有数据按Y排序。只比较y距离小于等于已知最小距离的点对。经典方法:按Y排序,分成两部分,递归调用。合并时只比较距离分界线不超过已知最小距离的点对。
实际证明500万数据以下,我的改进算法明显优于经典算法;1000万数据时,略强于经典算法。

核心源码节选:

double Dis(const CPT& pt1,const CPT& pt2)
{return sqrt((double) (pt1.x-pt2.x)*(pt1.x-pt2.x)+(pt1.y-pt2.y)*(pt1.y-pt2.y)+(pt1.z-pt2.z)*(pt1.z-pt2.z) );
}void InitData(CPT* pts,int iNum)
{srand(time(NULL));for( int i = 0 ; i < iNum ; i++){pts[i].x = rand()%10000;pts[i].y = rand()%10000;pts[i].z = rand()%10000;}
}double Fun1(CPT* pts,const int iNum)
{double dMinDis = 10000*10000 ;for(int i = 0 ; i < iNum ; i++ )for( int j = i+1 ; j < iNum ; j++ ){const double d = Dis(pts[i] , pts[j]);if( d < dMinDis){dMinDis = d ;}}return dMinDis;
}class CCmpY
{
public:bool operator()(const CPT& pt1,const CPT& pt2){return pt1.y < pt2.y ;}
};double Fun2(CPT* pts,const int iNum)
{std::sort(pts,pts+iNum,CCmpY() );double dMinDis = 10000*10000 ;for(int i = 0 ; i < iNum ; i++ )for( int j = i+1 ; j < iNum ; j++ ){const double d = Dis(pts[i] , pts[j]);if( d < dMinDis){dMinDis = d ;}if( abs(pts[i].y - pts[j].y )> dMinDis ){break;}}return dMinDis;
}double Fun3(CPT* pts,const int iNum)
{std::sort(pts,pts+iNum,CCmpY() );if( iNum < 100 ){return Fun1(pts,iNum);}const int iMid = iNum/2 ;const double dMin1 = Fun3(pts,iMid);const double dMin2 = Fun3(pts+iMid,iNum-iMid);double dMinDis = min(dMin1,dMin2) ;for(int i = iMid-1 ; i >= 0 ; i-- )//左集合{if( abs(pts[i].y - pts[iMid].y ) > dMinDis ){break;}for( int j = iMid ; j < iNum ; j++ )//右集合{const double d = Dis(pts[i] , pts[j]);if( d < dMinDis){dMinDis = d ;}if( abs(pts[i].y - pts[j].y )> dMinDis ){break;}}}return dMinDis;
}

测试环境

似乎是WinXP VS2005(VC8)

下载

可通过以下链接下载测试数据,exe,源码(VS2005,VC8)
https://download.csdn.net/download/he_zhidan/10887801

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/89039.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【无标题】C语言学习笔记5--循环结构和选择结构

5.1 运算符和运算符 5.1.1逻辑运算发 运算符说明结核性举例&&与运算&#xff0c;双目&#xff0c;对应数学中的“且”左结合1&&0、(9>3)&&(b>a与运算&#xff0c;双目&#xff0c;对应数学中的“且”左结合1!非运算&#xff0c;单目&#xff0…

MySQL(6)LOCK和MVCC

一、锁的分类 按照锁的属性&#xff1a;读锁、写锁、共享锁、排它锁、悲观锁、乐观锁 按照锁的范围&#xff1a;表锁、页锁、间隙锁、临键锁、行锁 按照锁的作用&#xff1a;意向锁、意向共享锁、意向排它锁、IS锁、IX锁 二、MySQL为什么要有锁 锁是计算机协调多个进程或线程并…

vue-cli2项目运行时中断解决方案记录

目录 前情提要解决方案1. cmd设置【结论&#xff1a;不得行】2. npm cache verify 清除缓存【结论&#xff1a;不得行】3. 删除node_modules&#xff0c;再重新npm i安装【结论&#xff1a;不得行】4. 安装cross-env和increase-memory-limit【结论&#xff1a;不得行】5. 在sci…

Redis的安装与基本使用

文章目录 Linux 环境下安装Redis下载Redis 安装包解压安装包安装Redis进入redis安装包下编译并且安装到指定目录下 启动redis配置远程访问找到Redis.config文件 Windows 环境下安装Redis说明官方提供方式安装或启用WSL2在WSL&#xff08;Ubuntu&#xff09;上安装Redis启动Redi…

26052-2022 硬质合金管状焊条 思维导图

声明 本文是学习GB-T 26052-2022 硬质合金管状焊条. 而整理的学习笔记,分享出来希望更多人受益,如果存在侵权请及时联系我们 1 范围 本文件规定了硬质合金管状焊条的技术要求、试验方法、检验规则、标志、包装、运输、贮存、随行文 件和订货单内容。 本文件适用于钢件的表…

Hive【Hive(二)DML】

启动 hive 命令行&#xff1a; hive DML 数据操作 1、数据导入 1.1、向表中装载数据&#xff08;load&#xff09; 语法&#xff1a; hive> load data [local] inpath 数据的path [overwrite] into table student [partition (partcol1val1,…)];&#xff08;1&#x…

修炼k8s+flink+hdfs+dlink(一:安装flink)

一&#xff1a;standalone的ha环境部署。 创建目录&#xff0c;上传安装包。 mkdir /opt/app/flink 上传安装包到本目录。 tar -zxvf flink-1.13.6-bin-scala_2.12.tgz配置参数。 在flink-conf.yaml中添加zookeeper配置 jobmanager.rpc.address: node01 high-availability: …

R语言学习笔记

R语言学习笔记 一.准备环境二.认识控制台三.R包四.数据结构1.向量Vector1.1创建向量1.2访问向量中的数据1.3向量的循环补齐 2.矩阵matrix2.1创建矩阵2.2访问矩阵中的数据 3数组Array3.1创建数组3.2访问数组中的数据 4.数据框Dataframe4.1创建数据框4.2访问数据框中的数据 5因子…

函数指针与回调函数

(1) 1分钟学会函数指针与回调函数&#xff08;上&#xff09;_哔哩哔哩_bilibili #include <QCoreApplication> #include <stdio.h>int sum(int x,int y){ return xy; } int difference(int x,int y){ return x-y; } int product(int x,int y){ return x*y; }in…

GitLab多人开发步骤

目录 一、基于develop创建自己的feature分支二、提交代码格式简易版 三、提交到远程仓库四、提交合并请求 一、基于develop创建自己的feature分支 使用git checkout切换到develop分支 git checkout develop 基于develop分支创建feature分支 格式&#xff1a; git checkout -b …

进阶JS-作用域和作用域链

作用域规定了变量能够被访问的“范围”&#xff0c;离开了这个“范围”变量便不能被访问。 作用域分为&#xff1a;局部作用域、全局作用域 局部作用域&#xff08;Local&#xff09; 局部作用域分为函数作用域和块作用域。 函数作用域&#xff1a;在函数内部声明的变量只能…

MYSQL操作命令大全

mysql命令大全 命令说明示例mysql -h ip -P port -u user -p password远程连接mysqlmysql -h192.168.0.201 -P3306 -uroot -p123mysql -u user -p password本地连接mysqlmysql -uroot -p123select * from information_schema.TABLES where TABLE_SCHEMA(select database()) 查…

完全指南:mv命令用法、示例和注意事项 | Linux文件移动与重命名

文章目录 mv命令使用指南1. 简介什么是mv命令&#xff1f;mv命令的作用和功能是什么&#xff1f; 2. 基本用法基本语法格式如何移动文件&#xff1f;如何重命名文件&#xff1f;如何移动和重命名目录&#xff1f; 3. 高级用法使用通配符进行批量移动和重命名使用选项进行文件移…

数据库设计步骤

一、数据库设计重点&#xff1a;有效存储、高效访问 1、为什么要进行数据库设计&#xff1a; &#xff08;1&#xff09;减少数据冗余&#xff1b; &#xff08;2&#xff09;避免数据异常维护&#xff1b; &#xff08;3&#xff09;节约存储空间&#xff1b; &#xff08;4&…

【CMU15-445 Part-12】Query Execution I

Part12-Query Execution I Processing Models Processing Model主要指的是明确如何去执行一个查询计划&#xff08;top 2 bottom or bottom 2 top,operator之间的传递&#xff09;。 Iterator Model (volcano model/pipeline model);每个算子实现一个Next( )&#xff0c;父…

vim缓存-交换文件

Catf1agCTF靶场 web swp 题目链接&#xff1a;http://catf1ag.cn/ 个人博客&#xff1a;https://sword-blogs.com/ 题目考点&#xff1a; vim在编辑文档的过程中如果异常退出&#xff0c;会产生缓存文件 vim 交换文件名 参考文章&#xff1a;vim手册 https://yianwillis.…

React 全栈体系(十六)

第八章 React 扩展 五、Context 1. 代码 /* index.jsx */ import React, { Component } from react import ./index.css//创建Context对象 const MyContext React.createContext() const {Provider,Consumer} MyContext export default class A extends Component {state …

Docker最基本使用

1 安装&#xff1a; sudo apt-get -y install docker.io测试&#xff1a; sudo docker run hello-world成功&#xff1a; Hello from Docker! This message shows that your installation appears to be working correctly.2 查看 查看已有镜像&#xff1a; sudo docker i…

Linux进程概念

文章目录 前言一、操作系统1、概念2、设计OS的目的3、总结 二、进程1、基本概念2、查看进程2.1 使用ps axj命令2.2 通过 /proc 系统文件夹查看 3、通过系统调用获取进程标示符3.1 getpid() 系统调用3.2 getppid() 系统调用 4、通过系统调用fork创建子进程4.1 使用fork创建子进程…

Linux 常用基本命令

1.走近Linux系统 开机登录 开机会启动许多程序。它们在Windows叫做"服务"&#xff08;service&#xff09;&#xff0c;在Linux就叫做"守护进程"&#xff08;daemon&#xff09;。 关机 在linux领域内大多用在服务器上&#xff0c;很少遇到关机的操作。毕竟…