login
a(n) = phi(A003401(n)).
2

%I #15 Sep 25 2024 00:37:34

%S 1,1,2,2,4,2,4,4,4,8,8,16,8,8,8,16,16,16,16,32,16,32,32,32,64,32,32,

%T 32,64,64,64,64,64,64,64,128,128,256,128,128,128,128,128,128,128,256,

%U 256,256,256,256,256,512,256,256,256,512,512,512,512,1024,512,512,512,512

%N a(n) = phi(A003401(n)).

%C All terms are powers of 2. - _Jianing Song_, Sep 28 2018

%H Amiram Eldar, <a href="/A094269/b094269.txt">Table of n, a(n) for n = 1..1000</a>

%F a(n) = 2^A319821(n). - _Amiram Eldar_, Sep 25 2024

%t Do[If[IntegerQ[Log[2, EulerPhi[n]]], Print[n];ta[[u]]=n;u=u+1], {n, 1, 10000}] EulerPhi[ta]

%o (PARI) for(n=1, 1000, my(i=eulerphi(n)); if(omega(2*i)==1, print1(i, “, “))) \\ _Jianing Song_, Sep 28 2018

%Y Cf. A000010, A003401, A319821.

%K nonn

%O 1,3

%A _Labos Elemer_, Jun 02 2004