login
Product_{d|n, A007431(d) > 0} prime(A007431(d)), where A007431 is the Möbius transform of Euler's totient function.
3

%I #5 Sep 05 2018 17:17:21

%S 2,2,4,4,10,4,22,12,28,10,46,16,62,22,100,84,94,28,118,100,484,46,146,

%T 144,530,62,1036,484,206,100,218,1596,2116,94,5170,784,298,118,3844,

%U 3900,334,484,358,2116,25900,146,394,7056,3322,530,8836,3844,466,1036,23690,42108,13924,206,538,10000,554,218,240548

%N Product_{d|n, A007431(d) > 0} prime(A007431(d)), where A007431 is the Möbius transform of Euler's totient function.

%H Antti Karttunen, <a href="/A318838/b318838.txt">Table of n, a(n) for n = 1..4096</a>

%F a(n) = product_{d|n} A008578(1+A007431(d)).

%F For all n >= 1, A056239(a(n)) = A000010(n).

%o (PARI)

%o A007431(n) = sumdiv(n,d,moebius(n/d)*eulerphi(d));

%o A318838(n) = { my(m=1); fordiv(n,d,if((A007431(d)!=0),m *= prime(A007431(d)))); (m); };

%Y Cf. A000010, A007431, A318839 (rgs-transform).

%Y Cf. also A318836.

%K nonn

%O 1,1

%A _Antti Karttunen_, Sep 05 2018