OFFSET
1,1
COMMENTS
EXAMPLE
24 is in the sequence since its aliquot bi-unitary divisors are 1, 2, 3, 4, 6, 8, 12 and 24 and 1 + 2 + 3 + 4 + 6 + 8 = 24.
MATHEMATICA
f[n_] := Select[Divisors[n], Function[d, CoprimeQ[d, n/d]]]; bdiv[m_] := Select[Divisors[m], Last@Intersection[f@#, f[m/#]] == 1 &]; subtr = If[#1 < #2, Throw[#1], #1 - #2] &; selDivs[n_] := Catch@Fold[subtr, n, Drop[bdiv[n], -2]]; a = {}; Do[ If[selDivs[n] == 0, AppendTo[a, n]; Print[n]], {n, 2, 10^6}]; a (* after Alonso del Arte at A194472 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Oct 13 2017
STATUS
approved