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”).

A326136
a(n) = sigma(n) - sigma(A028234(n)), where sigma is the sum of divisors of n, and A028234 gives n without any occurrence of its smallest prime factor.
3
0, 2, 3, 6, 5, 8, 7, 14, 12, 12, 11, 24, 13, 16, 18, 30, 17, 26, 19, 36, 24, 24, 23, 56, 30, 28, 39, 48, 29, 48, 31, 62, 36, 36, 40, 78, 37, 40, 42, 84, 41, 64, 43, 72, 72, 48, 47, 120, 56, 62, 54, 84, 53, 80, 60, 112, 60, 60, 59, 144, 61, 64, 96, 126, 70, 96, 67, 108, 72, 96, 71, 182, 73, 76, 93, 120, 84, 112, 79, 180, 120, 84, 83
OFFSET
1,2
FORMULA
a(n) = A000203(n) - A000203(A028234(n)).
From Amiram Eldar, Dec 21 2024: (Start)
a(n) = A000203(n) - A326135(n).
Sum_{k=1..n} a(k) ~ (zeta(2)/2) * (1 - c) * n^2, where c is defined in the corresponding formula in A326135. (End)
PROG
(PARI)
A028234(n) = { my(f = factor(n)); if (#f~, f[1, 1] = 1); factorback(f); }; \\ From A028234
A326136(n) = (sigma(n) - sigma(A028234(n)));
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Antti Karttunen, Jun 08 2019
STATUS
approved