OFFSET
1,1
COMMENTS
Differs from A054979 first at a(44), since 4900 is in this sequence but not in A054979. - R. J. Mathar, Jun 02 2020
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
36 is a term since there is a single way in which 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]] == 2]]; Select[Range[10^4], ezQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 27 2020
STATUS
approved