OFFSET
1,1
COMMENTS
First differs from A318100 at n = 49: 4900 is a term that is not an exponential pseudoperfect number.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
36 is a term since its exponential divisors, {6, 12, 18, 36}, can be partitioned into 2 disjoint sets whose sum is equal: 6 + 12 + 18 = 36.
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]]]}]]; eDivs[n_] := Module[{d = Rest[Divisors[n]]}, Select[d, expDivQ[n, #] &]]; ezQ[n_] := Module[{d = eDivs[n], sum, x}, sum = Plus @@ d; If[sum < 2*n || OddQ[sum], False, CoefficientList[Product[1 + x^i, {i, d}], x][[1 + sum/2]] > 0]]; Select[Range[10^4], ezQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 27 2020
STATUS
approved