def B(n,p):
    return sum([Bernoulli(p) for k in range(n)])