魷魚遊戲玻璃橋模擬

魷魚遊戲玻璃橋

18關,16人

各玩家遵守規則,依次前行。

模擬各人存活概率。

import numpy as np

class gamer:
    def __init__(self):
        # 表示該遊戲者已經走過了x塊玻璃
        self.x = 0
    
    def choose(self):
        # 表示該遊戲者做出選擇,0表示選擇左邊,1表示選擇右邊
        self.choice = np.random.randint(0,2)
    
    def go(self,glass):
        for i in range(0,glass.n):
            if glass.walked[i]==1:
                self.x += 1
            else:
                glass.walked[i]=1
                self.choose()
                if self.choice==glass.which[i]:
                    self.x += 1
                else:
                    break
        self.final = self.x #該玩家最終通過了的玻璃數量
        return glass

class glass:
    def __init__(self,n):
        self.n = n
        # 生成n個0或1的隨機數,0表示左邊爲安全的玻璃,1則表示右邊
        self.which = np.random.randint(0,2,size=n)
        # 該序列表示每一塊玻璃是否被走過,初始置0表示未被走過
        self.walked = np.zeros(n)

n=18
k=16

#%%
p = np.zeros(k)
sim = 1000000
for s in range(0,sim):
    if s%10000==0:
        print(s)
    g = glass(n)
    gamers = [gamer() for _ in range(0,k)]
    for i in range(0,k):
        g = gamers[i].go(g)
    for i in range(0,k):
        if gamers[i].final==n:
            p[i] += 1

import matplotlib.pyplot as plt
plt.figure(figsize=(10,5))
plt.bar(range(1,k+1),p/sim)
for a,b in zip(range(1,k+1),p/sim):
    plt.text(a, b,'%.2f%%'%(b*100), ha = 'center',va = 'bottom',fontsize=8)
plt.show()

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章