OFFSET
1,1
COMMENTS
All the squarefree abundant numbers (A087248) are also modified exponential abundant numbers (A379029). This sequence lists the terms of A379029 that are not in A087248.
The numbers of terms that do no exceed 10^k, for k = 3, 4, ..., are 14, 211, 2090, 21236, 212744, 2123071, 21235175, 212450318, ... . Apparently, the asymptotic density of this sequence exists and equals 0.0212... .
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := DivisorSum[e + 1, p^(# - 1) &]; mesigma[1] = 1; mesigma[n_] := Times @@ f @@@ FactorInteger[n]; nsmeAbQ[n_] := !SquareFreeQ[n] && mesigma[n] > 2*n; Select[Range[2400], nsmeAbQ]
PROG
(PARI) is(n) = {my(f=factor(n)); if(issquarefree(f), 0, prod(i=1, #f~, sumdiv(f[i, 2]+1, d, f[i, 1]^(d-1))) > 2*n); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Dec 14 2024
STATUS
approved