| 5 | 1/1 | 返回列表 |
| 查看: 1038 | 回復(fù): 5 | |||
| 本帖產(chǎn)生 2 個 程序強帖 ,點擊這里進行查看 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[交流]
Project Euler 49 歐拉工程 49 題
|
|||
|
以3330為公差的等差數(shù)列1487, 4817, 8147在兩個方面比較特殊: (1)每一項都是四位數(shù)的素數(shù) (2)任一項都可以通過其他項再排列得到 已知沒有一位,兩位或者三位數(shù)的三個素數(shù)能夠展現(xiàn)出上述性質(zhì),但是還有一個由4位素數(shù)組成的數(shù)列滿足上述性質(zhì)。 請問把這個數(shù)列中的三個數(shù)依次連接組成的12位的數(shù)字是多少? |
» 搶金幣啦!回帖就可以得到:
+1/263
+3/164
+1/127
+1/87
+1/83
+1/78
+1/76
+2/68
+1/38
+1/38
+1/35
+1/20
+1/17
+1/17
+1/11
+1/6
+2/6
+1/6
+1/2
+1/2
|
python # -*- coding: cp936 -*- #Project Euler 49 歐拉工程 49 題 #等差3330數(shù)列 #(1)每一項都是四位數(shù)的素數(shù) #(2)任一項都可以通過其他項再排列得到 #由4位素數(shù)組成的三個數(shù) from math import sqrt def isprime(p): #驗證素數(shù),素數(shù)則返回素數(shù)本身,合數(shù)則返回False k=1 for i in xrange(2,int(sqrt(p))+1): if p%i==0: k=0 return False break if k: return p for x in xrange(1001,3338): #9999-6660=3339,x為三個數(shù)中最小的,不大于此 if isprime(x) and isprime(x+3330) and isprime(x+6660) and set(str(x))==set(str(x+3330))==set(str(x+6660)): print x,x+3330,x+6660 1487 4817 8147 2969 6299 9629 [ Last edited by tieer on 2011-9-9 at 11:41 ] |
|
結(jié)果 ['148748178147', '296962999629'] 代碼 #!/usr/bin/env python def generatePrimes(n):#generate all prime numbers less than a given integer n, just take 0.6 seconds for the case n equals one million isprimes=[True]*n for i in range(2,n): if isprimes[i]: for j in range(2*i,n,i): isprimes[j] = False primes=[i for i in range(3,n,2) if isprimes[i]] primes.insert(0,2) return primes def euler49(): p1=generatePrimes(10000) p2=[i for i in p1 if i>1000] res=[] for i in p2: for j in p2: if j>i and set(str(i))==set(str(j)): k=2*j-i if set(str(k))==set(str(i)) and k in p2: res.append(''.join([str(i),str(j),str(k)])) return res if __name__ == "__main__": print euler49() PS:樓上怎么知道公差是3330呢? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程322分 +7 | 哈哈哈吼吼吼哈 2026-04-01 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 求生物學調(diào)劑 +8 | 15172915737 2026-04-01 | 8/400 |
|
|
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +10 | GENJIOW 2026-03-31 | 11/550 |
|
|
[考研] 332求調(diào)劑 +8 | Lyy930824@ 2026-03-29 | 8/400 |
|
|
[考研] 26考研調(diào)劑 +3 | Wnz.20030617 2026-04-01 | 3/150 |
|
|
[考研] 324求調(diào)劑 +10 | hanamiko 2026-03-26 | 12/600 |
|
|
[考研] 材料調(diào)劑 +11 | 一樣YWY 2026-03-31 | 11/550 |
|
|
[考研] 考研調(diào)劑 +11 | Amber00 2026-03-31 | 11/550 |
|
|
[考研] 086000生物與醫(yī)藥298調(diào)劑求助 +4 | 元元青青 2026-03-31 | 6/300 |
|
|
[考研] 318求調(diào)劑 +8 | 七憶77 2026-04-01 | 8/400 |
|
|
[考研] 0855機械初試280求調(diào)劑 +3 | kazenotori 2026-03-31 | 3/150 |
|
|
[考研] 0710生物學求調(diào)劑! +6 | 敘述文 2026-03-31 | 6/300 |
|
|
[考研] 一志愿南師大0703化學 275求調(diào)劑 +5 | Ripcord上岸 2026-03-27 | 5/250 |
|
|
[考研] 一志愿大連理工大學,機械工程學碩,341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-31 | 4/200 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|