抽數期望值計算機 - 手遊

2022/12/10 11:32 更新：已排除了溢位問題，抽數設定無限大應該都可以正常運行

----------------------寫在前面----------------------

1. 本文含有大量數學算式，範圍涵蓋條件機率、二項式定理等數學觀念。
若對數學有強烈過敏或排斥反應，建議直接END看結論或是左轉離開。
2. 本文不含任何統計學觀念，完全只有機率。
3. 本文不考慮任何因抽卡玄學帶來的機率變動影響，也一概不討論玄學對機率之影響。
4. 文末提供之工作表在使用上很簡單，完全不知道這篇在講什麼也可以快樂(?)使用。
5. 本文忽略保底影響。
6. 本文為了計算方便，會將抽卡行為理想化，請自行斟酌下修實際機率。

----------------------------------------------------

大家好，無聊的數學家又來荼毒......分享一些有趣東西給大家了
這件事情本來是一個單純想法
「單抽中PU五星的機率是0.008，那是否1000抽就可期望抽中8隻PU五星？」
「實際上，1000抽中8隻以上PU五星的機率有多少？」

對於一個無課/福袋課玩家而言，石頭是有限的。
一年從營運收到的免費石，可能是幾百~一千多抽的程度。
那麼無課玩家就會想知道：
1. 假如我在一年內想抽到X個五星，用現有的抽數，是否太樂觀了？
2. 我總共有Y抽，在目前的卡池抽到我想要的角色的機率有多少？
3. 一年過去了，我抽到了Z個五星，算是哪一個大洲的人？
以下問題將會一一解答，並提供計算機，大家可以自己算。
(不想看數學的可以END了)

----------------------------------------------------

一切問題從簡單的範例開始。
假如我們用一個紅白球箱子來模擬卡池。
箱子內有4個紅球，1個白球。
抽中紅球算中獎(真容易中獎)，抽完球放回箱中。

那麼抽中紅球機率=4/5=0.8，抽中白球機率=1/5=0.2。

假如只抽1次，抽中紅球機率自然就是0.8。
抽2次的話就開始有趣了，根據結果，會有四種情形：

P(紅紅) = 0.8*0.8 = 0.64
P(紅白) = 0.8*0.2 = 0.16
P(白紅) = 0.2*0.8 = 0.16
P(白白) = 0.2*0.2 = 0.04

當然 0.64 + 0.16 + 0.16 + 0.04 = 1。
抽中兩紅的機率，是P(紅紅) = 0.64。
抽中兩白的機率，是P(白白) = 0.04。
一紅一白的機率，是P(紅白)+P(白紅)=0.16+0.16=0.32。
由於抽中順序不影響其最終機率，所以通常會合併看待。

----------------------------------------------------

讓我們用另一種觀點來看待此事。
下面是一個簡單的算式：

(0.8+0.2)^2

我們當然知道這個答案是1，因為0.8+0.2=1，
且1不管自乘幾次，答案都是1。
但是我們假如直接展開算式，就可以發現與上面的紅白球機率對應：

(0.8+0.2)^2
= 0.8*0.8 + 0.8*0.2 + 0.2*0.8 + 0.2*0.2 = 1。

而(0.8+0.2)^3的機率，也正好與抽樣三次的情況對應：

(0.8+0.2)^3
= 0.8*0.8*0.8 + 0.8*0.8*0.2 + 0.8*0.2*0.8 + 0.8*0.2*0.2
+ 0.2*0.8*0.8 + 0.2*0.8*0.2 + 0.2*0.2*0.8 + 0.2*0.2*0.2
= P(紅紅紅) + P(紅紅白) + P(紅白紅) + P(紅白白)
+ P(白紅紅) + P(白紅白) + P(白白紅) + P(白白白)
= 1。

於此可以類推，抽1000次的機率分布可以由
(0.8+0.2)^1000
展開得出。

可是單純展開仍然無法減少我們的工作，
所以到了本文另一主角：二項式定理出場。

----------------------------------------------------

二項式定理

(x+y)^n 可以展開為 C(n,k)*x^k*y^(n-k) 的級數和。

礙於數學公式不好在批踢踢上打出來，直接引用wiki截圖：

https://i.imgur.com/zmirw96.png

Ref: Wiki二項式定理條目
https://reurl.cc/10keeG

