版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

本帖產生 3 個程序強帖，點擊這里進行查看

holmescn

金蟲 (正式寫手)

程序強帖: 37
應助: 1 (幼兒園)
金幣: 1918.8
散金: 275
紅花: 1
帖子: 699
在線: 102.6小時
蟲號: 913482
注冊: 2009-11-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

[交流] Euler 工程第三十題已有5人參與

又是一個指數(shù)的題啦！

說有3個數(shù)可以寫成各位數(shù)字的4次方的和:

1634 = 1^4 + 6^4 + 3^4 + 4^4
8208 = 8^4 + 2^4 + 0^4 + 8^4
9474 = 9^4 + 4^4 + 7^4 + 4^4

雖然 1 = 1^4, 但這不是一個求和, 所以這個不算.

這三個數(shù)的和為: 1634 + 8208 + 9474 = 19316

那么那些數(shù)可以寫成各位數(shù)字的5次方的和呢? 這些數(shù)的和又是多少?

回復此樓

» 猜你喜歡

08工科 320總分求調劑已經(jīng)有9人回復
生物學071000 329分求調劑已經(jīng)有5人回復
306求調劑已經(jīng)有5人回復
284求調劑已經(jīng)有4人回復
尋找調劑已經(jīng)有4人回復
一志愿中南化學（0703）總分337求調劑已經(jīng)有10人回復
324求調劑已經(jīng)有6人回復
一志愿西北大學，070300化學學碩，總分287，雙非一本，求調劑。已經(jīng)有3人回復
289求調劑已經(jīng)有7人回復
319求調劑已經(jīng)有4人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

Euler 工程第三十八題已經(jīng)有9人回復
Euler 工程第三十六題：已經(jīng)有18人回復
Euler 工程第三十五題：循環(huán)質數(shù) 已經(jīng)有16人回復
Euler 工程第廿九題：有多少不同的項? 已經(jīng)有30人回復
Euler 工程第廿六題：最長的循環(huán)節(jié) 已經(jīng)有9人回復
Euler 工程第廿四題：全排列的第100萬項已經(jīng)有19人回復
Euler 工程第廿三題：已經(jīng)有16人回復
Euler 工程第十八題：三角陣上最大的和已經(jīng)有12人回復
Euler 工程第十六題：2的1000次方的各項和已經(jīng)有14人回復
Euler 工程第十五題：從左上角到右下角有多少條路？已經(jīng)有5人回復
Euler 工程第14題：找最長的數(shù)列已經(jīng)有9人回復
Euler Project Q13 歐拉工程第十三題已經(jīng)有20人回復
Euler Project Q12 歐拉工程第十二題已經(jīng)有23人回復
Euler 工程第十一題：相鄰元素乘積最大已經(jīng)有10人回復
Euler Project Q7. 歐拉工程第七題已經(jīng)有14人回復
Euler 工程第六題：平方和與和的平方差多少？已經(jīng)有5人回復

1樓 2011-06-18 11:22:15

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

holmescn

金蟲 (正式寫手)

程序強帖: 37
應助: 1 (幼兒園)
金幣: 1918.8
散金: 275
紅花: 1
帖子: 699
在線: 102.6小時
蟲號: 913482
注冊: 2009-11-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

★ ★ ★
余澤成(金幣+3, 程序強帖+1): 鼓勵交流！ 2011-06-18 15:52:24

Python

CODE:

# coding: utf-8

r = 0

n = 9**5 * 5

for i in xrange(2, n):

    s = sum([int(x)**5 for x in str(i)])

    if i == s:

        r += i

        print i

print "sum =", r

結果:

引用回帖:

4150
4151
54748
92727
93084
194979
sum = 443839

[ Last edited by holmescn on 2011-6-18 at 14:02 ]

贊一下(1人)

回復此樓

2樓2011-06-18 14:00:45

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

huycwork

金蟲 (著名寫手)

程序強帖: 22
應助: 0 (幼兒園)
金幣: 953
散金: 663
紅花: 8
沙發(fā): 13
帖子: 1080
在線: 264.1小時
蟲號: 1257243
注冊: 2011-04-06
專業(yè): 金融學

★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
余澤成(金幣+3, 程序強帖+1): 鼓勵交流！ 2011-06-18 15:52:37

這個不是水仙花數(shù)嘛~
傳說中的恐怖O(10^n)問題。位數(shù)大一些就要篩數(shù)了。
咋一看，兩邊的解空間是一樣的，大約都是9*10^(n-1)，但是左邊的解空間是緊湊的，右邊的則是松散的，而且，右邊的解空間映射到左邊的范圍不是很大，篩數(shù)從右邊開始，這樣就需要給出一個函數(shù)的上下限，比如要求出最接近100和1000的a^3+b^3+c^3。對于這樣的線性規(guī)劃問題可以在10*10*n的時間內找到最優(yōu)解，然后調用歐拉24題給出的那種以字典序計數(shù)的排列函數(shù)來求解，應該效率會不錯。

