OFFSET
0,2
EXAMPLE
a(0) = 1 because 1 cannot be written as the sum of two terms in A003601.
PROG
(Magma) f:=func<n|DivisorSigma(1, n) mod DivisorSigma(0, n) eq 0>; b:=[k: k in [1..1000] |f(n)]; a:=[]; for n in [0..60] do k:=1; while #RestrictedPartitions(k, 2, Set(b)) ne n do k:=k+1; end while; Append(~a, k); end for; a;
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, Jan 25 2024
STATUS
approved