from random import *

def est_dope(p):
    if random() < p:
        return 1
    else:
        return 0

def test_positif(dope):
    if dope == 1:
        if random() < 0.94:
            return 1
    else:
        if random() < 0.04:
            return 1
    return 0

def simulations(p):
    nb_dopes = 0
    nb_tests_positifs = 0
    nb_dopes_positifs = 0
    
    for sportif in range(2500):
        dope = est_dope(p)
        positif = test_positif(dope)
        dope_positif = dope*positif
        
        nb_dopes = nb_dopes + dope
        nb_tests_positifs = nb_tests_positifs + positif
        nb_dopes_positifs = nb_dopes_positifs + dope*positif
    
    return nb_dopes_positifs/nb_tests_positifs