careening = 1
foehn = 0
bonus = {}
bonus['atk'] = .1*foehn
bonus['fire'] = bonus['atk']
bonus['hp'] = .2*careening + .1*foehn
trials = 1000
navy = {
'Galley': {'n': 1, 'atk': 150, 'hp': 2000, 'fire': 0, 'chains': 0},
'Galleon': {'n': 1, 'atk': 500, 'hp': 15000, 'fire': 0, 'chains': 0},
'Fireship': {'n': 8, 'atk': 0, 'hp': 4000, 'fire': 100, 'chains': 0},
'Caravel': {'n': 0, 'atk': 200, 'hp': 6000, 'fire': 0, 'chains': 75}
}
enemies = {
'Boat': {'n': 0, 'atk': 50, 'hp': 300, 'fire': 0, 'chains': 0},
'Trireme': {'n': 14, 'atk': 200, 'hp': 3000, 'fire': 0, 'chains': 0},
'Turtle Ship': {'n': 4, 'atk': 150, 'hp': 10000, 'fire': 0, 'chains': 0},
'Blast Ship': {'n': 0, 'atk': 800, 'hp': 8000, 'fire': 0, 'chains': 0}
}
A = ['atk', 'hp', 'fire', 'chains']
enemy_totals = {}
for a in A: enemy_totals[a] = 0
for unit in enemies:
for a in A: enemy_totals[a] += enemies[unit]['n'] * enemies[unit][a]
#print enemy_totals
navy_totals = {}
for a in A: navy_totals[a] = 0
for unit in navy:
for a in A: navy_totals[a] += navy[unit]['n'] * navy[unit][a] * (1+bonus[a] if a in bonus else 1)
#print navy_totals
from numpy.random import rand
def randomize(n): return n * (1 + 1.0/3 * (rand() - rand()))
def fight(d):
my_hp = navy_totals['hp']
enemy_hp = enemy_totals['hp']
fire, chains = 0, 0
i = 0
while enemy_hp > 0 and i < 50:
i += 1
fire += randomize(navy_totals['fire'])
atk = randomize(navy_totals['atk'])
my_dmg = (fire + atk)
enemy_dmg = max(randomize(enemy_totals['atk']) - chains, 0)
enemy_hp -= my_dmg
my_hp -= enemy_dmg
chains += randomize(navy_totals['chains'])
return my_hp > 0 and i < 50
print sum(map(fight, range(trials))), "/", trials
Y2FyZWVuaW5nID0gMQpmb2VobiA9IDAKCmJvbnVzID0ge30KYm9udXNbJ2F0ayddID0gLjEqZm9laG4KYm9udXNbJ2ZpcmUnXSA9IGJvbnVzWydhdGsnXQpib251c1snaHAnXSA9IC4yKmNhcmVlbmluZyArIC4xKmZvZWhuCgp0cmlhbHMgPSAxMDAwCgpuYXZ5ID0gewogICdHYWxsZXknOiAgIHsnbic6ICAgMSwgJ2F0ayc6IDE1MCwgJ2hwJzogIDIwMDAsICdmaXJlJzogICAwLCAnY2hhaW5zJzogICAwfSwKICAnR2FsbGVvbic6ICB7J24nOiAgIDEsICdhdGsnOiA1MDAsICdocCc6IDE1MDAwLCAnZmlyZSc6ICAgMCwgJ2NoYWlucyc6ICAgMH0sCiAgJ0ZpcmVzaGlwJzogeyduJzogICA4LCAnYXRrJzogICAwLCAnaHAnOiAgNDAwMCwgJ2ZpcmUnOiAxMDAsICdjaGFpbnMnOiAgIDB9LAogICdDYXJhdmVsJzogIHsnbic6ICAgMCwgJ2F0ayc6IDIwMCwgJ2hwJzogIDYwMDAsICdmaXJlJzogICAwLCAnY2hhaW5zJzogIDc1fQp9CmVuZW1pZXMgPSB7CiAgJ0JvYXQnOiAgICAgICAgeyduJzogIDAsICdhdGsnOiAgIDUwLCAnaHAnOiAgIDMwMCwgJ2ZpcmUnOiAwLCAnY2hhaW5zJzogMH0sCiAgJ1RyaXJlbWUnOiAgICAgeyduJzogMTQsICdhdGsnOiAgMjAwLCAnaHAnOiAgMzAwMCwgJ2ZpcmUnOiAwLCAnY2hhaW5zJzogMH0sCiAgJ1R1cnRsZSBTaGlwJzogeyduJzogNCwgJ2F0ayc6ICAxNTAsICdocCc6IDEwMDAwLCAnZmlyZSc6IDAsICdjaGFpbnMnOiAwfSwKICAnQmxhc3QgU2hpcCc6ICB7J24nOiAgMCwgJ2F0ayc6ICA4MDAsICdocCc6ICA4MDAwLCAnZmlyZSc6IDAsICdjaGFpbnMnOiAwfQp9CgpBID0gWydhdGsnLCAnaHAnLCAnZmlyZScsICdjaGFpbnMnXQplbmVteV90b3RhbHMgPSB7fQpmb3IgYSBpbiBBOiBlbmVteV90b3RhbHNbYV0gPSAwCmZvciB1bml0IGluIGVuZW1pZXM6CiAgZm9yIGEgaW4gQTogZW5lbXlfdG90YWxzW2FdICs9IGVuZW1pZXNbdW5pdF1bJ24nXSAqIGVuZW1pZXNbdW5pdF1bYV0KI3ByaW50IGVuZW15X3RvdGFscwpuYXZ5X3RvdGFscyA9IHt9CmZvciBhIGluIEE6IG5hdnlfdG90YWxzW2FdID0gMApmb3IgdW5pdCBpbiBuYXZ5OgogIGZvciBhIGluIEE6IG5hdnlfdG90YWxzW2FdICs9IG5hdnlbdW5pdF1bJ24nXSAqIG5hdnlbdW5pdF1bYV0gKiAoMStib251c1thXSBpZiBhIGluIGJvbnVzIGVsc2UgMSkKI3ByaW50IG5hdnlfdG90YWxzCmZyb20gbnVtcHkucmFuZG9tIGltcG9ydCByYW5kCmRlZiByYW5kb21pemUobik6IHJldHVybiBuICogKDEgKyAxLjAvMyAqIChyYW5kKCkgLSByYW5kKCkpKQpkZWYgZmlnaHQoZCk6CiAgbXlfaHAgPSBuYXZ5X3RvdGFsc1snaHAnXQogIGVuZW15X2hwID0gZW5lbXlfdG90YWxzWydocCddCiAgZmlyZSwgY2hhaW5zID0gMCwgMAogIGkgPSAwCiAgd2hpbGUgZW5lbXlfaHAgPiAwIGFuZCBpIDwgNTA6CiAgICBpICs9IDEKICAgIGZpcmUgKz0gcmFuZG9taXplKG5hdnlfdG90YWxzWydmaXJlJ10pCiAgICBhdGsgPSByYW5kb21pemUobmF2eV90b3RhbHNbJ2F0ayddKQogICAgbXlfZG1nID0gKGZpcmUgKyBhdGspCiAgICBlbmVteV9kbWcgPSBtYXgocmFuZG9taXplKGVuZW15X3RvdGFsc1snYXRrJ10pIC0gY2hhaW5zLCAwKQogICAgZW5lbXlfaHAgLT0gbXlfZG1nCiAgICBteV9ocCAtPSBlbmVteV9kbWcKICAgIGNoYWlucyArPSByYW5kb21pemUobmF2eV90b3RhbHNbJ2NoYWlucyddKQogIHJldHVybiBteV9ocCA+IDAgYW5kIGkgPCA1MApwcmludCBzdW0obWFwKGZpZ2h0LCByYW5nZSh0cmlhbHMpKSksICIvIiwgdHJpYWxzIA==