%I #17 Apr 22 2019 13:49:22
%S 0,1,1,-1,1,6,1,-5,-5,8,1,0,1,10,9,-13,1,-6,1,-2,11,14,1,-12,-19,16,
%T -23,-4,1,42,1,-29,15,20,13,-24,1,22,17,-22,1,54,1,-8,-21,26,1,-36,
%U -41,-32,21,-10,1,-42,17,-32,23,32,1,12,1,34,-31,-61,19,78,1,-14,27,74,1,-60,1,40,-51,-16,19,90,1,-62,-77,44,1,12,23,46
%N a(n) = A048250(n) - n, where A048250(n) is the sum of squarefree divisors of n.
%H Antti Karttunen, <a href="/A325313/b325313.txt">Table of n, a(n) for n = 1..16384</a>
%H Antti Karttunen, <a href="/A325313/a325313.txt">Data supplement: n, a(n) computed for n = 1..65537</a>
%F a(n) = A048250(n) - n.
%F a(n) = A325314(n) - A033879(n).
%F a(A228058(n)) = -A325319(n).
%t Array[DivisorSum[#, # &, SquareFreeQ] - # &, 86] (* _Michael De Vlieger_, Apr 21 2019 *)
%o (PARI)
%o A048250(n) = factorback(apply(p -> p+1,factor(n)[,1]));
%o A325313(n) = (A048250(n) - n);
%Y Cf. A033879, A048250, A126795, A228058, A325314, A325315, A325319, A325320.
%K sign
%O 1,6
%A _Antti Karttunen_, Apr 21 2019