版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

本帖產(chǎn)生 4 個程序強帖，點擊這里進行查看

holmescn

金蟲 (正式寫手)

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

[交流] Euler 工程第14題：找最長的數(shù)列已有6人參與

周末了，放個題出來玩玩：

定義一個正整數(shù)數(shù)列，其迭代公式為：
n = n/2 (當(dāng)n為偶數(shù))
n = 3n+1 (當(dāng)n為奇數(shù))

比如從n=13開始，計算這個數(shù)列得：
13 ->40->20->10->5->16->8->4->1
這個數(shù)列一共有10項。
這個數(shù)列是不是總收斂到1還是個沒有解決的問題（稱為Collatz Problem)
不過，我們并不是要解這個難題，而是要求在小于1百萬的所有起始數(shù)中，哪個數(shù)能產(chǎn)生最長的數(shù)列。

這里要注意的是，數(shù)列中間的項是可以大于1百萬的，要看提最后數(shù)列終止到1時候的長度。

PS：這東西不能用歸納法證出來嗎？
PS2：注意1分鐘原則啊。不過，大于1分鐘的程序也可以發(fā)上來，這樣才有的交流。

Here we go!

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有7人回復(fù)
272求調(diào)劑已經(jīng)有7人回復(fù)
352分化工與材料已經(jīng)有3人回復(fù)
0703化學(xué) 已經(jīng)有6人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有3人回復(fù)
266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑已經(jīng)有4人回復(fù)
一志愿南師大0703化學(xué) 275求調(diào)劑已經(jīng)有4人回復(fù)
085701環(huán)境工程求調(diào)劑已經(jīng)有3人回復(fù)
266求調(diào)劑已經(jīng)有10人回復(fù)

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

Euler 工程第四十四題已經(jīng)有4人回復(fù)
Euler 工程第四十二題: 三角詞已經(jīng)有4人回復(fù)
Euler 工程第四十一題已經(jīng)有5人回復(fù)
Euler 工程第三十八題已經(jīng)有9人回復(fù)
Euler 工程第三十七題已經(jīng)有6人回復(fù)
Euler 工程第三十六題：已經(jīng)有18人回復(fù)
Euler 工程第三十五題：循環(huán)質(zhì)數(shù) 已經(jīng)有16人回復(fù)
Euler 工程第三十二題：pandigital 數(shù) 已經(jīng)有3人回復(fù)
Euler 工程第三十一題: 換零錢已經(jīng)有10人回復(fù)
Euler 工程第三十題已經(jīng)有12人回復(fù)
Euler 工程第廿九題：有多少不同的項? 已經(jīng)有30人回復(fù)
Euler 工程第廿八題：旋轉(zhuǎn)矩陣對角線的和已經(jīng)有6人回復(fù)
Euler 工程第廿七題：系數(shù)的積已經(jīng)有15人回復(fù)
Euler 工程第廿六題：最長的循環(huán)節(jié) 已經(jīng)有9人回復(fù)
Euler 工程第廿五題：Fibonacci 數(shù)列第一個包含1000個數(shù)字的項已經(jīng)有3人回復(fù)
Euler 工程第廿四題：全排列的第100萬項已經(jīng)有19人回復(fù)
Euler 工程第廿三題：已經(jīng)有16人回復(fù)
Euler 工程第廿二題: 姓的總分已經(jīng)有13人回復(fù)
Euler 工程第廿題：100! 的各項和已經(jīng)有5人回復(fù)
Euler 工程第十九題：每月第一天是周日的天數(shù) 已經(jīng)有4人回復(fù)
Euler 工程第十八題：三角陣上最大的和已經(jīng)有12人回復(fù)
Euler 工程第十六題：2的1000次方的各項和已經(jīng)有14人回復(fù)
Euler 工程第十一題：相鄰元素乘積最大已經(jīng)有10人回復(fù)
Euler 工程第三題：尋找600851475143的最大質(zhì)因子已經(jīng)有18人回復(fù)

1樓 2011-05-20 16:41:24

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

huycwork

金蟲 (著名寫手)

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

★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
余澤成(金幣+3, 程序強帖+1): 辛苦了！ 2011-05-20 21:02:20

昨天寫好這個了，C++版本的是醬紫：

CODE:

#include 

enum {BUFSZ = 1000000};

size_t eular14(size_t bufsz = 1000000){

        size_t buf[BUFSZ+1];

        size_t max = 1, c, r;

        long long n;

        buf[1] = 1;

        for(int i = 2; i < BUFSZ + 1; ++i){

                n = i;

                c ^= c;

                while(1){

                        if(n < BUFSZ && buf[n])

                                break;

                        ++c;

                        if(n % 2){

                                n = 3*n + 1;

                        }else

                                n = n/2;

                }

                buf[i] = buf[n] + c;

                if(buf[i] > max){

                        max = buf[i];

                        r = i;

                }

        }

        return r;

}

int main(void){

        std::cout<
}

贊一下(2人)

回復(fù)此樓

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

