from probastat import binom


def E(n):
    X = binom(n,0.95)
    s = 100*n + X.proba_cdf(100)
    
    for k in range(1,n-99):
        s += (100*n-250*k)*X.proba(100+k)
        
    return s

n = 101
M = 10000

while True:
    if E(n) > M:
        M = E(n)
        n += 1
        print(n,M)
    else:
        print("Valeur maximale atteinte pour n =", n)
        break