facs = [1]*101 for i in range(1, 101): facs[i] = facs[i-1] * i count = 0 for n in range(23, 101): for r in range(n+1): x = facs[n] / (facs[r]*facs[n-r]) if x > 1000000: count += 1 print count