[ Last edited by huycwork on 2011-6-18 at 15:29 ]

贊一下(2人)

回復此樓

漩渦的中心有一塊空地，空空的。

3樓2011-06-18 15:22:44

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

libralibra

至尊木蟲 (著名寫手)

驃騎將軍

程序強帖: 40
應助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時
蟲號: 696514
注冊: 2009-02-05
專業(yè): 計算機軟件

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
xzhdty(金幣+2): 鼓勵交流！ 2011-06-18 22:06:27

我也是py

CODE:

 from mytictoc import tic,toc

tic()

# find upper bound

n = 1

while n*(9**5)>10**(n+1)-1:

    n += 1

# compute

print reduce(lambda x,y:x+y,[n for n in xrange(2,n*(9**5)) if n==reduce(lambda x,y:x+y,[int(c)**5 for c in str(n)])])

toc()

效率低下

CODE:

 443839

Elapsed time: 5.64833575 seconds

贊一下(2人)

回復此樓

matlab/VB/python/c++/Java寫程序請發(fā)QQ郵件:790404545@qq.com

4樓2011-06-18 16:31:38

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

holmescn

金蟲 (正式寫手)

程序強帖: 37
應助: 1 (幼兒園)
金幣: 1918.8
散金: 275
紅花: 1
帖子: 699
在線: 102.6小時
蟲號: 913482
注冊: 2009-11-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
dubo(金幣+1): 歡迎常來程序語言版討論 2011-06-19 15:30:25

引用回帖:

Originally posted by libralibra at 2011-06-18 16:31:38:
我也是py
[code] from mytictoc import tic,toc

tic()

# find upper bound
n = 1
while n*(9**5)>10**(n+1)-1:
n += 1

