|
%I #6 Mar 19 2018 22:14:40
%S 1,1,1,2,1,1,1,6,2,1,1,12,1,1,1,30,1,6,1,20,1,1,1,540,2,1,12,60,1,1,1,
%T 210,1,1,1,2520,1,1,1,1260,1,1,1,84,20,1,1,94500,2,6,1,140,1,540,1,
%U 18900,1,1,1,25200,1,1,60,2310,1,1,1,44,1,1,1,8731800,1,1,12,132,1,1,1,346500,168,1,1,39600,1,1,1,41580,1,1260
%N a(n) = Product_{d|n} A019565(d)^(1-A008966(n/d)).
%H Antti Karttunen, <a href="/A300830/b300830.txt">Table of n, a(n) for n = 1..8192</a>
%F a(n) = Product_{d|n} A019565(d)^(1-abs(A008683(n/d))).
%F a(n) = A293214(n) / (A300831(n)*A300832(n)).
%o (PARI)
%o A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
%o A300830(n) = { my(m=1); fordiv(n,d,if(!moebius(n/d),m *= A019565(d))); m; };
%Y Cf. A008683, A008966, A019565.
%Y Cf. also A293214, A300831, A300832, A300833.
%K nonn
%O 1,4
%A _Antti Karttunen_, Mar 16 2018
|
