%I #26 Dec 21 2023 01:15:34
%S 0,0,1,0,1,2,1,0,3,2,1,4,1,2,7,0,1,6,1,4,9,2,1,8,5,2,9,4,1,14,1,0,13,
%T 2,11,12,1,2,15,8,1,18,1,4,21,2,1,16,7,10,19,4,1,18,15,8,21,2,1,28,1,
%U 2,27,0,17,26,1,4,25,22,1,24,1,2,35,4,17,30,1,16,27,2,1,36,21,2,31,8,1,42,19
%N a(n) = n - phi(2*n), where phi() is the Euler totient A000010().
%D M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 840.
%D Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 24.
%H Antti Karttunen, <a href="/A176095/b176095.txt">Table of n, a(n) for n = 1..16384</a>
%H Antti Karttunen, <a href="/A176095/a176095.txt">Data supplement: n, a(n) computed for n = 1..65537</a>
%F a(n) = n - A062570(n).
%F a(2^k) = 0, k >= 0. - _Michel Lagneau_, Dec 17 2010
%F a(A000040(k)) = 1, k >= 2. - _Michel Lagneau_, Dec 17 2010, corrected by _Antti Karttunen_, May 19 2021
%F a(2^m*A000040(k)) = 2^m, m >= 1, k >= 2. - _Michel Lagneau_, Dec 17 2010
%F Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = 1 - 8/Pi^2 = 0.1894305... . - _Amiram Eldar_, Dec 21 2023
%e a(1) = 1 - phi(2) = 0;
%e a(2) = 2 - phi(2*2) = 2 - 2 = 0;
%e a(3) = 3 - phi(2*3) = 3 - 2 = 1;
%e If n = (2^m)*p, with m = 3 and p = 7, then n = 2^3 * 7 = 56, and a(56) = 2^3 = 8.
%p A176095 := proc(n)
%p n-numtheory[phi](2*n) ;
%p end proc:
%p seq(A176095(n),n=1..60) ;
%t Table[n-EulerPhi[2n],{n,0,100}] (* _Harvey P. Dale_, Jul 24 2011 *)
%o (PARI) A176095(n) = (n-eulerphi(n+n)); \\ _Antti Karttunen_, May 19 2021
%Y Cf. A000010, A000040, A051953, A062570.
%K nonn,easy
%O 1,6
%A _Michel Lagneau_, Apr 08 2010
%E Offset corrected; entry corrected and edited by _Michel Lagneau_, Apr 25 2010
|