OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..3000
EXAMPLE
48 is a term since it is not exponentially odd number (48 = 2^4 * 3 and 4 is even), so not all of its divisors are bi-unitary, and it is the sum of a subset of its bi-unitary divisors: 8 + 16 + 24 = 48.
MATHEMATICA
f[n_] := Select[Divisors[n], Function[d, CoprimeQ[d, n/d]]]; bdiv[m_] := Select[Divisors[m], Last@Intersection[f@#, f[m/#]] == 1 &]; bPspQ[n_] := Module[{d = Most @ bdiv[n], x}, SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n] > 0]; expOddQ[n_] := AllTrue[Last /@ FactorInteger[n], OddQ]; Select[Range[1000], ! expOddQ[#] && bPspQ[#] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jun 30 2020
STATUS
approved
