Math 板


LINE

※ 引述《ntpuisbest (阿龍)》之銘言: : 如題 : 魷魚遊戲第五關玻璃橋 : 到達終點總共要經過18塊玻璃 : 而每次的經過都是兩片玻璃二選一 : 選對了就是強化玻璃 : 選錯了就是掉下去 : 假設選對選錯的機率都是二分之一 : 然後選手總共20位好了 : 再假設每位選手都有超凡記憶力 : 都有辦法記得自己前面的人經過哪些玻璃 : 而且可以趨吉避凶 : 那麼在不互相殘殺的狀況下 : 20位選手的期望通關人數是多少呢 : 我覺得這個問題很複雜 : 因為加入了人有記憶性這個條件後 : 感覺只有用程式模擬 : 配合大數法則才有可能算出來? : 但是程式感覺rule也不太好寫 從 Linearity of Expectation 我們可以知道, 20 通關人數的期望值 = Σ P(第 i 位選手安穩地站在第18塊上) i=1 讓隨機變數 X_i 代表第i位選手最終所能安穩站的玻璃 我舉幾個例子方便大家理解: X_1 = 0 的意思是第1位選手死在第1塊玻璃上 X_1 = 2 的意思是第1位選手死在第3塊玻璃上 X_3 = 4 的意思是第3位選手死在第5塊玻璃上 X_5 = 18的意思是第5位選手成功地站在第18塊玻璃上 20 所以通關人數期望值就是 Σ P(X_i=18) i=1 怎麼算這些 P(X_i=x) 呢? P(X_1=x) = (1/2)^(x+1) , 0<=x<=17 1/2^18 , x =18 0 ,otherwise 幸運的是 X_1,X_2, ... 有 Markovianity,所以可以用下面這方式算出來: P(X_i=x) = Σ P(X_{i-1}=y)P( X_i=x | X_{i-1}=y) 0<=y<=18 對於 i>=2, 我們有 (i) y <= 17 Pr( X_i=x | X_{i-1} = y) = (1/2)^(x-y) , y+1<=x<=17 (1/2)^(17-y) , x = 18 0 , otherwise (ii) y = 18 Pr( X_i=x | X_{i-1} = y) = 1 , x= 18 0 , otherwise 我用程式算出來是 11 。 實際上,這其實只要算到第18個人就可以了。 因為第19跟第20個人從前面的人的結果絕對可以安全通過。 然後我通過計算發現看起來有下面這個情況 for 1<=i<=18 , P(X_i=18) + P(X_{19-i}=18) = 1 弄得更General一點,如果總共有 n 塊玻璃要通過, for 1<=i<=n, P(X_i=n) + P(X_{n+1-i}=n) = 1 如果上面那個猜想是對的,那應該有個不錯的解讀方式可以省下這些計算 這就交給大家想了 -- 因為有大家的支持,才有角卷綿芽的Sololive https://i.imgur.com/CbO6fr2.jpg 直到台灣時間 2021/11/12 (星期五) 下午 9:59 為止都可以在SPWN觀看喔! SPWN連結:https://spwn.jp/events/21101201-jpwatame --



※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 98.45.135.233 (美國)
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Math/M.1634872344.A.D16.html
1F:推 fragmentwing: 推公布程式算 良心碼農 10/22 12:21
2F:推 ntpuisbest : 拜讀一下,感恩 10/22 13:21
3F:推 llrabel : 讚讚 跟我設的隨機變數一樣 10/22 13:31
4F:→ llrabel : 不過原文下已經有人給出另一個「不錯的解讀方式」 10/22 13:33
5F:→ llrabel : 改成用第 i 塊玻璃有沒有死人設隨機變數 Y_i 10/22 13:36
6F:→ llrabel : Y_i = 1 表示有死, Y_i = 0 表示沒有 秒殺 10/22 13:38
你們有誰方便寫一篇用這方式的分析嗎?XD 然後我 Python的程式碼放在這邊,有興趣的可以看看 M = 18 # 玻璃數量 T = 20 # 參加者人數 # 玻璃的數量和參加者人數也可以改成其他的數字 dp = [[0]*(M+1) for _ in range(T+1)] # dp[i][x] = P(X_i=x) for x in range(M): dp[1][x] = 0.5**(x+1) dp[1][M] = 0.5**M for t in range(2,T+1): #calculating dp[t][.] # dp[t][M] for y in range(0,M): dp[t][M] += 0.5**(M-1-y)*dp[t-1][y] dp[t][M] += dp[t-1][M] for x in range(0,M): for y in range(0,x): dp[t][x] += 0.5**(x-y)*dp[t-1][y] result = 0 for t in range(1,T+1): result+=dp[t][M] print(result)
7F:推 fragmentwing: 你文章死在玻璃上的舉例中間兩個是不是打錯啦 10/22 14:02
是第一個打錯XD ※ 編輯: arrenwu (98.45.135.233 美國), 10/22/2021 14:09:31
8F:→ fragmentwing: 原來如此 看懂了 很像國小的以上超過陷阱題XD 10/22 14:24







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燈, 水草

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

TOP