每天学一点Python之NumPy基础11质因数
off999 2024-12-11 16:24 14 浏览 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 中的生成器
相关推荐
- 用Python编制生成4位数字字母混合验证码
-
我们登录一些网站、APP的时候经常会有验证码，这个为了防止有人不停的去试探密码，还有发送短信验证之前，输入验证码就可以减少误点，错误操作等等。可以提高安全性，我们可以生成数字，也可以生成字母，也可...
- Python电子发票管理工具4:前后端业务逻辑实现
-
用一系列文章介绍如何用python写一个发票管理小工具。在前面的文章中前端页面和后端框架已经实现，本文将介绍功能实现的代码。数据库操作使用sqlalchemy操作sqlite数据库。sqlalchem...
- 【代码抠图】4行Python代码帮你消除图片背景
-
在修图工具满天飞的年代其实仍然还有很多人不会扣图（比如我），在很多需要去除某些照片上面的背景的时候就会很难受，所以今天就给不会扣图的小伙伴们来带一个简单的代码扣图教程，只需要4行代码，不用再多了。准备...
- Python 3.14 的 UUIDv6/v7/v8 上新，别再用 uuid4 () 啦!
-
Python3.14重磅更新！UUIDv6/v7/v8强势来袭，别再用uuid4()啦！为什么说UUID升级是2025年Python开发者的必学技能？在当今互联网应用中，UU...
- 殊途同归 python 第 4 节:有趣的键值对(字典)
-
字典数据的突出特点就是“键”和“值”，前文已经简单介绍过，本文来聊聊关于字典的几个高级玩法。1.函数打包后，通过键来调用globalf1,f2a={"k1":f1,"k2...
- 更有效地使用 Python Pandas 的 4 个技巧
-
一个简单而实用的指南照片由simonsun在Unsplash上拍摄Pandas是一个用于数据分析和操作任务的非常实用且功能强大的库。自2019年以来，我一直在使用Pandas，它始终能够为我...
- 4.python学习笔记-集合（python里面集合）
-
1.关于集合集合是一类元素无序不重复的数据结构，常用场景是元素去重和集合运算。python可以使用大括号{}或者set()函数创建集合，如果创建一个空集合必须用set()而不是{}，因为{}是用来表示...
- python生成4种UUID（python随机生成uuid）
-
总结了一份python生成4种UUID的代码：UUID用4种uuid生成方法：uuid1：基于时间戳由MAC地址、当前时间戳、随机数字。保证全球范围内的唯一性。但是由于MAC地址使用会带来安全问题...
- 你不知道的4种方法:python方法绘制扇形
-
1说明：=====1.1是问答中的我的一个回答。1.1因为问答中没有代码块的，所以我改为这里写文章，然后链接过去。1.24种方法：turtle法、OpenCV法、pygame法和matplot...
- 30天学会Python编程:4. Python运算符与表达式
-
4.1运算符概述4.1.1运算符分类Python运算符可分为以下几大类：4.1.2运算符优先级表4-1Python运算符优先级（从高到低）运算符描述示例**指数2**3→8~+-按位取...
- 这3个高级Python函数，不能再被你忽略了
-
全文共1657字，预计学习时长3分钟Python其实也可以带来很多乐趣。重新审视一些一开始并不被人们熟知的内置函数并没有想象中那么难，但为什么要这么做呢?今天，本文就来仔细分析3个在日常工作中或多或少...
- beautifulSoup4，一个超实用的python库
-
一.前言我们在学习python爬虫的时候，数据提取是一个常见的任务。我们一般使用正则表达式，lxml等提取我们需要的数据，今天我们介绍一个新的库beautifulSoup4，使用它您可以从HTML和...
- AI指导:打造第一个Python应用(4)（python ai开发）
-
眼瞅着迈过几个里程碑，与目标越来越近。尽管过程中照旧因返工而心焦，而欣喜与急躁比例，是喜悦运大于焦虑。从初次熟悉智能大模型，尝试编程起步，不定期进行复盘反思，这是小助手指导编程的第四篇。复盘以为记。需...
- wxPython 4.2.0终于发布了（wxpython安装教程）
-
wxPython是Python语言的跨平台GUI工具包。使用wxPython，软件开发人员可以为他们的Python应用程序创建真正的本地用户界面，这些应用程序在Windows、Ma...
- 《Python学习手册(第4版)》PDF开放下载，建议收藏
-
书籍简介如果你想动手编写高效、高质量并且很容易与其他语言和工具集成的代码，本书将快速地帮助你利用Python提高效率。本书基于Python专家的流程培训课程编写，内容通俗易懂。本书包含很多注释的例子和...
欢迎 你 发表评论:
- 一周热门
- 最近发表
- 标签列表
-
- python计时 (73)
- python安装路径 (56)
- python类型转换 (93)
- python进度条 (67)
- python吧 (67)
- python字典遍历 (54)
- python的for循环 (65)
- python格式化字符串 (61)
- python静态方法 (57)
- python列表切片 (59)
- python面向对象编程 (60)
- python 代码加密 (65)
- python串口编程 (60)
- python读取文件夹下所有文件 (59)
- java调用python脚本 (56)
- python操作mysql数据库 (66)
- python获取列表的长度 (64)
- python接口 (63)
- python调用函数 (57)
- python多态 (60)
- python匿名函数 (59)
- python打印九九乘法表 (65)
- python赋值 (62)
- python异常 (69)
- python元祖 (57)