OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
48 is in the sequence since its bi-unitary divisors are 1, 2, 3, 6, 8, 16, 24, 48 and 48 = 8 + 16 + 24.
MATHEMATICA
f[n_] := Select[Divisors[n], Function[d, CoprimeQ[d, n/d]]]; bdiv[m_] := Select[Divisors[m], Last@Intersection[f@#, f[m/#]] == 1 &]; a = {}; n = 0; While[n < 1000, n++; d = Most[bdiv[n]]; c = SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n]; If[c > 0, AppendTo[a, n]]]; a (* after T. D. Noe at A005835 and Michael De Vlieger at A188999 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Sep 27 2017
STATUS
approved