|
%I #7 Jan 18 2019 20:54:37
%S 1,1,0,0,1,1,1,0,0,0,1,0,1,0,0,0,1,0,1,0,1,0,1,0,0,0,0,0,1,1,1,0,1,0,
%T 0,0,1,0,1,0,1,1,1,0,0,0,1,0,0,0,1,0,1,0,0,0,1,0,1,0,1,0,0,0,0,1,1,0,
%U 1,0,1,0,1,0,0,0,0,1,1,0,0,0,1,0,0,0,1,0,1,0,0,0,1,0,0,0,1,0,0,0,1,1,1,0,0,0,1,0,1,0,1,0,1,1,0
%N a(n) = 1 if n is squarefree and A079559(A048675(n)) = 1, otherwise a(n) = 0.
%C Characteristic function for A283476, the range of A283475.
%H Antti Karttunen, <a href="/A323511/b323511.txt">Table of n, a(n) for n = 1..65537</a>
%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>
%F a(n) = A008966(n) * A079559(A048675(n)).
%F a(A283475(n)) = a(A283476(n)) = 1.
%o (PARI)
%o A048675(n) = my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; \\ From A048675
%o A036987(n) = !bitand(n,1+n);
%o A053644(n) = { my(k=1); while(k<=n, k<<=1); (k>>1); }; \\ From A053644
%o A053645(n) = (n-A053644(n));
%o A079559(n) = if(!n,1,(1-A036987(1+n))*A079559(A053645(1+n)));
%o A323511(n) = (issquarefree(n) && A079559(A048675(n)));
%Y Cf. A008966, A048675, A079559, A283475, A283476.
%Y Cf. also A323512.
%K nonn
%O 1
%A _Antti Karttunen_, Jan 18 2019
|
