OFFSET
1,1
COMMENTS
The exponential version of A171641.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
MATHEMATICA
dQ[n_, m_] := (n>0&&m>0 &&Divisible[n, m]); expDivQ[n_, d_] := Module[ {ft = FactorInteger[n]}, And@@MapThread[dQ, {ft[[;; , 2]], IntegerExponent[ d, ft[[;; , 1]]]} ]]; ediv[n_] := Module[ {d=Rest[Divisors[n]]}, Select[ d, expDivQ[n, #]&]]; esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[ Last[#]]}] &) /@ FactorInteger[n]; seq={}; Do[s=esigma[n]; If[OddQ[s] || s<=2n, Continue[]]; div = ediv[n]; If[Coefficient[Times @@ (1 + x^div) // Expand, x, s/2] == 0, AppendTo[seq, n]], {n, 1, 10000}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jan 11 2019
STATUS
approved