むるおか君のPythonパチンコシミュレーション

Pythonを使ってパチンコの新台解析やってます。たまに自分が処理したい作業のGUIも作ったり

P緋弾のアリア4の天井期待値&シミュレーション

# -*- coding: utf-8 -*-
"""
Created on Mon Jan 25 16:40:50 2021

@author: HOST1191
"""
import random
import statistics 
import collections
import matplotlib.pyplot as plt
import seaborn as sns
import csv
#デンチュー振り分け
def d_furi():
    d_furir = random.randint(1,400)
    d_dedama =  0
    d_kakuhen =  0
    if  d_furir  >=  1  and  d_furir <= 199:
        d_dedama +=  360
        d_kakuhen   =  1
    else:
        d_dedama  += 1200
        d_kakuhen  =  1
 
    return  d_dedama,d_kakuhen

#ヘソの振り分け
def h_furi():
    h_furir = random.randint(1,100)
    h_dedama =  0
    h_kakuhen =  0
    if  h_furir  >=  1  and  h_furir  <= 50:
        h_dedama +=  360
        h_kakuhen   =  1
    else:
        h_dedama  +=  360
        h_kakuhen   =  0
    return  h_dedama,h_kakuhen  

#ライトニグバレットチャンス
def l_furi():    
    l_kakuhen =  0
    l_tyusen = random.random()
    if l_tyusen <= (1/4):
        l_kakuhen  =  2
    return  l_kakuhen

#バレットアタックストックタイム
def v_renchan():
    v_st  =  0
    stock =  0
    while  v_st  <  38:
        v_tyusen  =  random.random()
        v_st  += 1
        if  v_tyusen  <=  (1/73):
            stock  += 1
        if  stock  >  4:
            stock  =  4
    return  stock
        
        
#連荘
def renchan(r_kakuhen):
    wk_st    = 0
    wk_jitan  = 0
    r_dedama  =  0
    wk_dedama =  0
    wk_ren  =  0
    end  =  0
#パラメータ説明
#r_kakuhen =  0  確変ではない
#          =  1  確変
#          =  2  ライトニングバレットチャンス 
    while  end  < 1:
        if  r_kakuhen ==  1:
            r_tyusen = random.random()
            wk_st  +=  1
            if r_tyusen <= (1/73):
#               ライトニングバレットチャンス発生時
                l_kakuhen = l_furi()
                if l_kakuhen == 2:
                    v_stock  =  1
                    while  v_stock  > 0:
                        v_stock -= 1
                        wk_ren += 1
                        r_dedama,r_kakuhen = d_furi()
                        wk_dedama  +=  r_dedama
                        v_stock +=  v_renchan()
                else:
                    wk_ren  +=  1
                    r_dedama,r_kakuhen = d_furi()
                    wk_dedama  +=  r_dedama
                wk_st  =  0
            if  wk_st  >=  100:
                end  +=  1
#       確変じゃないとき
        else:
            r_tyusen  =  random.random()
            wk_jitan  += 1
            if r_tyusen <= (1/199.1):
                r_dedama,r_kakuhen  =  h_furi()
                wk_ren  +=  1
                wk_dedama  +=  r_dedama
                wk_jitan  =  0
            if  wk_jitan  >=  54:
                end  +=  1
    return  wk_ren,wk_dedama

sikou  =  100000
#連荘回数
ren1 = []
#ラッシュの出玉数
dedama1 = []
for  sim  in  range(1,sikou+1):
#初当たりの振り分け処理天井考慮  
    sim3  =  0
    end1  =  0
    while  end1  <  1:
        t_tyusen  =  random.random()
        sim3  +=  1
        if  sim3  >=  500:
            d_dedama,d_kakuhen  =  d_furi()
            ren,dedama  =  renchan(d_kakuhen)
            ren  +=  1
            dedama  +=  d_dedama
            end1  +=  1
        else:
            if  t_tyusen  <=  (1/199.1):
                h_dedama,h_kakuhen  =  h_furi()
                ren,dedama  =  renchan(h_kakuhen)
                ren  +=  1
                dedama  +=  h_dedama
                end1  +=  1
    dedama1.append(dedama)
    ren1.append(ren)
#連荘平均数
m1 = statistics.mean(ren1)
m2 = statistics.mean(dedama1)
#連荘と出玉の整理
ren1.sort() 
dedama1.sort() 
#print ("連荘数=",ren1)
#print ("出玉数=",dedama1)
print('天井考慮連荘平均: {0:.2f}'.format(m1))
print('天井考慮出玉平均: {0:.1f}'.format(m2))

