n = 8
p = 0.25

print(n*p)
print(sqrt(n*p*(1-p)))

print(c(dbinom(n, n, p), sum(dbinom(3:n, n, p))))