2樓2011-05-20 17:05:54

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

wangww2011

木蟲 (著名寫手)

程序強帖: 13
應(yīng)助: 11 (小學(xué)生)
金幣: 4023.1
散金: 2709
紅花: 18
沙發(fā): 1
帖子: 1915
在線: 1537.1小時
蟲號: 772953
注冊: 2009-05-17
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
余澤成(金幣+1): 說到！ 2011-05-20 21:02:42
微塵、夢想(程序強帖+1): 很好，歡迎常來！ 2011-05-21 19:18:45

給個迭代算法吧，0.05s 先看結(jié)果：

CODE:

max=837799  count=525

elapsed time=0.050000 seconds.

代碼為

CODE:

#include 

#include 

#include 

#define TIMERSTART clock_t start_time,stop_time;double elapsed_time;start_time = clock();

#define TIMERSTOP stop_time = clock();elapsed_time=(double)(stop_time-start_time)/CLOCKS_PER_SEC;printf("elapsed time=%f seconds.\n",elapsed_time);

#define N 1000001

static int count[N];

int euler14(long long n){

  int result;

  if(n0){

    return count[n];

  }

  if(n%2) {

    n=3*n+1;

  }else {

    n=n/2;

  }

  result=euler14(n);

  if(n
    count[n]=result;

  } 

  return result+1;

}

int main(void){

  int i=0;

  int max_count,max;

  TIMERSTART;

  count[1]=1;

  max_count=0;

  for(i=N-1;i>1;i--){

    if(count[i]==0){

       count[i]=euler14(i);

    }

    if(count[i]>max_count){

      max_count=count[i];

      max=i;

    }

  }

  printf("max=%d  count=%d\n",max,max_count);

  TIMERSTOP;

  return 0;

}

[ Last edited by wangww2011 on 2011-5-21 at 14:14 ]

贊一下(3人)

回復(fù)此樓

3樓2011-05-20 17:53:37

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

huycwork

金蟲 (著名寫手)

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

★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
余澤成(金幣+3): 鼓勵交流！ 2011-05-20 21:03:01

這個用歸納法證明是證明不出來的
前提條件可以選2->1
但是后面無法作出歸納假設(shè)：就算你假設(shè)an -> ak，歸納的下一步驟就是a(n+1)沒法推導(dǎo)a(k+1)，你看代碼就知道了，你無法知道中間產(chǎn)生的鏈條有多長，實際原因是無法知道中間經(jīng)過的元素有哪些，當(dāng)然也就無法歸納證明。

贊一下(2人)

回復(fù)此樓

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

4樓2011-05-20 18:03:01

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

微塵、夢想

木蟲 (知名作家)

程序強帖: 6
應(yīng)助: 2 (幼兒園)
貴賓: 0.353
金幣: 4757.9
散金: 3089
紅花: 31
沙發(fā): 247
帖子: 8788
在線: 1125小時
蟲號: 1203290
注冊: 2011-02-14
專業(yè): 制造系統(tǒng)與自動化

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
余澤成(金幣+1): 鼓勵交流！ 2011-05-20 21:03:21

這不是著名的3x+1問題嗎？以前寫過這個程序，也試著證過這道題，看著挺簡單，要真證起來，才發(fā)現(xiàn)根本不可能，否則它也不是世界難題了……

贊一下(2人)

回復(fù)此樓

任風(fēng)云變幻，我笑對人生！

5樓2011-05-20 18:27:20

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

libralibra

至尊木蟲 (著名寫手)

驃騎將軍

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

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
jjdg(金幣+2): 辛苦了 2011-05-21 00:51:59
微塵、夢想(程序強帖+1): 2011-05-21 19:19:39

matlab的,13s多哦

CODE:

 %% Which starting number, under one million, produces the longest chain?

%   n = n/2 (even); n = 3*n+1 (odd) till n==1

% like starting from 13, we generate the following sequence (contains 10 terms)

%  13 → 40 → 20 → 10 → 5 → 16 → 8 → 4 → 2 → 1

function result = euler14()

tic;

n = 0;

len = 0;

big = 0;

result = 0;

for i=1:1e6

    n = i;

    len = 0;

    while n~=1

        if mod(n,2)==0

            n = n/2;

        else

            n = n*3+1;

        end

        len = len+1;

    end

    len = len+1; % n==1

    if len>big

        big = len;

        result = i;

    end

end

toc;

end

結(jié)果

CODE:

 % Elapsed time is 13.644453 seconds.

% ans =

%       837799

贊一下(3人)

回復(fù)此樓

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

6樓2011-05-20 21:23:11

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

holmescn

金蟲 (正式寫手)

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

微塵、夢想(程序強帖+1): 2011-05-21 19:19:52

matlab果然是不給力啊，F(xiàn)ortran版的只要0.75s

CODE:

