careening = 1
foehn = 18
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': 9, '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': 26, 'atk': 200, 'hp': 3000, 'fire': 0, 'chains': 0},
'Turtle Ship': {'n': 13, 'atk': 150, 'hp': 10000, 'fire': 0, 'chains': 0},
'Blast Ship': {'n': 7, '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
Y2FyZWVuaW5nID0gMQpmb2VobiA9IDE4Cgpib251cyA9IHt9CmJvbnVzWydhdGsnXSA9IC4xKmZvZWhuCmJvbnVzWydmaXJlJ10gPSBib251c1snYXRrJ10KYm9udXNbJ2hwJ10gPSAuMipjYXJlZW5pbmcgKyAuMSpmb2VobgoKdHJpYWxzID0gMTAwMAoKbmF2eSA9IHsKICAnR2FsbGV5JzogICB7J24nOiAgIDEsICdhdGsnOiAxNTAsICdocCc6ICAyMDAwLCAnZmlyZSc6ICAgMCwgJ2NoYWlucyc6ICAgMH0sCiAgJ0dhbGxlb24nOiAgeyduJzogICAxLCAnYXRrJzogNTAwLCAnaHAnOiAxNTAwMCwgJ2ZpcmUnOiAgIDAsICdjaGFpbnMnOiAgIDB9LAogICdGaXJlc2hpcCc6IHsnbic6ICAgOSwgJ2F0ayc6ICAgMCwgJ2hwJzogIDQwMDAsICdmaXJlJzogMTAwLCAnY2hhaW5zJzogICAwfSwKICAnQ2FyYXZlbCc6ICB7J24nOiAgIDAsICdhdGsnOiAyMDAsICdocCc6ICA2MDAwLCAnZmlyZSc6ICAgMCwgJ2NoYWlucyc6ICA3NX0KfQplbmVtaWVzID0gewogICdCb2F0JzogICAgICAgIHsnbic6ICAwLCAnYXRrJzogICA1MCwgJ2hwJzogICAzMDAsICdmaXJlJzogMCwgJ2NoYWlucyc6IDB9LAogICdUcmlyZW1lJzogICAgIHsnbic6IDI2LCAnYXRrJzogIDIwMCwgJ2hwJzogIDMwMDAsICdmaXJlJzogMCwgJ2NoYWlucyc6IDB9LAogICdUdXJ0bGUgU2hpcCc6IHsnbic6IDEzLCAnYXRrJzogIDE1MCwgJ2hwJzogMTAwMDAsICdmaXJlJzogMCwgJ2NoYWlucyc6IDB9LAogICdCbGFzdCBTaGlwJzogIHsnbic6ICA3LCAnYXRrJzogIDgwMCwgJ2hwJzogIDgwMDAsICdmaXJlJzogMCwgJ2NoYWlucyc6IDB9Cn0KCkEgPSBbJ2F0aycsICdocCcsICdmaXJlJywgJ2NoYWlucyddCmVuZW15X3RvdGFscyA9IHt9CmZvciBhIGluIEE6IGVuZW15X3RvdGFsc1thXSA9IDAKZm9yIHVuaXQgaW4gZW5lbWllczoKICBmb3IgYSBpbiBBOiBlbmVteV90b3RhbHNbYV0gKz0gZW5lbWllc1t1bml0XVsnbiddICogZW5lbWllc1t1bml0XVthXQojcHJpbnQgZW5lbXlfdG90YWxzCm5hdnlfdG90YWxzID0ge30KZm9yIGEgaW4gQTogbmF2eV90b3RhbHNbYV0gPSAwCmZvciB1bml0IGluIG5hdnk6CiAgZm9yIGEgaW4gQTogbmF2eV90b3RhbHNbYV0gKz0gbmF2eVt1bml0XVsnbiddICogbmF2eVt1bml0XVthXSAqICgxK2JvbnVzW2FdIGlmIGEgaW4gYm9udXMgZWxzZSAxKQojcHJpbnQgbmF2eV90b3RhbHMKZnJvbSBudW1weS5yYW5kb20gaW1wb3J0IHJhbmQKZGVmIHJhbmRvbWl6ZShuKTogcmV0dXJuIG4gKiAoMSArIDEuMC8zICogKHJhbmQoKSAtIHJhbmQoKSkpCmRlZiBmaWdodChkKToKICBteV9ocCA9IG5hdnlfdG90YWxzWydocCddCiAgZW5lbXlfaHAgPSBlbmVteV90b3RhbHNbJ2hwJ10KICBmaXJlLCBjaGFpbnMgPSAwLCAwCiAgaSA9IDAKICB3aGlsZSBlbmVteV9ocCA+IDAgYW5kIGkgPCA1MDoKICAgIGkgKz0gMQogICAgZmlyZSArPSByYW5kb21pemUobmF2eV90b3RhbHNbJ2ZpcmUnXSkKICAgIGF0ayA9IHJhbmRvbWl6ZShuYXZ5X3RvdGFsc1snYXRrJ10pCiAgICBteV9kbWcgPSAoZmlyZSArIGF0aykKICAgIGVuZW15X2RtZyA9IG1heChyYW5kb21pemUoZW5lbXlfdG90YWxzWydhdGsnXSkgLSBjaGFpbnMsIDApCiAgICBlbmVteV9ocCAtPSBteV9kbWcKICAgIG15X2hwIC09IGVuZW15X2RtZwogICAgY2hhaW5zICs9IHJhbmRvbWl6ZShuYXZ5X3RvdGFsc1snY2hhaW5zJ10pCiAgcmV0dXJuIG15X2hwID4gMCBhbmQgaSA8IDUwCnByaW50IHN1bShtYXAoZmlnaHQsIHJhbmdlKHRyaWFscykpKSwgIi8iLCB0cmlhbHMg