OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
fun[p_, e_] := Module[{b = IntegerDigits[e, 2]}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; isigma[1] = 1; isigma[n_] := Times @@ fun @@@ FactorInteger[n]; infabQ[n_] := isigma[n] > 2*n; idivs[x_] := If[x == 1, 1, Sort @ Flatten @ Outer[Times, Sequence @@ (FactorInteger[x] /. {p_, m_Integer} :> p^Select[Range[0, m], BitOr[m, #] == m &])]]; infwQ[n_] := infabQ[n] && Module[{d = Most @ idivs[n]}, SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n] == 0]; pow2Q[n_] := n == 2^IntegerExponent[n, 2]; seq = {}; Do[If[!pow2Q[DivisorSigma[0, n]] && infwQ[n], AppendTo[sm n]], {n, 1, 10^5}]; s
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jun 30 2020
EXTENSIONS
More terms from Amiram Eldar, Mar 25 2023
STATUS
approved