for i in range(33333, 55555+1): sum_del = 0 for d1 in range(2, int(i**0.5)+1): if i%d1 == 0: d1_norm = True for j in range(2, int(d1**0.5)+1): if d1%j == 0: d1_norm = False d2 = i//d1 d2_norm = True for j in range(2, int(d2**0.5)+1): if d2%j == 0: d2_norm = False if d2 == d1 and d1_norm and d2_norm: sum_del += d1 else: if d1_norm: sum_del += d1 elif d2_norm: sum_del += d2 if sum_del > 250 and i%sum_del == 0: print(i, sum_del)