Modeling3

import random


def get_random():
    x_ = random.uniform(0, 1)
    result = 0
    while x_ > 0:
        result += 1
        x_ -= 0.25
    return result


counts = []
tmp = {1: 0, 2: 0, 3: 0, 4: 0}
n = 0
for _ in range(0, 1000):
    n += 1
    x = get_random()
    tmp[x] += 1
    if 3 in tmp.values():
        counts.append(n)
        tmp = {1: 0, 2: 0, 3: 0, 4: 0}
        n = 0


print("Мат ожидание: {}".format(float(sum(counts)) / max(len(counts), 1)))