C_and_CPP 板


LINE

最近有個熱門的討論話題 就是計算費氏數列的複雜度到底是 O(1) 還是 O(n) 剛好我前幾天在看 wiki 嘗試 compiler 的一些東西的時候 https://zh.wikipedia.org/wiki/%E5%B0%BE%E8%B0%83%E7%94%A8 也遇到一些有趣的 O(1) 還是 O(n) 的問題 覺得很有趣所以就分享上來 我也有把問題丟在 stackoverflow 上面問 沒想到上面的反應也蠻熱烈的 https://stackoverflow.com/questions/54686395 讓我不小心賺到了一些 reputation,大概比我回答十個問題還多 --- 不管複雜度是 O(1) 或是 O(n),也不管 lookup table 到底要不要存在月球上 費氏數列遞迴的形式都是這條:F(x) = F(x-1) + F(x-2), F(1) = F(2) = 1 不過今天要討論的是一個更簡單的問題 F(x) = F(x-1)+1, F(0) = 0 學過中學以上歸納法的人類都能知道 F(x) = x 而學過 C++ 的人都可以把這個式子轉換變成程式碼 int Identity(int i) { if (i == 1) return 1; else return Identity(i-1)+1; } 上面的 wiki 說明了這個並不是有效的 tail recursion 形式 理論上應該不會變成 for loop 會產生 O(n) memory, O(n) runtime 的程式 (PS: 如果是 for loop,應該是 O(1), memory O(n) runtime) 為了驗證,我用了 gcc 8.2.1 編譯看看,結果大出意外 % g++ a.cpp -c -O2 && objdump -d a.o Disassembly of section .text: 0000000000000000 <_Z8Identityi>: 0: 89 f8 mov %edi,%eax 2: c3 ret Linux 下 x64 的第一個整數是 %edi,回傳值放在 %eax 等等,所以我以為 O(n) 的問題,真的可以在 O(1) 解出來嗎(誤) 難道 compiler 做了一些數學計算,推出 F(x) = x 了嗎 該不會在這個大 AI 時代,compiler 也要內建高中生等級的 super AI 了吧 該不會我下次升級到 gcc 9 的時候,我的 compiler 就會跑去當 Youtuber 了吧?! --- Sorry 扯遠了,回歸正題 我一開始的想法是 1. gcc 知道了 negative i 會撞到 UB,因此可以隨便回傳任何值 (PS: negative i 不會變成 infinite loop 的 UB,而是 overflow) 2. positive i 的情形 gcc 經過了某些數學推導算出 F(x) = x 但是怎麼看都覺得太神奇,感覺不會有人實做這種東西 總之經過 stackoverflow 一番討論之後 看起來的結論如下 首先,當代的 gcc 不只可以化簡基本的 tail recursion 就連上面那個形式都可以(可以去看 stackoverflow 的一些討論) 雖然詳細上原理我不太明白,但是應該、大約、好像會變成這樣子 int Identity(int i) { int ans = 0; for (; i != 0; i--, ans++); return ans; } 接著我猜測這個形式對 compiler 來說應該比較好化簡了 因為這種 for loop 非常常見,應該有機會做某種化簡得到 F(x) = x 順帶一題,直接寫這個程式的話,gcc 是可以化簡 O(n) -> O(1) 的 如果 i != 0 改成 i >= 0 的話,gcc 會變成 return i > 0 ? i : 0; 真的很厲害 --- 另外 stackoverflow 裡面有人直接挖出 gcc code 來解釋 但是其實我不是 compiler 專家 所以我這篇主要還是單純分享一些我的觀察啦 如果這個版有人能做出更淺顯易懂,又更完整的解釋的話就太好了 謝謝大家~ -- Time waits for no one. ↑ (。A。)ハァ --



