[中譯] ProjectEuler 362 Squarefree factors

看板puzzle作者 (企鵝靈)時間14年前 (2011/12/19 19:31), 編輯推噓3(300)
留言3則, 3人參與, 最新討論串1/1
※ [本文轉錄自 chyrliin 信箱] 作者: babufong (嗶嗶) 看板: puzzle 標題: [中譯] ProjectEuler 362 Squarefree factors 時間: Sun Dec 11 15:30:32 2011 362. Squarefree factors http://projecteuler.net/problem=362 來看看 54 這個數字 54 可以拆成 7 種全然不同的一個或多個因數相乘的組成方法,而這些因數都大於 1 如:54 , 2 X 27 , 3 X 18 , 6 X 9 , 3 X 3 X 6 , 2 X 3 X 9 和 2 X 3 X 3 X 3 如果我們要求這些因數都要是「無平方數因數」的數所組成,那就只剩下 2 種組法: 3 X 3 X 6 和 2 X 3 X 3 X 3 Fsf(n) 為 n 有幾種上述條件的方法去組成,所以 Fsf(54) = 2 S(n) 為 ΣFsf(k) , k = 2 到 n 已知 S(100) = 193 請求出 S(10000000000) -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 125.224.9.90

12/11 16:55, , 1F
好酷
12/11 16:55, 1F

12/12 11:54, , 2F
OEIS有數列 (http://oeis.org/A050320 ) 問題在於 10^10...
12/12 11:54, 2F

12/17 19:49, , 3F
如果要跑暴力破解法 應該不會那麼難
12/17 19:49, 3F

12/17 19:49, , 4F
不過要跑一分鐘以內 就有點難度
12/17 19:49, 4F

12/17 19:51, , 5F
OEIS中已經有說同樣的prime signature 其a(n)是一樣的
12/17 19:51, 5F

12/17 19:52, , 6F
所以不必每一個數都去算其a(n), 同樣的prime signature
12/17 19:52, 6F

12/17 19:54, , 7F
算一就好, 並存到表格內(應該是Hash),遇到同樣的prime
12/17 19:54, 7F

12/17 19:54, , 8F
signature, 再從hash table取出即可
12/17 19:54, 8F

12/17 19:57, , 9F
素數的枚舉,也不用枚舉到10^10即可, 枚舉到10^8即可
12/17 19:57, 9F

12/17 19:58, , 10F
大如10^8的素數,假若為p, p*1~p*100其a(n)等於101~101*100
12/17 19:58, 10F

12/17 20:00, , 11F
所以只要算出小於(10^10)/1,(10^10)/2,...,(10^10)/100
12/17 20:00, 11F

12/17 20:00, , 12F
的素數個數即可
12/17 20:00, 12F

12/17 20:03, , 13F
我的方法, 如果是用C++加上現在主流的PC(core i7)
12/17 20:03, 13F

12/17 20:04, , 14F
應該是十分鐘以內可以跑出結果
12/17 20:04, 14F
※ Deleted by: puzzlez (123.194.45.3) 12/18/2011 20:31:07 ※ 發信站: 批踢踢實業坊(ptt.cc) ※ 轉錄者: chyrliin (180.176.101.150), 時間: 12/19/2011 19:31:11

12/19 19:41, , 15F
感謝企鵝靈解救文章!
12/19 19:41, 15F

12/19 20:16, , 16F
結果真的是帕索不小心刪掉的...XD
12/19 20:16, 16F

12/19 20:35, , 17F
= = 不過為什麼企鵝靈找得到文章呢?我都找不到.....
12/19 20:35, 17F
文章代碼(AID): #1Exo212N (puzzle)