Euler #53

from math import factorial
def combination(n, r):
    return int(factorial(n)/(factorial(r)*factorial(n-r)))

answer = 0
nums = []
for a in range(1, 101):
    print(a)
    nums.append(a)
    for b in range(1, a):
        if combination(a, b) > 10**6:
            answer += 1
print(answer)