Program euler14

    Implicit None

    Integer :: maxLen = 0, l

    Integer(8) :: n

    Integer :: start, maxStart

    Real :: t1, t2

    Call CPU_Time(t1)

    Do start = 2, 1000000

        n = start

        l = 0

        Do While (n .gt. 1) 

            If(mod(n, 2) .eq. 0) Then

                n = n / 2

            Else

                n = 3*n + 1

            EndIf

            l = l + 1

        EndDo

        l = l + 1

        If(l .gt. maxLen) Then

            maxLen = l

            maxStart = start

        EndIf

    EndDo

    Call CPU_Time(t2)

    Print *, "Result:", maxStart, maxLen, t2-t1

End Program euler14

結(jié)果同libralibra同學(xué)。不過恐怖的是，中間項居然會大于40億啊。32位整數(shù)都保存不了了。最后總共有525項。這東西用Mathematica的NestList展開會是什么樣？

現(xiàn)在的計算機速度太快了，以至于我們都不想去找一個更好的解法了，不過，這東西的解空間應(yīng)該還能減小，我看到2到3000以內(nèi)都還是單調(diào)上漲的呢。

贊一下(2人)

回復(fù)此樓

7樓2011-05-21 13:40:36

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

wmc_1979

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1232.8
散金: 30
帖子: 128
在線: 110.9小時
蟲號: 1243434
注冊: 2011-03-24
專業(yè): 運籌學(xué)

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

我計算出是 837799
項數(shù) 524

[ Last edited by wmc_1979 on 2011-5-21 at 20:38 ]

贊一下(1人)

回復(fù)此樓

8樓2011-05-21 20:04:41

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

holmescn

金蟲 (正式寫手)

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

★
余澤成(金幣+1): 謝謝參與交流！ 2011-05-22 11:53:05

引用回帖:

Originally posted by wmc_1979 at 2011-05-21 20:04:41:
我計算出是 837799
項數(shù) 524

[ Last edited by wmc_1979 on 2011-5-21 at 20:38 ]

當(dāng)n=1的時候，就退出循環(huán)了，所以認(rèn)為最后計數(shù)的時候沒有數(shù)1

贊一下(1人)

回復(fù)此樓

9樓2011-05-22 09:08:54

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

chyanog

金蟲 (小有名氣)

程序強帖: 1
應(yīng)助: 12 (小學(xué)生)
金幣: 808.1
紅花: 5
帖子: 85
在線: 52.8小時
蟲號: 1540955
注冊: 2011-12-17
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算

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

Mathematica和中的遞歸可以很方便的“記憶化”，很簡潔但還是沒有編譯后的循環(huán)快：

CODE:

c@1 = 1; c@n_ := c@n = 1 + If[EvenQ@n, c[n/2], c[3 n + 1]]; 

Ordering[c /@ Range@1*^6, -1] // AbsoluteTiming

cf = Compile[{{n0, _Integer}},

   Block[{n = N@n0},

    Catch@Do[If[n == 1, Throw@i];

      n = If[FractionalPart[n/2] == 0, n/2, 3*n + 1],

      {i, 1*^8}]

    ], RuntimeAttributes -> Listable,

   CompilationTarget -> "C", RuntimeOptions -> "Speed"

   ];

Ordering[cf@Range@1*^6, -1] // AbsoluteTiming

贊一下

回復(fù)此樓

10樓2013-10-05 01:49:15

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 holmescn 的主題更新

返回列表

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

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

[考研] 272求調(diào)劑 +7	腳滑的守法公民 2026-03-27	7/350	2026-03-27 17:23 by laoshidan
[考研] 340求調(diào)劑 +4	jhx777 2026-03-27	4/200	2026-03-27 16:50 by sanrepian
[考博] 26申博 +3	加油沖��！ 2026-03-26	3/150	2026-03-27 15:38 by cls512
[考研] 085600，材料與化工321分，求調(diào)劑 +9	大饞小子 2026-03-27	9/450	2026-03-27 14:30 by mmm just
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +3	Zuhui0306 2026-03-25	3/150	2026-03-27 14:09 by shangxh
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 327求調(diào)劑 +7	prayer13 2026-03-23	7/350	2026-03-26 20:48 by 不吃魚的貓
[考研] 求調(diào)劑 +8	Auroracx 2026-03-22	8/400	2026-03-26 19:55 by 不吃魚的貓
[考研] 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +8	y7czhao 2026-03-26	8/400	2026-03-26 19:51 by 不吃魚的貓
[考研] 材料277求調(diào)劑 +5	min3 2026-03-24	5/250	2026-03-26 15:13 by zzll406
[考研] 085600 材料與化工 329分求調(diào)劑 +9	Mr. Z 2026-03-25	9/450	2026-03-26 10:36 by baoball
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 303求調(diào)劑 +6	藍山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8	皮皮皮123456 2026-03-22	8/400	2026-03-23 23:47 by 熱情沙漠
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO

24小時熱門版塊排行榜

[交流] Euler 工程 第14題：找最長的數(shù)列 已有6人參與

» 猜你喜歡

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

[交流] Euler 工程第14題：找最長的數(shù)列已有6人參與