|
%I #6 Mar 19 2018 22:14:47
%S 1,1,1,1,1,2,1,1,1,2,1,3,1,2,2,1,1,6,1,3,2,2,1,5,1,2,1,3,1,180,1,1,2,
%T 2,2,15,1,2,2,5,1,540,1,3,6,2,1,7,1,10,2,3,1,14,2,5,2,2,1,1575,1,2,6,
%U 1,2,756,1,3,2,900,1,35,1,2,10,3,2,1260,1,7,1,2,1,7875,2,2,2,5,1,44100,2,3,2,2,2,11,1,30,6,21,1,396,1,5,1800
%N a(n) = Product_{d|n, d<n} A019565(d)^[moebius(n/d) = +1].
%H Antti Karttunen, <a href="/A300831/b300831.txt">Table of n, a(n) for n = 1..8192</a>
%F a(n) = A293214(n) / (A300830(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 A300831(n) = { my(m=1); fordiv(n,d,if((d < n)&&(1==moebius(n/d)), m *= A019565(d))); m; };
%Y Cf. A008683, A008966, A019565.
%Y Cf. also A293214, A300830, A300832, A300833.
%K nonn
%O 1,6
%A _Antti Karttunen_, Mar 16 2018
|
