OFFSET
1,2
COMMENTS
Here, m' denotes the arithmetic derivative of m (A003415).
Not the same as A211781 since this sequence does not contain 225, 252, etc.
All prime numbers p are in the sequence since p' = 1 | p.
Numbers k such that k' | k. - The original definition of the sequence.
Sequence consists of 0, primes, and the prime powers of the form p^p (A051674, that together with 0 give the only fixed points of A003415). This can be seen from theorems 4-6 given in the Ufnarovski & Ã…hlander paper. - Antti Karttunen, May 17 2025
LINKS
Antti Karttunen & R. J. Mathar, Table of n, a(n) for n = 1..20001
Victor Ufnarovski and Bo Ã…hlander, How to Differentiate a Number, J. Integer Seqs., Vol. 6, 2003, #03.3.4.
EXAMPLE
MAPLE
q:= n-> is(irem(n, n*add(i[2]/i[1], i=ifactors(n)[2]))=0):
select(q, [$2..300])[]; # Alois P. Heinz, Oct 11 2021
PROG
(PARI) ad(n) = vecsum([n/f[1]*f[2]|f<-factor(n+!n)~]); \\ A003415
isok(k) = (0==k) || ((k>1) && !(k % ad(k))); \\ Michel Marcus, Oct 10 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Oct 09 2021
EXTENSIONS
a(1) = 0 inserted because of a new, more inclusive definition. - Antti Karttunen, May 17 2025
STATUS
approved