# compute
print reduce(lambda x,y:x+y,[n for n in xrange(2,n*(9** ...

libralibra兄, reduce 比sum 要快嗎?

贊一下(1人)

回復此樓

5樓2011-06-18 18:00:12

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

libralibra

至尊木蟲 (著名寫手)

驃騎將軍

程序強帖: 40
應助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時
蟲號: 696514
注冊: 2009-02-05
專業(yè): 計算機軟件

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
dubo(金幣+1): 歡迎常來程序語言版討論 2011-06-19 15:30:34
微塵、夢想(程序強帖+1): 2011-06-19 17:23:58

引用回帖:

Originally posted by holmescn at 2011-06-18 18:00:12:
libralibra兄, reduce 比sum 要快嗎?

不會吧,我覺得sum應該快吧.我檢測結果也是,還只測試了10^3-10^7的數(shù),要是規(guī)模大,估計lambda表達式更慢.
感覺reduce,map,zip,lambda..就是為了寫出來好看,少寫幾行代碼的.

另,matlab我也發(fā)現(xiàn)這個問題,第28題螺旋矩陣四角求和那個,后面六十幾還是八十幾有個類似的,找到4個數(shù)等差數(shù)列的通項,用sum(start:step:end),并沒有這樣快: start+(start+step)+(start+2*step)+(start+3*step)

CODE:

================================

REDUCE+LAMBDA: 1000

499500

Elapsed time: 0.00649915 seconds

SUM: 1000

499500

Elapsed time: 0.00453717 seconds

================================

REDUCE+LAMBDA: 10000

49995000

Elapsed time: 0.00579124 seconds

SUM: 10000

49995000

Elapsed time: 0.00484922 seconds

================================

REDUCE+LAMBDA: 100000

4999950000

Elapsed time: 0.03130733 seconds

SUM: 100000

4999950000

Elapsed time: 0.00988506 seconds

================================

REDUCE+LAMBDA: 1000000

499999500000

Elapsed time: 0.34010105 seconds

SUM: 1000000

499999500000

Elapsed time: 0.13641119 seconds

================================

REDUCE+LAMBDA: 10000000

49999995000000

Elapsed time: 3.26226713 seconds

SUM: 10000000

49999995000000

Elapsed time: 1.39761869 seconds

代碼

CODE:

from mytictoc import tic,toc

for i in xrange(3,8):

    print '================================'

    print 'REDUCE+LAMBDA: %d' % (10**i)

    tic()

    print reduce(lambda x,y:x+y,xrange(10**i))

    toc()

    print 'SUM: %d' % (10**i)

    tic()

    print sum(xrange(10**i))

    toc()

[ Last edited by libralibra on 2011-6-18 at 22:00 ]

贊一下(3人)

回復此樓

matlab/VB/python/c++/Java寫程序請發(fā)QQ郵件:790404545@qq.com

6樓2011-06-18 21:55:14

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

holmescn

金蟲 (正式寫手)

程序強帖: 37
應助: 1 (幼兒園)
金幣: 1918.8
散金: 275
紅花: 1
帖子: 699
在線: 102.6小時
蟲號: 913482
注冊: 2009-11-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
dubo(金幣+1): 歡迎常來程序語言版討論 2011-06-19 15:30:40

喔，那這最后就可以這么寫了：

CODE:

print sum([n for n in xrange(2, n) if n == sum([int(c)**5 for c in str(n)])])

是吧，

贊一下(1人)

回復此樓

» 本帖已獲得的紅花（最新10朵）

stone19111

7樓2011-06-18 22:51:13

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

qinghuoly

木蟲 (正式寫手)

程序強帖: 1
應助: 140 (高中生)
金幣: 4753.1
紅花: 10
帖子: 885
在線: 918.1小時
蟲號: 966225
注冊: 2010-03-09
性別: GG
專業(yè): 金屬有機化學

★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
dubo(金幣+1): 歡迎常來程序語言版討論 2011-06-19 15:30:47
微塵、夢想(金幣+2): 請使用BBcode代碼，詳見維基百科BBcode 2011-06-19 17:25:25

上我代碼，scheme語言

[define [ans30]
  [define N 1e7]
  [define [d x]
[apply +
         [map [lambda [n] [expt n 5]]
            [map string->number
                  [map string
                        [string->list [number->string x]]]]]]]
  [define [fun n l flag]
[if [> n [add1 N]]
      [if flag
         [cons [- n 1] l]
         l]
      [fun [add1 n]
         [if flag
               [cons [- n 1] l]
               l]
         [= n [d n]]]]]

  [apply + [fun 2 '[] #f]]]
;end of code

;答案：443839

;符合的數(shù)為：(194979 93084 92727 54748 4151 4150)

;N取值為1e7時cpu time: 76078 real time: 77312 gc time: 10089
;N取值為1e6時cpu time: 7453 real time: 7672 gc time: 638

[ Last edited by qinghuoly on 2011-6-19 at 12:14 ]

贊一下(3人)

回復此樓

天地為帳，日月為燈，風雷為號角，云虹為旗令，山川為陣圖，草木為兵卒。運陰陽五行為謀，策古今興替為略。

8樓2011-06-19 12:11:32

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

holmescn

金蟲 (正式寫手)

程序強帖: 37
應助: 1 (幼兒園)
金幣: 1918.8
散金: 275
紅花: 1
帖子: 699
在線: 102.6小時
蟲號: 913482
注冊: 2009-11-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

Originally posted by qinghuoly at 2011-06-19 12:11:32:
上我代碼，scheme語言

[define [ans30]
  [define N 1e7]
  [define [d x]
[apply +
         [map [lambda [n] [expt n 5]]
            [map string->number
                  [ma ...

原來 scheme用的是方括號啊, 真是不習慣啊! 還是Lisp的圓括號好看,哈哈

贊一下

回復此樓

9樓2011-06-19 17:02:07

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

whiterye

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 44.4
紅花: 1
帖子: 11
在線: 4.2小時
蟲號: 1320212
注冊: 2011-06-10
專業(yè): 計算機網(wǎng)絡

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

引用回帖:

Originally posted by holmescn at 2011-06-19 17:02:07:
原來 scheme用的是方括號啊, 真是不習慣啊! 還是Lisp的圓括號好看,哈哈

用Scheme 語言的人好少哦！

贊一下(1人)

回復此樓

10樓2011-06-19 22:43:56

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 holmescn 的主題更新

返回列表

普通表情龍兔虎貓高級回復 (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿北京化工大學070300 學碩336求調劑 +5	vv迷 2026-03-21	8/400	2026-03-22 14:20 by ColorlessPI
[考研] 化學調劑 +5	yzysaa 2026-03-21	5/250	2026-03-21 22:12 by peike
[考研] 化學工程321分求調劑 +18	大米飯！ 2026-03-15	22/1100	2026-03-21 20:20 by HH領袖
[考研] 266求調劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[基金申請] 學校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 0805材料320求調劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 332求調劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 材料學學碩080502 337求調劑-一志愿華中科技大學 +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 070300化學319求調劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 初始318分求調劑（有工作經(jīng)驗） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 華東師范大學-071000生物學-293分-求調劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 一志愿重慶大學085700資源與環(huán)境專碩，總分308求調劑 +3	墨墨漠 2026-03-18	3/150	2026-03-21 00:39 by JourneyLucky
[考研] 材料專業(yè)求調劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 材料學求調劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 生物學調劑招人�。�！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 302求調劑 +4	小賈同學123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學123
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 333求調劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客

24小時熱門版塊排行榜

[交流] Euler 工程 第三十題 已有5人參與

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

» 本帖已獲得的紅花（最新10朵）

[交流] Euler 工程第三十題已有5人參與

» 本主題相關價值貼推薦，對您同樣有幫助: