OFFSET
1,1
COMMENTS
For squarefree numbers k, csigma(k) = k, where csigma(k) is the sum of the coreful divisors of k (A057723). Thus, if m is a term (csigma(m) > 2*m) and k is a squarefree number coprime to k, then csigma(k*m) = csigma(k) * csigma(m) = k * csigma(m) > 2*k*m, so k*m is a coreful abundant number. Therefore, the sequence of coreful abundant numbers (A308053) can be generated from this sequence by multiplying with coprime squarefree numbers. The asymptotic density of the coreful abundant numbers can be calculated from this sequence (see comment in A308053).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
72 is a term since csigma(72) = 168 > 2 * 72, and 72 = 2^3 * 3^2 is powerful.
MATHEMATICA
f[p_, e_] := (p^(e+1)-1)/(p-1)-1; s[1] = 1; s[n_] := If[AllTrue[(fct = FactorInteger[n])[[;; , 2]], #>1 &], Times @@ f @@@ fct, 0]; seq={}; Do[If[s[n] > 2*n, AppendTo[seq, n]], {n, 1, 5000}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Sep 02 2022
STATUS
approved