Day 4, Statistics Binomial Distribution 2.py

# Enter your code here. Read input from STDIN. Print output to STDOUT
p,n=map(float,raw_input().split())
p/=100
q=1-p
def factorial(n):
    if n==0 or n==1:
        return 1
    else:
        return n*factorial(n-1)
    
def nCr(n,r):
    return factorial(n)/(factorial(r) * factorial(n-r))

firstSum=0
for i in range(0,3):
    ncr=nCr(n,i)
    temp=(p**i) * (q** (n-i))
    firstSum+=ncr*temp
    
print round(firstSum,3)

secondSum=0
for i in range(2,11):
    ncr=nCr(n,i)
    temp=(p**i) * (q** (n-i))
    secondSum+=ncr*temp

print round(secondSum,3)