import itertools import math vs = [1, 2, 4, 8, 16, 3, 6, 12, 24, 5, 10, 15, 20] for a in [1, 2, 4, 8, 16]: for b in [1, 3, 9, 27]: for c in [1, 5, 25]: vs.append(a * b * c) vs = sorted(set(vs)) # print(vs) for sz in range(10, len(vs) + 1): # print(sz) for c in itertools.combinations(vs, sz): l = list(c) tgt = sum(l) - math.lcm(*l) # print(l, tgt) if tgt in l: print(len(l)) print(*l) exit(0)