login
A391844
a(n) = 1 if the least prime not dividing the arithmetic derivative of n is equal to the least prime dividing n, otherwise 0. a(1) = 1 by convention.
7
1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1
OFFSET
1
PROG
(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A020639(n) = if(1==n, n, vecmin(factor(n)[, 1]));
A053669(n) = forprime(p=2, , if(n%p, return(p))); \\ From A053669
A391844(n) = if(n<2, n, (A053669(A003415(n))==A020639(n)));
CROSSREFS
Characteristic function of A391845.
Cf. also A370124.
Sequence in context: A379117 A185059 A179776 * A360130 A353570 A089802
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 16 2026
STATUS
approved