login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A379029
Modified exponential abundant numbers: numbers k such that A241405(k) > 2*k.
3
30, 42, 66, 70, 78, 102, 114, 120, 138, 150, 168, 174, 186, 210, 222, 246, 258, 270, 282, 294, 318, 330, 354, 366, 390, 402, 420, 426, 438, 462, 474, 498, 510, 534, 546, 570, 582, 606, 618, 630, 642, 654, 660, 678, 690, 714, 726, 750, 762, 770, 780, 786, 798, 822
OFFSET
1,1
COMMENTS
All the squarefree abundant numbers (A087248) are terms since A241405(k) = A000203(k) for a squarefree number k.
If k is a term and m is coprime to k them k*m is also a term.
The numbers of terms that do no exceed 10^k, for k = 2, 3, ..., are 5, 67, 767, 7595, 76581, 764321, 7644328, 76468851, 764630276, ... . Apparently, the asymptotic density of this sequence exists and equals 0.07646... .
LINKS
MATHEMATICA
f[p_, e_] := DivisorSum[e + 1, p^(# - 1) &]; mesigma[1] = 1; mesigma[n_] := Times @@ f @@@ FactorInteger[n]; meAbQ[n_] := mesigma[n] > 2*n; Select[Range[1000], meAbQ]
PROG
(PARI) is(n) = {my(f=factor(n)); prod(i=1, #f~, sumdiv(f[i, 2]+1, d, f[i, 1]^(d-1))) > 2*n; }
CROSSREFS
Subsequence of A005101.
Subsequences: A034683, A087248, A379030, A379031.
Similar sequences: A064597, A129575, A129656, A292982, A348274, A348604.
Sequence in context: A053858 A075819 A306217 * A034683 A328328 A302574
KEYWORD
nonn,easy,new
AUTHOR
Amiram Eldar, Dec 14 2024
STATUS
approved