careening = True
foehn = 0
atk_bonus = .1*foehn
hp_bonus = .2*careening + .1*foehn
trials = 10000
army = {
'Galley': {'n': 10, 'atk': 150, 'hp': 2000, 'fire': 0, 'chains': 0},
'Galleon': {'n': 1, 'atk': 500, 'hp': 15000, 'fire': 0, 'chains': 0},
'Fireship': {'n': 1, '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': 11, 'atk': 200, 'hp': 3000, 'fire': 0, 'chains': 0},
'Turtle Ship': {'n': 3, '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]
army_totals = {}
for a in A: army_totals[a] = 0
for unit in army:
for a in A: army_totals[a] += army[unit]['n'] * army[unit][a]
from numpy.random import rand
def randomize(n): return n * (1 + .25 * (rand() - rand()))
def fight(d):
my_hp = army_totals['hp'] * (1+hp_bonus)
enemy_hp = enemy_totals['hp']
fire, chains = 0, 0
i = 0
while enemy_hp > 0 and i < 50:
i += 1
chains += randomize(army_totals['chains'])
fire += randomize(army_totals['fire'])
atk = randomize(army_totals['atk'])
my_dmg = (fire + atk)*(1+atk_bonus)
enemy_dmg = max(randomize(enemy_totals['atk']) - chains, 0)
enemy_hp -= my_dmg
my_hp -= enemy_dmg
return my_hp > 0
print sum(map(fight, range(trials))), "/", trials
Y2FyZWVuaW5nID0gVHJ1ZQpmb2VobiA9IDAKCmF0a19ib251cyA9IC4xKmZvZWhuCmhwX2JvbnVzID0gLjIqY2FyZWVuaW5nICsgLjEqZm9laG4KCnRyaWFscyA9IDEwMDAwCgphcm15ID0gewogICdHYWxsZXknOiAgIHsnbic6ICAgMTAsICdhdGsnOiAxNTAsICdocCc6ICAyMDAwLCAnZmlyZSc6ICAgMCwgJ2NoYWlucyc6ICAgMH0sCiAgJ0dhbGxlb24nOiAgeyduJzogICAxLCAnYXRrJzogNTAwLCAnaHAnOiAxNTAwMCwgJ2ZpcmUnOiAgIDAsICdjaGFpbnMnOiAgIDB9LAogICdGaXJlc2hpcCc6IHsnbic6ICAgMSwgJ2F0ayc6ICAgMCwgJ2hwJzogIDQwMDAsICdmaXJlJzogMTAwLCAnY2hhaW5zJzogICAwfSwKICAnQ2FyYXZlbCc6ICB7J24nOiAgIDAsICdhdGsnOiAyMDAsICdocCc6ICA2MDAwLCAnZmlyZSc6ICAgMCwgJ2NoYWlucyc6ICA3NX0KfQplbmVtaWVzID0gewogICdCb2F0JzogICAgICAgIHsnbic6IDAsICdhdGsnOiAgIDUwLCAnaHAnOiAgIDMwMCwgJ2ZpcmUnOiAwLCAnY2hhaW5zJzogMH0sCiAgJ1RyaXJlbWUnOiAgICAgeyduJzogMTEsICdhdGsnOiAgMjAwLCAnaHAnOiAgMzAwMCwgJ2ZpcmUnOiAwLCAnY2hhaW5zJzogMH0sCiAgJ1R1cnRsZSBTaGlwJzogeyduJzogMywgJ2F0ayc6ICAxNTAsICdocCc6IDEwMDAwLCAnZmlyZSc6IDAsICdjaGFpbnMnOiAwfSwKICAnQmxhc3QgU2hpcCc6ICB7J24nOiAwLCAnYXRrJzogIDgwMCwgJ2hwJzogIDgwMDAsICdmaXJlJzogMCwgJ2NoYWlucyc6IDB9Cn0KCkEgPSBbJ2F0aycsICdocCcsICdmaXJlJywgJ2NoYWlucyddCmVuZW15X3RvdGFscyA9IHt9CmZvciBhIGluIEE6IGVuZW15X3RvdGFsc1thXSA9IDAKZm9yIHVuaXQgaW4gZW5lbWllczoKICBmb3IgYSBpbiBBOiBlbmVteV90b3RhbHNbYV0gKz0gZW5lbWllc1t1bml0XVsnbiddICogZW5lbWllc1t1bml0XVthXQoKYXJteV90b3RhbHMgPSB7fQpmb3IgYSBpbiBBOiBhcm15X3RvdGFsc1thXSA9IDAKZm9yIHVuaXQgaW4gYXJteToKICBmb3IgYSBpbiBBOiBhcm15X3RvdGFsc1thXSArPSBhcm15W3VuaXRdWyduJ10gKiBhcm15W3VuaXRdW2FdCgpmcm9tIG51bXB5LnJhbmRvbSBpbXBvcnQgcmFuZApkZWYgcmFuZG9taXplKG4pOiByZXR1cm4gbiAqICgxICsgLjI1ICogKHJhbmQoKSAtIHJhbmQoKSkpCmRlZiBmaWdodChkKToKICBteV9ocCA9IGFybXlfdG90YWxzWydocCddICogKDEraHBfYm9udXMpCiAgZW5lbXlfaHAgPSBlbmVteV90b3RhbHNbJ2hwJ10KICBmaXJlLCBjaGFpbnMgPSAwLCAwCiAgaSA9IDAKICB3aGlsZSBlbmVteV9ocCA+IDAgYW5kIGkgPCA1MDoKICAgIGkgKz0gMQogICAgY2hhaW5zICs9IHJhbmRvbWl6ZShhcm15X3RvdGFsc1snY2hhaW5zJ10pCiAgICBmaXJlICs9IHJhbmRvbWl6ZShhcm15X3RvdGFsc1snZmlyZSddKQogICAgYXRrID0gcmFuZG9taXplKGFybXlfdG90YWxzWydhdGsnXSkKICAgIG15X2RtZyA9IChmaXJlICsgYXRrKSooMSthdGtfYm9udXMpCiAgICBlbmVteV9kbWcgPSBtYXgocmFuZG9taXplKGVuZW15X3RvdGFsc1snYXRrJ10pIC0gY2hhaW5zLCAwKQogICAgZW5lbXlfaHAgLT0gbXlfZG1nCiAgICBteV9ocCAtPSBlbmVteV9kbWcKICByZXR1cm4gbXlfaHAgPiAwCnByaW50IHN1bShtYXAoZmlnaHQsIHJhbmdlKHRyaWFscykpKSwgIi8iLCB0cmlhbHMg