login
A326184
a(n) = sigma(n) - A057521(n), where A057521 gives the powerful part of n, and sigma gives the sum of divisors of n.
3
0, 2, 3, 3, 5, 11, 7, 7, 4, 17, 11, 24, 13, 23, 23, 15, 17, 30, 19, 38, 31, 35, 23, 52, 6, 41, 13, 52, 29, 71, 31, 31, 47, 53, 47, 55, 37, 59, 55, 82, 41, 95, 43, 80, 69, 71, 47, 108, 8, 68, 71, 94, 53, 93, 71, 112, 79, 89, 59, 164, 61, 95, 95, 63, 83, 143, 67, 122, 95, 143, 71, 123, 73, 113, 99, 136, 95, 167, 79, 170, 40, 125, 83
OFFSET
1,2
LINKS
FORMULA
a(n) = A000203(n) - A057521(n).
PROG
(PARI)
A057521(n) = { my(f=factor(n)); prod(i=1, #f~, if(f[i, 2]>1, f[i, 1]^f[i, 2], 1)); }; \\ From A057521.
A326184(n) = (sigma(n)-A057521(n));
CROSSREFS
Cf. also A326187.
Sequence in context: A236165 A049855 A286868 * A326053 A296083 A294285
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 11 2019
STATUS
approved