si999 - 経済統計練習帳

import numpy as np
 
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("松本代官")
number=21
life=[0]
age=[7]
px=[0]
for i in range(number):
    life.append(100)
    px.append(i+1)
for i in range(6):
    for j in range(i+1):
        age.append(6-i)
 
q=[]
for a1 in [1,2,3,4.5,6]:
    for i in range(number):
        h=0
        if life[i+1]>50:
           h=h+1
        if age[i+1]==a1: 
           h=h+1
        if h==2:
           q.append(i+1)
    q1=np.random.choice(q)
    life[q1]=0
 
de=[]
for i in range(number):
    if life[i+1]>50:
       de.append(i+1)
 
 
for j in range(15):
    d2=np.random.choice(de)
    de.remove(d2)
    px[d2]=j+1
 
 
for i in range(6):
    life.append(100)
    px.append(16+i)
    age.append(0)
 
number=number+6
print(px)