orobability of controlling >50% of an N-size committee

committee_takeover.py

from decimal import Decimal
import math

def calculate_summation(N, p):
    p = Decimal(p)
    result = Decimal(0)
    for k in range(math.floor(N/2), N+1):
        term = Decimal(math.comb(N, k)) * (p**k) * ((1-p)**(N-k))
        result += term
    return float(result)

# Testing the function
N = 9000
p = 0.40

result = calculate_summation(N, p)
print(result)