A323511 a(n) = 1 if n is squarefree and A079559(A048675(n)) = 1, otherwise a(n) = 0.

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, 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, 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

Offset

1

Comments

Characteristic function for A283476, the range of A283475.

Links

Antti Karttunen, Table of n, a(n) for n = 1..65537

Index entries for characteristic functions

Formula

a(n) = A008966(n) * A079559(A048675(n)).

a(A283475(n)) = a(A283476(n)) = 1.

Prog

(PARI)
A048675(n) = my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; \\ From A048675
A036987(n) = !bitand(n, 1+n);
A053644(n) = { my(k=1); while(k<=n, k<<=1); (k>>1); }; \\ From A053644
A053645(n) = (n-A053644(n));
A079559(n) = if(!n, 1, (1-A036987(1+n))*A079559(A053645(1+n)));
A323511(n) = (issquarefree(n) && A079559(A048675(n)));

Crossrefs

Cf. A008966, A048675, A079559, A283475, A283476.

Cf. also A323512.

Sequence in context: A267683 A267869 A068434 * A284391 A127015 A281114

Adjacent sequences: A323508 A323509 A323510 * A323512 A323513 A323514

Keyword

nonn

Author

Antti Karttunen, Jan 18 2019

Status

approved