%I #4 Jul 19 2021 18:30:38

%S 1,1,2,2,4,1,6,5,7,3,10,4,12,5,7,12,16,11,18,8,11,9,22,13,23,11,24,12,

%T 28,-1,30,27,19,15,23,26,36,17,23,23,40,1,42,20,32,21,46,34,47,41,31,

%U 24,52,45,39,33,35,27,58,14,60,29,46,58,47,5,66,32,43,11,70,59,72,35,64,36,59,7,78,58,77,39,82,22,63,41

%N a(n) = n - A342001(n).

%F a(n) = n - A342001(n) = n - (A003415(n) / A003557(n)).

%o (PARI)

%o A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));

%o A003557(n) = (n/factorback(factorint(n)[, 1]));

%o A342001(n) = (A003415(n) / A003557(n));

%o A346244(n) = (n - A342001(n));

%Y Cf. A003415, A003557, A342001, A346245 (positions of negative terms).

%K sign

%O 1,3

%A _Antti Karttunen_, Jul 19 2021