OFFSET
1,2
EXAMPLE
MATHEMATICA
seq[len_, kmax_] := Module[{s = Table[0, {len}], c = 0, k = 1, ind}, While[k < kmax && c < len, ind = DivisorSum[k, 1 &, DivisorSigma[0, DivisorSigma[1, #]] == DivisorSigma[0, #] &]; If[ind <= len && s[[ind]] == 0, c++; s[[ind]] = k]; k++]; s]; seq[20, 10^6] (* Amiram Eldar, May 05 2023 *)
PROG
(Magma) f:=func<n|#Divisors(DivisorSigma(1, n)) eq #Divisors(n)>; a:=[]; for n in [1..37] do k:=1; while #[d:d in Divisors(k)|f(d)] ne n do k:=k+1; end while; Append(~a, k); end for; a;
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, May 03 2023
STATUS
approved