Re: [問題] 質數_巢狀迴圈_菲絲恩

作者: APM99 (血統純正台北人)   2017-08-10 12:50:28
: i=j=1
: for i in range(2,100,1):
: for j in range(2,int(i/j)+1):
: if(not i%j):
: break
: if j>i**0.5:
: print('%d is prime'%(i))
這裡用的數學方法是大家高中都學過的
Sieve of Eratosthenes(wiki中文:埃拉托斯特尼篩法)
其中的一個引理
想判斷 i = 29 能不能被 7 整除
只需要判斷4項 ,即 29/2 29/3 29/4 29/5 能不能整除
驗證 29/6 已經沒意義了,因為 7*5 已經大於29了
i/j是一種等分的概念.
#回到python虛擬碼
i=j=1
for i in range(2,100,1):
for j in range(2,(無條件進位到整數的i/j) +1):
if(not i%j):
break
if j>i**0.5:
print('%d is prime'%(i))
#
1. 灰色 +1 是python range特性
2. 淺藍色可能只是想打出2,不然跟等分法並不連貫~,
3. 紅色部分
一般都是取巧用四捨五入法,結果python沒印出 5..
才發現
python3 中 round(2.5) = 2
python2 中 round(2.5) = 3
可見
https://stackoverflow.com/questions/10825926/python-3-x-rounding-behavior
作者: APM99 (血統純正台北人)   2017-08-10 12:52:00
py3用四捨五入法得要 round(i/j + 0.0000000000001) 才行
作者: CaptainH (Cannon)   2017-08-10 12:59:00
你是不是誤會了什麼 要判斷29是不是被7整除 只要他媽的除一下就好四捨入+0.5好就 加一個0.000...1不知道搞什麼鬼或者用ceil/floor 根本沒這問題
作者: bruce0209 (士賢)   2017-08-10 13:12:00
要判斷29是不是被7整除 我也思考了很久在說啥…
作者: Django (Cython)   2017-08-10 15:47:00
?????????
作者: nknuukyo (我無所能因敵成體)   2017-08-10 16:00:00
謝謝~wiki裡有一段python3.6的引碼,不曉得是否方便說明對應到原程式的關係^^"
作者: bruce0209 (士賢)   2017-08-10 16:20:00
樓上是說sieve of Eratosthenes的wiki嗎?如果是的話 wiki裡面那個是用刪除法的 和這邊不太一樣你可以看他wiki裡面附的jpg應該就知道他在做什麼了是說...有程式碼了怎麼不自己trace看看
作者: stucode   2017-08-10 19:34:00
... 先不論原本的lemma是什麼 你的code真的是問題滿載
作者: bruce0209 (士賢)   2017-08-10 19:37:00
postive?????? 豆頁女子痛...推薦可以看一下clear code相關書.......
作者: stucode   2017-08-10 19:44:00
isprime就return n 那while isprime:是?每次call一開始j都是1 那n/j是在表現什麼?
作者: CaptainH (Cannon)   2017-08-10 20:43:00
愈寫愈錯XDD 平鋪直敘的算法也能寫這麼醜XDD

Links booklink

Contact Us: admin [ a t ] ucptt.com