from math import *
from random import *

def S(n,p):
    s=0
    for i in range(n):
        a=random()
        if a<=p:
            s=s+1
    return s

def P(n,p,N):
    b=0
    for i in range(N):
        if abs(S(n,p)-n*p)>=sqrt(n):
            b=b+1
    p=b/N
    return p