from random import random
from math import ceil, log

def exponential(p):
    y = random()
    x = ceil( -log(1-y) / p )
    return(x)


p = 1/6
r = 3

for i in range(0,5):
    y = 0
    for j in range(0,r):
        y += exponential(p)
    print(i, y)