A371084 a(n) = 1 if n is a nonsquarefree number not divisible by p^p for any prime p, otherwise 0.

0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1

Offset

1

Links

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

Index entries for characteristic functions

Formula

a(n) = (1-A008966(n)) * A359550(n).

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} (1 - 1/p^p) - 6/Pi^2 = 0.1140631325... . - Amiram Eldar, Mar 12 2024

Prog

(PARI)
A359550(n) = { my(f = factor(n)); prod(k=1, #f~, (f[k, 2]<f[k, 1])); };
A371084(n) = (!issquarefree(n) && A359550(n));

Crossrefs

Characteristic function of A371083.

Cf. A008966, A059956, A359550.

Sequence in context: A353477 A044939 A373973 * A353626 A347244 A366289

Adjacent sequences: A371081 A371082 A371083 * A371085 A371086 A371087

Keyword

nonn

Author

Antti Karttunen, Mar 12 2024

Status

approved