import numpy.random as rd

n = 5

# Simultation d'une partie

def partie(n) :
    L = rd.binomial(1,.5,n)
    s = 0
    for k in L :
        s += k
    if s == 1 or s == n-1 :
        return 1
    else :
        return 0
        
for k in range(50) :
    print(partie(n))

# temps d'attente d'une victoire
        
def X(n) :
    k = 1
    while partie(n) == 0 :
        k += 1
    return k
    
for k in range(50) :
    print(X(n))

# Approximation des espértances et variances
    
def E(n) :
    s = 0
    for k in range(500) :
        s += X(n)
    return s/500
    
def V(n) :
    s = 0
    for k in range(500) :
        s += X(n)**2
    return s/500 - E(n)**2
    
for k in range(10) :
    print(E(n),V(n))