OFFSET
1,1
COMMENTS
EXAMPLE
24 is in the sequence since its aliquot infinitary divisors are 1, 2, 3, 4, 6, 8, 12 and 24 and 1 + 2 + 3 + 4 + 6 + 8 = 24.
MATHEMATICA
idivs[x_] := If[x == 1, 1, Sort@ Flatten@ Outer[Times, Sequence @@ (FactorInteger[ x ] /. {p_, m_Integer} :> p^Select[Range[0, m], BitOr[m, #] == m &])]]; subtr = If[#1 < #2, Throw[#1], #1 - #2] &; selDivs[n_] := Catch@Fold[subtr, n, Drop[idivs[n], -2]]; s= {}; Do[If[selDivs[n] == 0, AppendTo[s, n]], {n, 2, 10^6}]; s(* after Alonso del Arte at A194472 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Sep 14 2019
STATUS
approved