login
A334405
Pseudoperfect numbers k such that there is a subset of divisors of k whose sum is 2*k and for each d in this subset k/d is also in it.
4
6, 28, 36, 60, 84, 90, 120, 156, 210, 216, 240, 252, 270, 300, 312, 330, 336, 352, 396, 420, 468, 480, 496, 504, 540, 546, 552, 576, 588, 594, 600, 616, 624, 630, 648, 660, 672, 714, 720, 756, 760, 780, 784, 792, 816, 840, 864, 888, 900, 924, 960, 972, 1000
OFFSET
1,1
COMMENTS
Includes all the perfect numbers (A000396).
The McCormack and Zelinsky preprint shows that no terms are 2 (mod 3), and also that no terms are 3 (mod 4). That paper also asks if there are infinitely many odd terms. Empirically, odd terms are much rarer than even terms. - Joshua Zelinsky, Feb 28 2024
LINKS
Tim McCormack and Joshua Zelinsky, Weighted Versions of the Arithmetic-Mean-Geometric Mean Inequality and Zaremba's Function, arXiv:2312.11661 [math.NT], 2023. Mentions this sequence.
FORMULA
36 is a term since {1, 2, 3, 12, 18, 36} is a subset of its divisors whose sum is 72 = 2 * 36, and for each divisor d in this subset 36/d is also in it: 1 * 36 = 2 * 18 = 3 * 12 = 36.
MATHEMATICA
seqQ[n_] := Module[{d = Divisors[n]}, nd = Length[d]; divpairs = If[EvenQ[nd], d[[1 ;; nd/2]] + d[[-1 ;; nd/2 + 1 ;; -1]], Join[d[[1 ;; (nd - 1)/2]] + d[[-1 ;; (nd + 3)/2 ;; -1]], {d[[(nd + 1)/2]]}]]; SeriesCoefficient[Series[Product[1 + x^divpairs[[i]], {i, Length[divpairs]}], {x, 0, 2*n}], 2*n] > 0]; Select[Range[1000], seqQ]
CROSSREFS
Subsequence of A005835.
A000396 is a subsequence.
Sequence in context: A272971 A117948 A252234 * A376997 A242344 A344588
KEYWORD
nonn
AUTHOR
Amiram Eldar, Apr 27 2020
STATUS
approved