|
%I #6 Jul 12 2019 07:21:12
%S 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,
%T 13,52,29,71,31,31,47,53,47,55,37,59,55,82,41,95,43,80,69,71,47,108,8,
%U 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
%N a(n) = sigma(n) - A057521(n), where A057521 gives the powerful part of n, and sigma gives the sum of divisors of n.
%H Antti Karttunen, <a href="/A326184/b326184.txt">Table of n, a(n) for n = 1..20000</a>
%F a(n) = A000203(n) - A057521(n).
%o (PARI)
%o A057521(n) = { my(f=factor(n)); prod(i=1, #f~, if(f[i, 2]>1, f[i, 1]^f[i, 2], 1)); }; \\ From A057521.
%o A326184(n) = (sigma(n)-A057521(n));
%Y Cf. A000203, A057521, A326185.
%Y Cf. also A326187.
%K nonn
%O 1,2
%A _Antti Karttunen_, Jul 11 2019
|
