4721.sagews

chances = range(97, 52, -2) + range(46, 39, -2) + range(28, 23, -2)
chances = range(97, 80, -2) + [74,72,70, 58,56,54]
chances = [float(c) / 100 for c in chances]

skills = [[31184, 15876, 189],     [51975, 26460, 315],     [72765, 34044, 441],     [93555, 47628, 567],     [187425, 28035, 315],   [262395, 39249, 441],     [337365, 50463, 567]]

skill = skills[3]

def cost(k, start, delta, acc):
    return start + (k//3) * delta + (k//3 - 1) * (k // 3 // 2) * acc

def f(n):
    if n == 0:
        return vector([0, 0, 0])
    return (f(n-1) + vector([1 if n == 1 else 0, cost(n-1, skill[0], skill[1], skill[2]), 1])) / chances[n-1]

def g(s, f):
    if f == s:
        return vector([0, 0, 0])
    else:
        return g(s, f-1) + vector([1, 5*cost(f-1, skill[0], skill[1], skill[2]), 1]) / chances[f-1]

def print_ench(n):
    for i in xrange(n, -1, -1):
        vf, vg = f(i), g(i, n)
        print 'Safe from %d: %.1f usual, %.1f safe, %.1f adena, %.1f clicks' % (i, vf[0], vg[0], vf[1]+vg[1], vf[2]+vg[2])

print_ench(30)

Safe from 30: 1822188.3 usual, 0.0 safe, 2268213750808.1 adena, 14764572.8 clicks Safe from 29: 437325.2 usual, 4.2 safe, 544382017155.4 adena, 3543500.6 clicks Safe from 28: 113704.6 usual, 8.0 safe, 141557505886.0 adena, 921316.1 clicks Safe from 27: 31837.3 usual, 11.6 safe, 39660797470.2 adena, 257976.8 clicks Safe from 26: 12734.9 usual, 14.1 safe, 15888737902.3 adena, 103199.2 clicks Safe from 25: 5348.7 usual, 16.5 safe, 6700326738.4 adena, 43353.2 clicks Safe from 24: 2353.4 usual, 18.7 safe, 2977456316.8 adena, 19085.9 clicks Safe from 23: 1082.6 usual, 20.9 safe, 1400318375.3 adena, 8790.8 clicks Safe from 22: 573.8 usual, 22.8 safe, 771044085.2 adena, 4669.8 clicks Safe from 21: 315.6 usual, 24.6 safe, 453583634.2 adena, 2579.5 clicks Safe from 20: 179.9 usual, 26.4 safe, 288083791.3 adena, 1481.6 clicks Safe from 19: 106.1 usual, 28.1 safe, 199556243.2 adena, 885.7 clicks Safe from 18: 64.7 usual, 29.7 safe, 151188307.6 adena, 551.9 clicks Safe from 17: 40.8 usual, 31.3 safe, 124057426.5 adena, 359.2 clicks Safe from 16: 26.5 usual, 32.8 safe, 108866823.9 adena, 245.0 clicks Safe from 15: 17.8 usual, 34.3 safe, 100462038.9 adena, 175.5 clicks Safe from 14: 12.3 usual, 35.8 safe, 95704344.4 adena, 132.2 clicks Safe from 13: 8.7 usual, 37.2 safe, 93294599.5 adena, 104.6 clicks Safe from 12: 6.4 usual, 38.6 safe, 92312597.6 adena, 86.8 clicks Safe from 11: 4.8 usual, 39.9 safe, 91930384.0 adena, 75.1 clicks Safe from 10: 3.7 usual, 41.2 safe, 92097496.8 adena, 67.3 clicks Safe from 9: 2.9 usual, 42.4 safe, 92613380.0 adena, 62.1 clicks Safe from 8: 2.3 usual, 43.7 safe, 93101947.7 adena, 58.6 clicks Safe from 7: 1.9 usual, 44.9 safe, 93729710.2 adena, 56.2 clicks Safe from 6: 1.7 usual, 46.1 safe, 94445671.1 adena, 54.7 clicks Safe from 5: 1.4 usual, 47.2 safe, 94986798.8 adena, 53.7 clicks Safe from 4: 1.3 usual, 48.3 safe, 95559281.9 adena, 53.1 clicks Safe from 3: 1.2 usual, 49.4 safe, 96148647.9 adena, 52.8 clicks Safe from 2: 1.1 usual, 50.5 safe, 96535981.0 adena, 52.7 clicks Safe from 1: 1.0 usual, 51.6 safe, 96924820.5 adena, 52.6 clicks Safe from 0: 0.0 usual, 52.6 safe, 97310614.3 adena, 52.6 clicks