ClearAll[ZumkellerQ] ZumkellerQ[n_] := Module[{d = Divisors[n], t, ds, x}, ds = Total[d]; If[Mod[ds, 2] == 1, False , t = CoefficientList[Product[1 + x^i, {i, d}], x]; t[[1 + ds/2]] > 0 ] ]; i = 1; res = {}; While[Length[res] < 220, r = ZumkellerQ[i]; If[r, AppendTo[res, i]]; i++; ]; res i = 1; res = {}; While[Length[res] < 40, r = ZumkellerQ[i]; If[r, AppendTo[res, i]]; i += 2; ]; res