careening = True
foehn = 31
atk_bonus = .1*foehn
hp_bonus = .2*careening + .1*foehn
trials = 10000
army = {
'Galley': {'n': 0, 'atk': 150, 'hp': 2000, 'fire': 0, 'chains': 0},
'Galleon': {'n': 10, 'atk': 500, 'hp': 15000, 'fire': 0, 'chains': 0},
'Fireship': {'n': 10, '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': 116, 'atk': 200, 'hp': 3000, 'fire': 0, 'chains': 0},
'Turtle Ship': {'n': 58, 'atk': 150, 'hp': 10000, 'fire': 0, 'chains': 0},
'Blast Ship': {'n': 30, '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
Y2FyZWVuaW5nID0gVHJ1ZQpmb2VobiA9IDMxCgphdGtfYm9udXMgPSAuMSpmb2VobgpocF9ib251cyA9IC4yKmNhcmVlbmluZyArIC4xKmZvZWhuCgp0cmlhbHMgPSAxMDAwMAoKYXJteSA9IHsKICAnR2FsbGV5JzogICB7J24nOiAgIDAsICdhdGsnOiAxNTAsICdocCc6ICAyMDAwLCAnZmlyZSc6ICAgMCwgJ2NoYWlucyc6ICAgMH0sCiAgJ0dhbGxlb24nOiAgeyduJzogICAxMCwgJ2F0ayc6IDUwMCwgJ2hwJzogMTUwMDAsICdmaXJlJzogICAwLCAnY2hhaW5zJzogICAwfSwKICAnRmlyZXNoaXAnOiB7J24nOiAgIDEwLCAnYXRrJzogICAwLCAnaHAnOiAgNDAwMCwgJ2ZpcmUnOiAxMDAsICdjaGFpbnMnOiAgIDB9LAogICdDYXJhdmVsJzogIHsnbic6ICAgMCwgJ2F0ayc6IDIwMCwgJ2hwJzogIDYwMDAsICdmaXJlJzogICAwLCAnY2hhaW5zJzogIDc1fQp9CmVuZW1pZXMgPSB7CiAgJ0JvYXQnOiAgICAgICAgeyduJzogMCwgJ2F0ayc6ICAgNTAsICdocCc6ICAgMzAwLCAnZmlyZSc6IDAsICdjaGFpbnMnOiAwfSwKICAnVHJpcmVtZSc6ICAgICB7J24nOiAxMTYsICdhdGsnOiAgMjAwLCAnaHAnOiAgMzAwMCwgJ2ZpcmUnOiAwLCAnY2hhaW5zJzogMH0sCiAgJ1R1cnRsZSBTaGlwJzogeyduJzogNTgsICdhdGsnOiAgMTUwLCAnaHAnOiAxMDAwMCwgJ2ZpcmUnOiAwLCAnY2hhaW5zJzogMH0sCiAgJ0JsYXN0IFNoaXAnOiAgeyduJzogMzAsICdhdGsnOiAgODAwLCAnaHAnOiAgODAwMCwgJ2ZpcmUnOiAwLCAnY2hhaW5zJzogMH0KfQoKQSA9IFsnYXRrJywgJ2hwJywgJ2ZpcmUnLCAnY2hhaW5zJ10KZW5lbXlfdG90YWxzID0ge30KZm9yIGEgaW4gQTogZW5lbXlfdG90YWxzW2FdID0gMApmb3IgdW5pdCBpbiBlbmVtaWVzOgogIGZvciBhIGluIEE6IGVuZW15X3RvdGFsc1thXSArPSBlbmVtaWVzW3VuaXRdWyduJ10gKiBlbmVtaWVzW3VuaXRdW2FdCgphcm15X3RvdGFscyA9IHt9CmZvciBhIGluIEE6IGFybXlfdG90YWxzW2FdID0gMApmb3IgdW5pdCBpbiBhcm15OgogIGZvciBhIGluIEE6IGFybXlfdG90YWxzW2FdICs9IGFybXlbdW5pdF1bJ24nXSAqIGFybXlbdW5pdF1bYV0KCmZyb20gbnVtcHkucmFuZG9tIGltcG9ydCByYW5kCmRlZiByYW5kb21pemUobik6IHJldHVybiBuICogKDEgKyAuMjUgKiAocmFuZCgpIC0gcmFuZCgpKSkKZGVmIGZpZ2h0KGQpOgogIG15X2hwID0gYXJteV90b3RhbHNbJ2hwJ10gKiAoMStocF9ib251cykKICBlbmVteV9ocCA9IGVuZW15X3RvdGFsc1snaHAnXQogIGZpcmUsIGNoYWlucyA9IDAsIDAKICBpID0gMAogIHdoaWxlIGVuZW15X2hwID4gMCBhbmQgaSA8IDUwOgogICAgaSArPSAxCiAgICBjaGFpbnMgKz0gcmFuZG9taXplKGFybXlfdG90YWxzWydjaGFpbnMnXSkKICAgIGZpcmUgKz0gcmFuZG9taXplKGFybXlfdG90YWxzWydmaXJlJ10pCiAgICBhdGsgPSByYW5kb21pemUoYXJteV90b3RhbHNbJ2F0ayddKQogICAgbXlfZG1nID0gKGZpcmUgKyBhdGspKigxK2F0a19ib251cykKICAgIGVuZW15X2RtZyA9IG1heChyYW5kb21pemUoZW5lbXlfdG90YWxzWydhdGsnXSkgLSBjaGFpbnMsIDApCiAgICBlbmVteV9ocCAtPSBteV9kbWcKICAgIG15X2hwIC09IGVuZW15X2RtZwogIHJldHVybiBteV9ocCA+IDAKcHJpbnQgc3VtKG1hcChmaWdodCwgcmFuZ2UodHJpYWxzKSkpLCAiLyIsIHRyaWFscyA=