27 lines
435 B
Text
27 lines
435 B
Text
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
|