|
%I #12 May 15 2022 23:49:19
%S 0,1,1,0,1,0,1,1,0,0,1,0,1,0,0,0,1,0,1,0,0,0,1,0,0,0,1,0,1,0,1,1,0,0,
%T 0,0,1,0,0,0,1,0,1,0,0,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,1,0,0,0,0,0,1,0,
%U 0,0,1,0,1,0,0,0,0,0,1,0,0,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,1,0,0,0,1,0,1,0,0,0,1,0,1,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,1
%N a(n) = 1 if n is a prime power with an odd exponent, otherwise 0.
%H Antti Karttunen, <a href="/A353817/b353817.txt">Table of n, a(n) for n = 1..65539</a>
%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>
%H <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>
%F a(n) = A010055(n) - A353800(n).
%F a(n) >= A010051(n).
%F a(n) >= A354028(n).
%t a[n_] := If[n > 1 && Length[(f = FactorInteger[n])] == 1 && OddQ[f[[1, 2]]], 1, 0]; Array[a, 100] (* _Amiram Eldar_, May 15 2022 *)
%o (PARI) A353817(n) = (isprimepower(n)%2);
%Y Characteristic function of A246551.
%Y Cf. A010051, A010055, A353800, A354028.
%K nonn
%O 1
%A _Antti Karttunen_, May 15 2022
|
