ｑ２２２ｘｘ - 経済統計練習帳

import numpy as np
 
def death(anx):
 
    z=[]
    for i in range(number):
        h=0
        if life[i]>60 :
            h=h+1
        if a[i]==anx :
           h=h+1
        if h==2 : 
           z.append(i)
 
    k=np.random.choice(z)
    life[k]=0
 
p = ["大御所", "将軍", "大老","美濃藩主","伊勢藩主","信濃藩主"]
p.append("老中")
p.append("駿河老中")
p.append("美濃守護")
p.append("伊勢守護")
p.append("信濃守護")
p.append("勘定奉行")
p.append("町奉行")
p.append("寺社奉行")
p.append("外国奉行")
p.append("大垣代官")
p.append("岐阜代官")
p.append("松坂代官")
p.append("熊野代官")
p.append("長野代官")
p.append("松本代官")
c=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]
a=[6,5,5,4,4,4,3,3,3,3,2,2,2,2,2,1,1,1,1,1,1]
life=[100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100]
number=21
 
q=np.zeros((15,21))
 
q[0][0]=100
q[0][1]=50
q[0][2]=50
 
q[1][1]=100
q[1][3]=50
q[1][4]=50
q[1][5]=50
 
q[2][2]=100
q[2][3]=50
q[2][4]=50
q[2][5]=50
 
q[3][3]=100
q[3][6]=50
q[3][7]=50
q[3][8]=25
 
q[4][4]=100
q[4][6]=50
q[4][7]=50
q[4][9]=25
 
q[5][5]=100
q[5][6]=50
q[5][7]=50
q[5][10]=25
 
q[6][6]=100
q[6][8]=50
q[6][9]=50
q[6][10]=50
 
q[7][7]=100
q[7][8]=50
q[7][9]=50
q[7][10]=50
 
q[8][8]=100
q[8][11]=50
q[8][12]=50
q[8][13]=50
q[8][14]=50
q[8][15]=25
q[8][16]=25
 
q[9][9]=100
q[9][11]=50
q[9][12]=50
q[9][13]=50
q[9][14]=50
q[9][17]=25
q[9][18]=25
 
q[10][10]=100
q[10][11]=50
q[10][12]=50
q[10][13]=50
q[10][14]=50
q[10][19]=25
q[10][20]=25
 
 
sig=[]
for tr in range(10):
 
    for a1 in range(1,7):
        death(a1)
 
 
    pol=[]
    for i in range(number):
        if life[i]>50:pol.append(i)
 
    for tr in range(15):
        maxp=-999
        sxs=0
        for p1 in pol:
            c1=c[p1]
            q1=q[tr][c1]
            if q1>maxp:sxs=p1
            if q1>maxp:maxp=q1
        c[sxs]=tr
        pol.remove(sxs)
 
    for i in range(6):
        a.append(0)
        life.append(100)
        c.append(i+15)
 
    number=number+6
    for i in range(number):
        a[i]=a[i]+1
 
    memo=[]
    for i in range(21):
        for j in range(number):
           h=0
           if life[j]>60:h=h+1
           if c[j]==i:h=h+1    
           if h==2:memo.append(j)
    sig.append(memo)
 
print(sig)
 
for sp in sig:
    sx=0
    for c1 in sp:
        if c1==21:print(p[sx])
        sx=sx+1