簡單用上面的例子類推
(0.8+0.2)^3
= C(3,0)*(0.8)^3*(0.2)^0
+ C(3,1)*(0.8)^2*(0.2)^1
+ C(3,2)*(0.8)^1*(0.2)^2
+ C(3,3)*(0.8)^0*(0.2)^3

= 1*(0.8)^3
+ 3*(0.8)^2*(0.2)
+ 3*(0.8)*(0.2)^2
+ 1*(0.2)^3

= 1*P(三紅) + 3*P(兩紅) + 3*P(一紅) + 1*P(無紅)

= 1。



那麼我們把中獎機率換成代數，當抽1000次時，就會變成
(不中+中)^1000
= C(1000,0)*P(不中) +
C(1000,1)*P(中一個) +
C(1000,2)*P(中兩個) +
C(1000,3)*P(中三個) +
... +
C(1000,1000)*P(全中)

至此，準備工作完成。

----------------------------------------------------

回到原本問題

1000抽中8隻以上PU五星的機率有多少？

此機率可以看成
P(1000抽中8隻以上PU五星)
= 1 - P(1000抽中0~7隻PU五星)
= 1 - [ C(1000,0)*P(不中) +
C(1000,1)*P(中一個) +
C(1000,2)*P(中兩個) +
... +
C(1000,7)*P(中七個)]

至於 C(1000,N)*P(中N個) 怎麼算？
假如1000抽中1個，表示另外999次都不中
機率就是 C(1000,1)*(抽中)^1*(不中)^999
中二就是 C(1000,2)*(抽中)^2*(不中)^998
中三就是 C(1000,3)*(抽中)^3*(不中)^997
...
以此類推。

完整算式列一下

P(1000抽中8隻以上PU五星)
= 1 - P(1000抽中0~7隻PU五星)
= 1 - [ C(1000,0)*P(不中) +
C(1000,1)*P(中一個) +
C(1000,2)*P(中兩個) +
C(1000,3)*P(中三個) +
C(1000,4)*P(中四個) +
C(1000,5)*P(中五個) +
C(1000,6)*P(中六個) +
C(1000,7)*P(中七個)]

= 1 - [ C(1000,0)*(抽中)^0*(不中)^1000 +
C(1000,1)*(抽中)^1*(不中)^999 +
C(1000,2)*(抽中)^2*(不中)^998 +
C(1000,3)*(抽中)^3*(不中)^997 +
C(1000,4)*(抽中)^4*(不中)^996 +
C(1000,5)*(抽中)^5*(不中)^995 +
C(1000,6)*(抽中)^6*(不中)^994 +
C(1000,7)*(抽中)^7*(不中)^993]

再把實際機率: 抽中=0.008，不中=0.992 代入算式中
大功告成了，算式裡面已經沒有任何未知數了。
於是我們可以把以上過程編進Google工作表中，
然後填上中獎機率與抽數，得出想要的結果。

---------------END請回到這裡看結論---------------

於是乎，經由工作表的努力，得出的結果如下圖

https://i.imgur.com/VnQPWlM.png

查表可以得知 P(1000抽中8隻以上PU五星) = 54.76%。
也就是說如果你的運氣是排名前54.76%的玩家，
你在任何卡池加總丟1000抽，可以得到8隻以上的PU五星。
(一半以上的玩家可以期待，1000抽中8隻以上PU五星呢)
如果連抽歪五星也想計入呢？把中獎率從0.008調成0.01即可。

另外這個表格也可以計算單個卡池的期望值
假如你想抽奧伯龍池，準備了200抽，
想知道這樣抽中寶1~5的機率有多少
把機率設為0.008，抽數設為200即可。

工作表如下，請自行另存副本/下載取用。

https://reurl.cc/7jbRk9

另外我的公式有寫活，如果覺得中10隻太少想計算更多隻，
選取最後一行往右拖拉即可。

----------------------寫在最後----------------------

1. 此表格有理想化的地方，
在於一年內玩家能領的石頭不是一開始就能領到，而是依時間增加的。
所以途中抽到一半沒石頭的情況，本表無法顧及。
但是判斷單個卡池的期望值不受此影響，可以放心使用。

2. 此表格可套用於其他遊戲卡池，如果該卡池的機率公正的話。


--