[go: up one dir, main page]

login
A373058 revision #5

A373058
The sum of the aliquot coreful divisors of the nonsquarefree numbers.
1
2, 6, 3, 6, 14, 6, 10, 18, 5, 12, 14, 30, 36, 30, 22, 15, 42, 7, 10, 26, 24, 42, 30, 21, 62, 34, 96, 15, 38, 70, 39, 42, 66, 30, 46, 90, 14, 33, 80, 78, 126, 98, 58, 39, 90, 11, 62, 30, 42, 126, 66, 60, 102, 70, 216, 21, 74, 30, 114, 51, 78, 150, 78, 82, 126, 13
OFFSET
1,1
COMMENTS
A coreful divisor d of n is a divisor that is divisible by every prime that divides n (see also A307958).
The positive terms of A336563: if k is a squarefree number (A005117) then the only coreful divisor of k is k itself, so k has no aliquot coreful divisors.
The number of the aliquot coreful divisors of the n-th nonsquarefree number is A368039(n).
FORMULA
a(n) = A336563(A013929(n)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = (A065487 - 1)/(1-1/zeta(2))^2 = 1.50461493205911656114... .
MATHEMATICA
f[p_, e_] := (p^(e + 1) - 1)/(p - 1) - 1; s[1] = 0; s[n_] := Times @@ f @@@ FactorInteger[n] - n; Select[Array[s, 300], # > 0 &]
PROG
(PARI) lista(kmax) = {my(f); for(k = 1, kmax, f = factor(k); if(!issquarefree(f), print1(prod(i = 1, #f~, (f[i, 1]^(f[i, 2]+1) - 1)/(f[i, 1] - 1) - 1) - k, ", "))); }
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, May 21 2024
STATUS
proposed