A326184 a(n) = sigma(n) - A057521(n), where A057521 gives the powerful part of n, and sigma gives the sum of divisors of n.

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

Antti Karttunen, Table of n, a(n) for n = 1..20000

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. A000203, A057521, A326185.

Cf. also A326187.

Sequence in context: A236165 A049855 A286868 * A326053 A296083 A294285

Adjacent sequences: A326181 A326182 A326183 * A326185 A326186 A326187

Keyword

nonn

Author

Antti Karttunen, Jul 11 2019

Status

approved