每天学一点Python之NumPy基础11质因数
off999 2024-12-11 16:24 12 浏览 0 评论
质因数 (http://en.wikipedia.org/wiki/Prime_factor) 是指正好整除一个整数而不留余数的质数。对于大数来说，寻找质因数几乎是不可能的。
因此，质因数在密码学中得到了应用。然而，使用正确的算法--费马因式分解法（http://en.wikipedia.org/wiki/Fermat%27s_factorization_method）和 NumPy--对于小数来说，因式分解变得相对容易。
其原理是根据下式将一个数 N 分解成两个数 c 和 d：
我们可以递归地应用因式分解，直到得到所需的质因数。
这需要四个步骤，让我们来看看是如何完成的。
算法要求我们对 a 的值进行多次试验。
首先 我们将创建一个包含试验值的数组。创建一个 NumPy 数组可以省去循环。
不过，在创建数组时要注意，不要创建占用内存过大的数组。在我的系统中，一百万个元素的数组大小似乎恰到好处：
a = np.ceil(np.sqrt(n))
lim = min(n, LIM)
a = np.arange(a, a + lim)
b2 = a ** 2 - n我们使用 ( ceil ) 方法按元素顺序返回输入值的上限。
第二步 获取 b 数组的小数部分。现在我们要检查 ( b ) 是否是正方形。使用名为 ( modf ) 的 NumPy 方法来获取 ( b ) 的小数部分。
数组的小数部分：
fractions = np.modf(np.sqrt(b2))[0]第三步 查找 0 分数。调用 ( where ) NumPy 方法查找零分数的索引，其中分数部分为 0：
indices = np.where(fractions == 0)第四步 找出第一个出现的零分数。
首先，使用上一步的索引数组调用 NumPy 方法 ( take )，获取零分率的值。然后使用 NumPy 方法 ( ravel ) 将数组网格化：
a = np.ravel(np.take(a, indices))[0]这一行有点复杂，但确实展示了两种有用的方法。如果写成
a = a[indices][0]你现在需要知道发生了什么？
在本节中，我们应用了一组有用且有趣的 NumPy 方法，这些方法的说明如下：
Method Name | Method Description |
ceil() | 计算数组元素的上限（参见 http://docs.scipy.org/doc/numpy/reference/generated/numpy.ceil.html） |
modf() | 返回浮点数的小数部分和积分部分（参见 http://docs.scipy.org/doc/numpy/reference/generated/ numpy.modf.html） |
where() | 根据条件返回数组索引（参见 http://docs.scipy.org/doc/numpy/reference/generated/numpy.where.html） |
ravel() | 返回一个扁平化数组（参见 http://docs.scipy.org/doc/numpy/reference/generated/numpy.ravel.html） |
take() | 从数组中提取元素（参见 http://docs.scipy.org/doc/numpy/reference/generated/numpy.take.html） |
- 上一篇：Python 初学者练习:验证哥德巴赫猜想
- 下一篇：详解 Python 中的生成器
相关推荐
- 独家 | 5 个Python高级特性让你在不知不觉中成为Python高手
-
你已经使用Python编程了一段时间，编写脚本并解决各种问题。是你的水平出色吗？你可能只是在不知不觉中利用了Python的高级特性。从闭包（closure）到上下文管理器（contextmana...
- Python装饰器
-
Python装饰器是一种用于修改函数或类的行为的特殊语法。它们允许在不修改原始代码的情况下，通过将函数或类作为参数传递给另一个函数来添加额外的功能。装饰器本质上是一个函数，它接受一个函数作为参数，并返...
- 中高阶Python常规用法--上下文管理器
-
Python以简单性和通用性著称，是一种深受全球开发人员喜爱的编程语言。它提供了大量的特性和功能，使编码成为一种愉快的体验。在这些功能中，一个经常被新手忽视的强大工具是上下文管理器。上下文管理器是高...
- Python小案例67- 装饰器
-
Python装饰器是一种用于修改函数或类的行为的特殊语法。它们允许在不修改原始代码的情况下，通过将函数或类作为参数传递给另一个函数来添加额外的功能。装饰器本质上是一个函数，它接受一个函数作为参数，并返...
- python常用的语法糖
-
概念Python的语法糖（SyntacticSugar）是指那些让代码更简洁、更易读的语法特性，它们本质上并不会增加新功能，但能让开发者更高效地编写代码。推导式写法推导式是Python最经典的...
- python - 常用的装饰器 decorator 有哪些?
-
python编程中使用装饰器(decorator)工具，可以使代码更简洁清晰，提高代码的重用性，还可以为代码维护提供方便。对于python初学者来说，根据装饰器(decorator)的字面意思并不...
- python数据缓存怎么搞 ?推荐一个三方包供你参考，非常简单好用。
-
1.数据缓存说明数据缓存可以说也是项目开发中比不可少的一个工具，像我们测试的系统中，你都会见到像Redis一样的数据缓存库。使用缓存数据库的好处不言而喻，那就是效率高，简单数据直接放在缓存中...
- 用于时间序列数据的Graphite监视工具
-
结合第三方工具，Graphite为IT性能监控提供了许多好处。本文介绍其核心组件，包括Carbon、Whisper以及安装的基本准则。Graphite监视工具可实时或按需，大规模地绘制来自多个来源的时...
- Python3+pygame实现的坦克大战
-
一、显示效果二、代码1.说明几乎所有pygame游戏，基本都遵循一定的开发流程，大体如下：初始化pygame创建窗口while循环检测以及处理事件（鼠标点击、按键等）更新UI界面2.代码创建一个m...
- Python之鸭子类型:一次搞懂with与上下文装饰器
-
引言在鸭子类型的理念的基础之上，从关注类型，转变到关注特性和行为。结合Python中的魔法函数的体系，我们可以将自定义的类型，像内置类型一样被使用。今天这篇文章中，接着该话题，继续聊一下with语法块...
- Python必会的50个代码操作
-
学习Python时，掌握一些常用的程序操作非常重要。以下是50个Python必会的程序操作，主要包括基础语法、数据结构、函数和文件操作等。1.HelloWorldprint("Hello,...
- 一文掌握Python 中的同步和异步
-
同步代码（Sync）同步就像在一个流水线上工作，每个任务都等待前一个任务完成。示例：机器A切割钢板→完成后，机器B钻孔→完成后，机器C上色。在Python中，同步代码看起来像这样：im...
- python 标注模块timeit: 测试函数的运行时间
-
在Python中，可以使用内置的timeit模块来测试函数的运行时间。timeit模块提供了一个简单的接口来测量小段代码的执行时间。以下是使用timeit测试函数运行时间的一般步骤：导入...
- Python带你找回童年的万花尺
-
还记得小时候的万花尺吧？这么画：一点也不费脑筋，就可以出来这么多丰富多彩的复杂几何图形。具体而言，可以用万花尺玩具（如图2-1所示）来绘制数学曲线。这种玩具由两个不同尺寸的塑料齿轮组成，一大一小。小的...
- Python 时间模块深度解析:从基础到高级的全面指南
-
直接上干货一、时间模块核心类介绍序号类名说明1datetime.datetime表示一个具体的日期和时间，结合了日期和时间的信息。2datetime.date表示一个具体的日期。3datetime.t...
欢迎 你 发表评论:
- 一周热门
- 最近发表
- 标签列表
-
- python计时 (54)
- python安装路径 (54)
- python类型转换 (75)
- python进度条 (54)
- python的for循环 (56)
- python串口编程 (60)
- python写入txt (51)
- python读取文件夹下所有文件 (59)
- java调用python脚本 (56)
- python操作mysql数据库 (66)
- python字典增加键值对 (53)
- python获取列表的长度 (64)
- python接口 (63)
- python调用函数 (57)
- python qt (52)
- python人脸识别 (54)
- python斐波那契数列 (51)
- python多态 (60)
- python命令行参数 (53)
- python匿名函数 (59)
- python打印九九乘法表 (65)
- centos7安装python (53)
- python赋值 (62)
- python异常 (69)
- python元祖 (57)