A308053 Coreful abundant numbers: numbers k such that csigma(k) > 2*k, where csigma(k) is the sum of the coreful divisors of k (A057723).

72, 108, 144, 200, 216, 288, 324, 360, 400, 432, 504, 540, 576, 600, 648, 720, 756, 784, 792, 800, 864, 900, 936, 972, 1000, 1008, 1080, 1152, 1188, 1200, 1224, 1296, 1368, 1400, 1404, 1440, 1512, 1568, 1584, 1600, 1620, 1656, 1728, 1764, 1800, 1836, 1872, 1936

Offset

1,1

Comments

Analogous to A005101 as A307958 is analogous to A000396.

The asymptotic density of this sequence is Sum_{n>=1} f(A356871(n)) = 0.0262215..., where f(n) = (6/(Pi^2*n))*Product_{prime p|n}(p/(p+1)). - Amiram Eldar, Sep 02 2022

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Example

72 is in the sequence since its coreful divisors are 6, 12, 18, 24, 36, 72, whose sum is 168 > 2 * 72.

Mathematica

f[p_, e_] := (p^(e+1)-1)/(p-1)-1; a[1]=1; a[n_] := Times @@ (f @@@ FactorInteger[n]); s={}; Do[If[a[n] > 2n, AppendTo[s, n]], {n, 1, 2000}]; s

Prog

(PARI) rad(n) = factorback(factorint(n)[, 1]); \\ A007947
s(n) = rad(n)*sigma(n/rad(n)); \\ A057723
isok(k) = s(k) > 2*k; \\ Michel Marcus, May 11 2019
(PARI) isok(k) = {my(f=factor(k)); prod(i = 1, #f~, (f[i, 1]^(f[i, 2]+1)-1) / (f[i, 1]-1)-1) > 2*k}; \\ Amiram Eldar, Sep 02 2022

Crossrefs

Cf. A005101, A057723, A307958.

A339940 and A356871 are subsequences.

Subsequence of A129575.

Sequence in context: A102562 A357460 A216426 * A356871 A359280 A307758

Adjacent sequences: A308050 A308051 A308052 * A308054 A308055 A308056

Keyword

nonn

Author

Amiram Eldar, May 10 2019

Status

approved