为什么在共线三点计数算法中必须使用 GCD 归一化斜率?

梦雪吖_2420

梦雪吖_2420

2026-10-09

704人浏览

原创

为什么在共线三点计数算法中必须使用 GCD 归一化斜率?

本文详解 gcd 在判断三点共线问题中的关键作用:它通过将方向向量约分为最简整数比,避免浮点精度误差与重复表示问题,确保斜率比较的数学严谨性与计算可靠性。

本文详解 gcd 在判断三点共线问题中的关键作用:它通过将方向向量约分为最简整数比,避免浮点精度误差与重复表示问题,确保斜率比较的数学严谨性与计算可靠性。

在解决“给定平面上若干整数坐标点,统计有多少组三点共线”这一经典计算几何问题时,一个常见思路是:对每个基准点 (P_i),枚举其余点 (P_j),计算向量 (\vec{P_iP_j} = (dx, dy)),并按其方向(即斜率)分组——因为所有与 (P_i) 共线的点,其相对于 (P_i) 的方向向量必成比例,即存在非零实数 (k) 使得 ((dx', dy') = k \cdot (dx, dy))。

然而,若直接用浮点数表示斜率(如 dy / dx),会引入两类致命缺陷:

  1. 浮点精度丢失:当坐标值很大时(例如接近 (10^9)),两个本质不同的分数(如 (999999997/999999998) 和 (999999998/999999999))在双精度浮点下可能被四舍五入为同一值,导致错误归并;
  2. 表示不唯一:向量 ((2,4))、((3,6))、((-1,-2)) 均代表相同方向,但若直接存储 (dy, dx) 或 dy/dx,它们无法自然哈希等价;而浮点除法还额外引入正负零、NaN 等边界问题。

✅ GCD 的作用正是消除上述问题:
对任意向量 ((dx, dy))(注意:需统一处理符号,通常约定让 dx 非负,若 dx == 0 则令 dy > 0),计算 (g = \gcd(|dx|, |dy|)),再将向量约简为最简整数比:

g = gcd(abs(dx), abs(dy))  # 实际代码中 gcd 已支持负数,但逻辑等价
norm_dx = dx // g
norm_dy = dy // g
# 进一步标准化符号:确保 norm_dx >= 0;若 norm_dx == 0,则 norm_dy > 0
if norm_dx <p>这样,所有共线向量均映射到<strong>唯一的、规范化的整数对</strong>,可安全用于字典键(<code>defaultdict[int]</code>)计数。</p><p>⚠️ 特别注意边界情况:</p>
  • 当 dx == 0(垂直线):此时 gcd(0, dy) = |dy|,约简后得 (sign(dy), 0) → 标准化为 (1, 0);
  • 当 dy == 0(水平线):约简后为 (0, 1);
  • 当 dx == dy == 0:不会发生(因点互异,且 j > i,两点不同)。

最终,对每个基准点 (P_i),若某斜率方向上有 k 个其他点,则这 k 个点两两与 (P_i) 构成共线三点组,组合数为 (\binom{k}{2})。累加所有方向即可得全局答案。

? 总结:GCD 不是“可选优化”,而是保障算法数学正确性的必要步骤——它将无限精度的方向关系,无损压缩为有限、唯一、可哈希的整数标识,彻底规避浮点陷阱与表示歧义。在涉及整数坐标的几何计数问题中,这种基于 GCD 的方向归一化是工业级实现的标准范式。

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
页面置换算法
页面置换算法

页面置换算法是操作系统中用来决定在内存中哪些页面应该被换出以便为新的页面提供空间的算法。本专题为大家提供页面置换算法的相关文章,大家可以免费体验。

2023.08.14

5396

4

FrankenPHP集成Laravel详细教程
FrankenPHP集成Laravel详细教程

本专题提供FrankenPHP集成Laravel的详细配置指南,全面解析运行原理、开发环境搭建、Caddyfile配置、Octane工作模式、数据库连接、队列任务、定时任务和生产环境优化,解决部署过程中常见的报错与兼容性问题。

2026.10.08

20

20

LLVM自定义Pass怎么写
LLVM自定义Pass怎么写

本专题聚焦LLVM自定义Pass开发,整理Pass类结构、run()方法、PreservedAnalyses、CMake构建、插件注册、-load-pass-plugin加载和测试用例编写流程。

2026.09.30

120

10

LLVM RISC-V参数配置教程
LLVM RISC-V参数配置教程

本专题介绍LLVM对RISC-V基础ISA和扩展的支持方式,涵盖RV32、RV64、标准扩展、实验性扩展、厂商扩展、-menable-experimental-extensions和版本差异。

2026.09.30

100

14

LLVM IR中间表示入门指南
LLVM IR中间表示入门指南

本专题整理LLVM IR的核心概念,包括中间表示作用、模块结构、函数、基本块、SSA形式、类型系统和常见语法,帮助新手理解LLVM编译流程中的关键层。

2026.09.30

80

12

PDF转图片方法
PDF转图片方法

需要把 PDF 页面用于上传、预览、分享或图片归档时,PDF 转图片方法专题整理 JPG/PNG 格式选择、逐页导出、清晰度设置、批量下载和结果检查等流程,帮助用户稳定完成 PDF 图片化处理。

2026.09.30

80

26

PixTV AI视频生成与无限画布创作
PixTV AI视频生成与无限画布创作

PixTV专题整理AI视频与视觉内容创作相关功能使用教程,涵盖AI生图、视频生成、无限画布、多模型创作、素材管理、声音音乐及视频剪辑等功能,帮助用户快速掌握PixTV从创意到成片的完整制作方法。

2026.09.29

100

15

Buffalo框架数据库开发全教程
Buffalo框架数据库开发全教程

本专题围绕Buffalo框架数据库开发,讲解database.yml多环境配置、soda与fizz迁移生成回滚、模型结构体标签、增删改查与条件查询、一对多与多对多关联、数据校验、回调钩子、事务处理及原生SQL执行能力。

2026.09.23

300

15

Buffalo框架路由与请求处理实操指南
Buffalo框架路由与请求处理实操指南

本专题讲解Buffalo框架路由与请求处理机制,涵盖路由注册与分组、资源路由、Handler编写规范、Context上下文方法、参数绑定、中间件编写挂载、Session与Cookie读写、Flash消息及错误页面定制方法。

2026.09.23

180

15

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
热门推荐
/
最新课程
phpStudy极速入门视频教程
phpStudy极速入门视频教程

共6课时 | 54.6万人学习

独孤九贱(4)_PHP视频教程
独孤九贱(4)_PHP视频教程

共89课时 | 133.4万人学习