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

Author

Antti Karttunen, Jul 11 2019

Keywords

• nonn

Links

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

Crossrefs

Cf. A000203, A057521, A326185.

Cf. also A326187.

Programs

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));

Formula

a(n) = A000203(n) - A057521(n).