OFFSET
1,1
COMMENTS
LINKS
Peter Luschny, Table of n, a(n) for n = 1..10000
EXAMPLE
The Zumkeller divisors of 80 are {20, 40, 80}, so 80 is a term.
The divisors of 81 are {1, 3, 9, 27, 81}, none of which is Zumkeller, so 81 is not a term.
MAPLE
with(NumberTheory):
isZumkeller := proc(n) option remember; local s, p, i, P;
s := SumOfDivisors(n);
if s::odd or s < n*2 then false else
P := mul(1 + x^i, i in Divisors(n));
is(0 < coeff(P, x, s/2)) fi end:
select(n -> ormap(isZumkeller, Divisors(n)), [seq(1..270)]);
MATHEMATICA
znQ[n_]:=Length[Select[{#, Complement[Divisors[n], #]}&/@Most[Rest[ Subsets[ Divisors[ n]]]], Total[#[[1]]]==Total[#[[2]]]&]]>0; zn=Select[Range[300], znQ] (* zn from A083207 *) ; Select[Range[270], IntersectingQ[Divisors[#], zn]&] (* James C. McMahon, Oct 23 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Oct 20 2024
EXTENSIONS
Incorrect comment removed by Peter Luschny, Dec 02 2024
STATUS
approved