def Seuil(p): n=1 S=0.5 while 1-S>=10**(-p): n=n+1 S=S+0.5**n return n,S