login
A368913
a(n) = 1 if there is no prime p such that p^p divides A342001(n), but for A003415(n) such a prime exists, otherwise 0. Here A003415(n) is the arithmetic derivative of n, and A342001(n) = A003415(n) / A003557(n).
4
0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0
OFFSET
1
COMMENTS
Question: What is the asymptotic mean of this sequence?
FORMULA
a(n) = A368914(n) - A368915(n).
For all n >= 0, a(A048103(n)) = a(A276086(n)) = 0.
PROG
(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A003557(n) = (n/factorback(factorint(n)[, 1]));
A342001(n) = (A003415(n) / A003557(n));
A359550(n) = { my(f = factor(n)); prod(k=1, #f~, (f[k, 2]<f[k, 1])); };
A368914(n) = ((n>1)&&A359550(A342001(n)));
A368915(n) = ((n>1)&&A359550(A003415(n)));
A368913(n) = (A368914(n)-A368915(n));
CROSSREFS
Characteristic function of A368903.
Sequence in context: A353470 A342753 A358752 * A354820 A188086 A105563
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 09 2024
STATUS
approved