limit = 12000 fractions = set() for d in xrange(2, limit+1): start_n = (d+3)/3 end_n = d/2 for n in xrange(start_n, end_n+1): fractions.add(float(n)/d) print len(fractions)