import numpy as np
from numpy.random import choice
from copy import copy
tier1 = [3185, 2937, 2897, 2670]
tier2 = [2536, 2397, 2153, 1976]
tier3 = [1914, 1573, 1369, 1333]
tier4 = [1315, 1265, 1240, 1237]
tier5 = [1172, 1106, 752, 740]

sims = []
for i in range(10000):
    tiers = {1: copy(tier1),
             2: copy(tier2),
             3: copy(tier3),
             4: copy(tier4),
             5: copy(tier5)}
    grps = {}
    for grp in ['a', 'b', 'c', 'd']:
        choices = {1: choice(tiers[1]),
                   2: choice(tiers[2]),
                   3: choice(tiers[3]),
                   4: choice(tiers[4]),
                   5: choice(tiers[5])}
        grps[grp] = list(choices.values())
        for x in choices:
            tiers[x].remove(choices[x])
    means = [np.mean(tuple(v)) for v in grps.values()]
    max_pts = max(means)
    means.remove(max_pts)
    avg = np.mean(means)
    sims.append(max_pts - avg)

diff = 2011 - np.mean([1677, 1685, 1779])
numerator = 0
for x in sims:
    if x >= diff:
        numerator += 1

print(numerator)