当前位置:
首页 > Python基础教程 >
-
电影里的代码之《机械姬》:筛法求质数(2)
代码的最后打印出来下面这个很奇怪的东西,目测是一本书的ISBN,上豆瓣查了一下,是Embodiment and the Inner Life,是关于思维、意识的内容的,和本片的主题息息相关。
1
|
ISBN =9780199226559[Finished in 0.1s] |
重点还是前面的两个函数实现的筛法求质数。首先介绍一下什么是筛法,筛法相传是古希腊的埃拉托斯特尼发明的一种检测素数的算法。筛法的思路非常简单,可以用下面的动图来描述。给定一个范围,首先用2去筛,把所有2的倍数都筛掉,然后再用3筛,用5筛,不断重复下去......
再来看代码
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
def sieve(n): / / 对n以内的数进行筛选,返回一个长度为n的布尔数组 #compute primes using sieve eratosthenes x = [ 1 ] * n / / 定义长度为n的布尔数组(实际上电影里用 1 和 0 来表示true和false了) x[ 1 ] = 0 / / 1 既不是素数也不是合数,设为 0 for i in range ( 2 ,n / 2 ): / / i从 2 开始一直到n / 2 j = 2 * i / / j从 2 倍i开始 while j < n: x[j] = 0 / / 把所有i的倍数筛除 j = j + i / / 下一个i的倍数 return x / / 返回数组 def prime(n,x): / / 求第n个素数,只需要在筛选好的布尔数组中找第n个标记为 1 的数就可以了 #Find nth prime i = 1 / / 初始化为 1 j = 1 while j < = n: / / 在布尔数组中寻找第n个标记为 1 的数 if x[i] = = 1 : j = j + 1 i = i + 1 return i - 1 / / 前面循环中i多加了一次,返回时需要减 1 |
栏目列表
最新更新
nodejs爬虫
Python正则表达式完全指南
爬取豆瓣Top250图书数据
shp 地图文件批量添加字段
爬虫小试牛刀(爬取学校通知公告)
【python基础】函数-初识函数
【python基础】函数-返回值
HTTP请求:requests模块基础使用必知必会
Python初学者友好丨详解参数传递类型
如何有效管理爬虫流量?
SQL SERVER中递归
2个场景实例讲解GaussDB(DWS)基表统计信息估
常用的 SQL Server 关键字及其含义
动手分析SQL Server中的事务中使用的锁
openGauss内核分析:SQL by pass & 经典执行
一招教你如何高效批量导入与更新数据
天天写SQL,这些神奇的特性你知道吗?
openGauss内核分析:执行计划生成
[IM002]Navicat ODBC驱动器管理器 未发现数据
初入Sql Server 之 存储过程的简单使用
这是目前我见过最好的跨域解决方案!
减少回流与重绘
减少回流与重绘
如何使用KrpanoToolJS在浏览器切图
performance.now() 与 Date.now() 对比
一款纯 JS 实现的轻量化图片编辑器
关于开发 VS Code 插件遇到的 workbench.scm.
前端设计模式——观察者模式
前端设计模式——中介者模式
创建型-原型模式