作者WingedDragon (黄金会死鸟-死後无法复活)
看板Python
标题[问题] 质数列表
时间Thu Jun 9 15:56:53 2016
下面程式是用来产生质数表
我大概知道是用 埃拉托斯特尼筛法
不过有些部分看不懂
--
from time import time
t1 = time()
def primes(n):
p = [true] * (n/2)
for i in xrange(3, int(n**0.5)+1, 2):
if p[i/2]:
p[i*i/2::i] = [False] * ((n-i*i-1)/(2*i)+1)
return [2] + [2*i+1 for i in xrange(1, n/2) if p[i]]
print len( primes(10**7) )
print time()-t1
--
我的问题主要有以下:
1. i/2 在奇数应该会变成浮点数, 为何 p[i/2] 不会错误
2. 黄色部分看不懂
来源:
http://tieba.baidu.com/p/1746377541 9楼回答
--
历代主角: 武藤
游戏---神抽
游城十代---强运 不动
游星---印卡 九十九
游马---搓牌
翼神龙 效果:
此卡不可特殊召唤...
神兽王 表示:同样三祭品 我免费炸半场外加三千打点
裁龙 表示:同样支一千 我能炸全场还不用扣血加攻
巨神兵 表示:听说我可以特召
天空龙 表示:我现在可以捏死原作狂特召的你
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.25.105
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Python/M.1465459017.A.285.html
※ 编辑: WingedDragon (140.112.25.105), 06/09/2016 15:57:53
1F:→ bigpigbigpig: i/2 在 Python 2 中是整数除法,5/2 = 2 而非 2.5 06/09 16:23
2F:→ bigpigbigpig: p与xrange的对映关系研究一下,黄色那行就是在实作 06/09 16:48
3F:→ bigpigbigpig: Sieve of Eratosthenes,把 i 的倍数全部删除 06/09 16:49
我主要是语法看不懂, python 3 刚接触, 很多语法和类 C 语言不同
p[i*i/2::i] 这句完全不知道是什麽意思
[False] * ((n-i*i-1)/(2*i)+1)
我只知道是要做出 ((n-i*i-1)/(2*i)+1) 数量的 False
我知道 埃式筛 的操作过程
不过语法看不懂, 所以黄色那句才完全看不懂
※ 编辑: WingedDragon (140.112.25.105), 06/11/2016 20:14:38
4F:→ bigpigbigpig: p[i*i/2::i]→从i*i/2到结尾,每隔i个间隔取值,即 06/11 23:55
5F:→ bigpigbigpig: index为i*i/2(=a),a+i,a+2*i,...,直到p的结尾, 06/11 23:56
6F:→ bigpigbigpig: 个人觉得这样的最佳化有点过头,找序列的slice算子 06/11 23:58
那 [False] * ((n-i*i-1)/(2*i)+1) 是如何和那些数字对应 ?
是一个 i 就产生一次, 还是每一个间隔产生一次 ?
※ 编辑: WingedDragon (140.112.4.192), 06/12/2016 21:10:16
7F:→ zerof: [False] 会把筛出来的位置都换成 False, 要下一次筛到 True 06/15 02:22
8F:→ zerof: 才会再执行 [False] * ____ 06/15 02:23