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

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

Pythonによる仕置人天井シミュレーション

# -*- coding: utf-8 -*-
"""
Created on Wed Jul  1 10:56:29 2020

@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,200)
    d_dedama =  0
    d_kakuhen =  0
    d_jitan  =  0
    if  d_furir  >=  1  and  d_furir  <= 33:
        d_dedama +=  900
        d_kakuhen   =  1
        d_jitan    =  1
    elif d_furir  >=  34  and  d_furir  <=  99:
        d_dedama  += 900
        d_kakuhen  =  1
        d_jitan   =  0
    elif  d_furir  >=  100  and  d_furir  <= 132:
        d_dedama  +=  270
        d_kakuhen  =  1
        d_jitan  =  1
    elif  d_furir  >=  133  and  d_furir  <= 198:
        d_dedama  =  270
        d_kakuhen =  1
        d_jitan  =   0
    elif  d_furir  ==  199:
        d_dedama  +=  900
        d_kakuhen  =  0
        d_jitan  =  2
    else:
        d_dedama  += 270
        d_kakuhen  =  0
        d_jitan  =  2
    return  d_dedama,d_kakuhen,d_jitan

#仕置人ヘソの振り分け
def h_furi():
    h_furir = random.randint(1,100)
    h_dedama =  0
    h_kakuhen =  0
    h_jitan  =  0
    if  h_furir  >=  1  and  h_furir  <= 99:
        h_dedama +=  270
        h_kakuhen   =  0
        h_jitan    =  3
    else:
        h_dedama  += 270
        h_kakuhen  =  1
        h_jitan  =  0
    return  h_dedama,h_kakuhen,h_jitan   
#連荘
def renchan(r_kakuhen,r_jitan):
    wk_jitan = 0
    wk_st    = 0
    r_dedama  =  0
    wk_dedama =  0
    wk_ren  =  0
    end  =  0
#パラメータ説明
#r_jitan  =  0  時短無し
#         =  1 時短50回
#         =  2  時短279回
#r_kakuhen =  0  確変ではない
#          =  1  確変
    while  end  < 1:
        if  r_kakuhen  ==  1:
            r_tyusen = random.random()
            wk_st  +=  1
#            dedama  =  dedama  -  0.5
            if r_tyusen <= (1/86.3):
                wk_ren  +=  1
                r_dedama,r_kakuhen,r_jitan = d_furi() 
                wk_dedama  +=  r_dedama
                wk_st  =  0
                wk_jitan  =  0
            else:
                if  wk_st  >  120:
                    if  r_jitan   !=  0:
                        r_kakuhen  =  0
                        wk_jitan   =  0
                    else:
                        end  +=  1
        #時短の処理
        else:
            #引き戻し抽選
            r_tyusen2 = random.random()
            wk_jitan  +=  1
            if  r_tyusen2  <=  (1/99.9):
                wk_ren  +=  1
                r_dedama,r_kakuhen,r_jitan = d_furi()
                wk_dedama  +=  r_dedama
                wk_st  =  0
                wk_jitan  =  0
                                
            else:
                if  r_jitan  ==  1:
                    if  wk_jitan  >=  54:
                        end  +=  1
                elif r_jitan  == 2:
                    if  wk_jitan  >=  663:
                        end  +=  1
                elif r_jitan  == 3:
                    if  wk_jitan  >=  44:
                        end  +=  1
                elif r_jitan  == 4:
                    if  wk_jitan  >= 383:
                        end  +=  1
                else:
                    print ("時短フラグ異常=",r_jitan)
                    end  +=  1
    return  wk_ren,wk_dedama

sikou  =  100000
#連荘回数
ren1 = []
#ラッシュの出玉数
dedama1 = []
no1  =  0
for  sim  in  range(1,sikou+1):
#初当たりの振り分け処理  
    t_tyusen = random.random()
    if  t_tyusen  <=  (94/100):
        h_dedama,h_kakuhen,h_jitan  =  h_furi()
        ren,dedama  =  renchan(h_kakuhen,h_jitan)
        ren  +=  1
        dedama  +=  h_dedama
    else:
        ren,dedama  = renchan(0,4)
        if  ren  == 0:
            no1  +=  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(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 = ["1000発以下","1000発~5000発","5000発~万発","万発over"]
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')


d = list(ren_bunpu.keys())
e = list(ren_bunpu.values())
sns.set()
sns.set_style("darkgrid", {'grid.linestyle': '--'})
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(xlim=(0,30), ylim=(0,50000))
ax.bar(d,e)
plt.show()

#天井期待値の計算
#試行回数
sikou  =  100000
#千円当たりの回転数
kaitenritu =  0
#作図用
sk = []
dk = []
sim1   =  0
sim2   =  0
atari = 0
suru  = 0
for  tenzyo  in  range(0,281,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:
                    soudedama = soudedama  -  0.5
                    if  t_tyusen  <=  1/(99.9):
                        t_dedama,t_kakuhen,t_jitan  =  d_furi()
                        ren,dedama  =  renchan(t_kakuhen,t_jitan)
                        ren  +=  1
                        soudedama  +=  t_dedama
                        soudedama  +=  dedama
                        suru  +=  1
                    else:
                        sim2  +=  1
                        if  sim2  >=  379:
                            soudedama  +=  0
                            suru +=  1
                else:
                    sim1  +=  1
                    soudedama  =  soudedama  -  (250 / kaitenritu)
                    if  t_tyusen  <=  1/(99.9):
                        t_dedama,t_kakuhen,t_jitan  =  h_furi()
                        ren,dedama  =  renchan(t_kakuhen,t_jitan)
                        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()