※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 118.160.89.176
※ 文章網址: https://webptt.com/m.aspx?n=bbs/C_and_CPP/M.1550249453.A.382.html
1F:推 KanzakiHAria: 你應該沒搞懂那個討論 02/16 06:33
2F:→ KanzakiHAria: O(1)是說費式數列有一個公式解 02/16 06:33
3F:→ KanzakiHAria: 可是裡面有開根號 所以實務上並不是O(1) 02/16 06:34
4F:→ KanzakiHAria: 開根號速度跟數字長度有關係 02/16 06:34
5F:→ KanzakiHAria: 那個作者非常智障的嗆人time it 實際上就不是O(1) 02/16 06:35
6F:→ KanzakiHAria: 那個人履歷蠻漂亮的 電機出身+待過微軟開發過VS前期 02/16 06:36
7F:→ KanzakiHAria: 看他講話好像連基本的計算機原理和演算法數學都不懂 02/16 06:36
8F:→ KanzakiHAria: 連我以前當助教的學生都可以電爆他了XD 02/16 06:37
9F:推 KanzakiHAria: 講錯了 不是開根號 是次方問題 02/16 06:41
10F:推 KanzakiHAria: 然後O(n)裡的n 一種是編碼長度 一種是input數量 02/16 06:47
11F:→ KanzakiHAria: 因為是編碼長度問題 所以實際上是O(lgn) 02/16 06:49
12F:推 KanzakiHAria: 不過說不定原作者是想表達C++有編譯時期運算技術 02/16 06:52
13F:→ KanzakiHAria: 所以不管n多大C++都會在編譯時期算好 02/16 06:52
14F:→ KanzakiHAria: 所以run-time是O(1)wwwwwwwwww 02/16 06:52
15F:推 KanzakiHAria: https://i.imgur.com/N6tFX0a.png 02/16 07:02
16F:推 alan23273850: 神奇給推! 02/16 08:57
17F:→ KanzakiHAria: 已經跟地球月球時間無關 而是根本上錯誤 02/16 08:57
18F:推 steve8625: 真有趣~~ 02/16 11:55
19F:推 b98901056: Vtuber compiler XD 02/16 12:44
20F:推 firejox: 以編碼長度來看,也不會是O(lg n) 02/16 20:02
21F:→ Sanvean: 我是不是錯過了什麼新聞? 求費氏數列的討論串 pAq 02/17 14:45
22F:推 Higana: https://i.imgur.com/nZmquHA.jpg 02/17 16:26
23F:→ Ommm5566: ㄟ 這是兩件事 1. O(logN) 是因為乘法有快速乘法logN 02/17 17:14
24F:→ Ommm5566: 2. turing machine來看編碼長度確實是logN 02/17 17:16
25F:→ Ommm5566: 然後巧合的是剛好這兩件事可以掛勾在一起 02/17 17:17
26F:推 Ommm5566: 這個討論串居很無聊,居然這麼多人關注。 02/17 17:20
27F:推 LPH66: 樓上的兩個 logN 是不同 N 吧? 02/17 20:02
28F:→ LPH66: 你的"快速乘法"的 log N 的 N 是編碼長度 02/17 20:02
29F:→ LPH66: turing machine 編碼的 logN 的 N 是數值本身 02/17 20:03
30F:推 Ommm5566: Fabonacci(X) 這個是編碼長度logX 所以放在tap上是logX 02/17 20:19
31F:→ Ommm5566: 然後公式解是一個const連乘X次 02/17 20:20
32F:→ Ommm5566: 因為有快速乘法所以時間是logX 02/17 20:21
33F:→ Ommm5566: 這題只是剛好快速乘法的行為跟二進為編碼直接有相關 02/17 20:21
34F:推 Ommm5566: 1000 是四位數編碼 100是三位數編碼 10是兩位數編碼 02/17 20:24
35F:→ Ommm5566: 放在tap上長度本來就是logN 02/17 20:24
36F:推 KanzakiHAria: 所以正確來說編碼長度N的費氏數列時間複雜度O(N) 02/17 21:32
37F:推 KanzakiHAria: 前面可能我表達不好 這個應該沒爭議了吧 02/17 22:36
38F:推 AstralBrain: O(N)次乘法, 但是乘法的時間不一定是O(1), 看你要怎 02/18 02:33
39F:→ AstralBrain: 麼算 02/18 02:33
40F:→ AstralBrain: 算出沒誤差的精確值至少是O(2^N), 因為答案本身就有 02/18 02:34
41F:→ AstralBrain: O(2^N)bit了 02/18 02:34
42F:推 AstralBrain: 啊..修正一下, 底不確定是不是2, 但總之是指數級成長 02/18 02:53
43F:→ gus2: 怎麼推文都在回討論串?本文重心是編譯器優化遞迴f(i)成i吧 02/18 03:13
44F:推 gus2: 有趣給推 02/18 03:15
45F:→ Caesar08: 請問Higana,這是在哪個社團,那麼有趣 XD 02/18 13:04
46F:推 asd456fgh778: 樓上 Python台灣 02/18 14:42
47F:推 Higana: 是的 但該篇他老早就刪了 剩一些後續討論這樣 02/20 01:44







like.gif 您可能會有興趣的文章
icon.png[問題/行為] 貓晚上進房間會不會有憋尿問題
icon.pngRe: [閒聊] 選了錯誤的女孩成為魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一張
icon.png[心得] EMS高領長版毛衣.墨小樓MC1002
icon.png[分享] 丹龍隔熱紙GE55+33+22
icon.png[問題] 清洗洗衣機
icon.png[尋物] 窗台下的空間
icon.png[閒聊] 双極の女神1 木魔爵
icon.png[售車] 新竹 1997 march 1297cc 白色 四門
icon.png[討論] 能從照片感受到攝影者心情嗎
icon.png[狂賀] 賀賀賀賀 賀!島村卯月!總選舉NO.1
icon.png[難過] 羨慕白皮膚的女生
icon.png閱讀文章
icon.png[黑特]
icon.png[問題] SBK S1安裝於安全帽位置
icon.png[分享] 舊woo100絕版開箱!!
icon.pngRe: [無言] 關於小包衛生紙
icon.png[開箱] E5-2683V3 RX480Strix 快睿C1 簡單測試
icon.png[心得] 蒼の海賊龍 地獄 執行者16PT
icon.png[售車] 1999年Virage iO 1.8EXi
icon.png[心得] 挑戰33 LV10 獅子座pt solo
icon.png[閒聊] 手把手教你不被桶之新手主購教學
icon.png[分享] Civic Type R 量產版官方照無預警流出
icon.png[售車] Golf 4 2.0 銀色 自排
icon.png[出售] Graco提籃汽座(有底座)2000元誠可議
icon.png[問題] 請問補牙材質掉了還能再補嗎?(台中半年內
icon.png[問題] 44th 單曲 生寫竟然都給重複的啊啊!
icon.png[心得] 華南紅卡/icash 核卡
icon.png[問題] 拔牙矯正這樣正常嗎
icon.png[贈送] 老莫高業 初業 102年版
icon.png[情報] 三大行動支付 本季掀戰火
icon.png[寶寶] 博客來Amos水蠟筆5/1特價五折
icon.pngRe: [心得] 新鮮人一些面試分享
icon.png[心得] 蒼の海賊龍 地獄 麒麟25PT
icon.pngRe: [閒聊] (君の名は。雷慎入) 君名二創漫畫翻譯
icon.pngRe: [閒聊] OGN中場影片:失蹤人口局 (英文字幕)
icon.png[問題] 台灣大哥大4G訊號差
icon.png[出售] [全國]全新千尋侘草LED燈, 水草

請輸入看板名稱,例如:Tech_Job站內搜尋

TOP