|
%I #12 Mar 19 2018 22:14:15
%S 0,1,1,0,1,2,1,4,3,2,1,4,1,0,5,0,1,6,1,14,3,8,1,8,3,12,9,16,1,14,1,10,
%T 11,8,3,12,1,0,15,-4,1,12,1,8,15,2,1,10,1,8,11,4,1,18,11,24,3,26,1,28,
%U 1,30,9,30,15,20,1,32,5,46,1,24,1,36,7,40,9,24,1,50,27,38,1,60,11,36,29,32,1,42,13,32,33,26,3,50,1
%N Difference between A244042 and its Möbius transform.
%H Antti Karttunen, <a href="/A300823/b300823.txt">Table of n, a(n) for n = 1..65537</a>
%F a(n) = A244042(n) - A300821(n).
%F a(n) = -Sum_{d|n, d<n} A008683(n/d)*A244042(d) = Sum_{d|n, d<n} A300821(d).
%F a(n) = A051953(n) - A300824(n).
%t f[n_] := FromDigits[IntegerDigits[n, 3] /. 2 -> 0, 3]; Table[f@ n - DivisorSum[n, MoebiusMu[n/#] f@ # &], {n, 97}] (* _Michael De Vlieger_, Mar 17 2018 *)
%o (PARI)
%o A244042(n) = fromdigits(apply(x->(x%2), digits(n, 3)), 3);
%o A300823(n) = -sumdiv(n,d,(d<n)*moebius(n/d)*A244042(d));
%Y Cf. A008683, A051953, A244042, A300821, A300824, A300825.
%K sign,base
%O 1,6
%A _Antti Karttunen_, Mar 14 2018
|