ren_bunpu = collections.Counter(ren1)
dedama_bunpu = collections.Counter(dedama1)
print("分布=",ren_bunpu)
#print("分布=",dedama_bunpu)
#print(ren_bunpu)
#出玉の円グラフ作表用定義
x1 = sum( x1 <=  1000 for x1 in dedama1)
x2 = sum(  1000 < x2 <=  5000 for x2 in dedama1) 
x3 = sum(  5000 < x3 <=  10000 for x3 in dedama1)
x4 = sum(  10000 < x4  for x4 in dedama1) 
Xgraph = ["{0:.1f}".format(x1/sikou*100),"{0:.1f}".format(x2/sikou*100),
          "{0:.1f}".format(x3/sikou*100),"{0:.1f}".format(x4/sikou*100)]
label = ["U1k","1k~5k","5k~10k","10kover"]
colors2 = ["0.3", "0.5", "0.7", "0.9", "1.1"]
plt.pie(Xgraph, labels=label, counterclock=False, startangle=90, colors=colors2,autopct="%1.1f%%")
plt.axis('equal')
print("x=",x1,x2,x3,x4)

d = list(ren_bunpu.keys())
e = list(ren_bunpu.values())
sns.set()
sns.set_style("darkgrid", {'grid.lines1tyle': '--'})
sns.set_context("paper", 1.0, {"lines.linewidth": 1})
sns.set_palette("winter_r", 2, 1)
sns.set('talk', 'darkgrid', 'dark', font_scale=1.0,
        rc={"lines.linewidth": 10, 'grid.linestyle': '--'})
fig, ax = plt.subplots(figsize = (10,5))
ax.set(xlabel ='ren', ylabel='occcurs', xlim=(0,30), ylim=(0,50000))
ax.bar(d,e)
plt.show()

#実戦シミュレーション
sikou  =  3000

for  sim  in  range(1,8):
#初当たりの振り分け処理天井考慮
    sim4  =  0
    sim3  =  0
#回転数
    kaiten1 = []
#出玉推移
    dedama1 = []
    soudedama  =  0
    while  sim4  <  sikou:
        t_tyusen  =  random.random()
        sim3  +=  1
        sim4  +=  1
        if  sim3  >=  500:
            d_dedama,d_kakuhen  =  d_furi()
            ren,dedama  =  renchan(d_kakuhen)
            dedama  +=  d_dedama
            soudedama  +=  dedama
            sim3  =  0
        else:
            if  t_tyusen  <=  (1/199.1):
                h_dedama,h_kakuhen  =  h_furi()
                ren,dedama  =  renchan(h_kakuhen)
                dedama  +=  h_dedama
                soudedama  +=  dedama
                sim3  =  0
            else:
                soudedama  =  soudedama  -   (250 / 18)
        dedama1.append(soudedama)
        kaiten1.append(sim4)
    plt.plot(kaiten1,dedama1,linewidth=1)
    plt.show()

#天井期待値の計算
#試行回数
sikou  =  10000
#千円当たりの回転数
kaitenritu =  0
#作図用
sk = []
dk = []
sim1   =  0
sim2   =  0
atari = 0
suru  = 0
for  tenzyo  in  range(0,500,50):
    for kaitenritu  in  range  (15,21):
        soudedama =  0
        for  a  in  range(1,sikou+1):
            suru  =  0
            sim1  =  0
            sim2  =  0
            while  suru  <  1: 
                t_tyusen = random.random()
#               天井到達時点
                if  sim1  ==  tenzyo:
                    t_dedama,t_kakuhen  =  d_furi()
                    ren,dedama  =  renchan(t_kakuhen)
                    ren  +=  1
                    soudedama  +=  t_dedama
                    soudedama  +=  dedama
                    suru  +=  1
                else:
                    sim1  +=  1
                    soudedama  =  soudedama  -  (250 / kaitenritu)
                    if  t_tyusen  <=  1/(199.1):
                        t_dedama,t_kakuhen  =  h_furi()
                        ren,dedama  =  renchan(t_kakuhen)
                        ren  +=  1
                        soudedama  +=  t_dedama
                        soudedama  +=  dedama
                        suru  +=  1
        sk.append(["{0:.0f}".format((soudedama/sikou)*4),tenzyo,kaitenritu])



print  ("作表用=",sk)
with open('data.csv', 'w') as file:
    writer = csv.writer(file, lineterminator='\n')
    writer.writerows(sk)
#plt.plot(atarikai,dedama1)
#plt.show()