OFFSET
1,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..250
Wikipedia, Arithmetic derivative
MAPLE
d:= n-> n*add(i[2]/i[1], i=ifactors(n)[2]):
q:= n-> not isprime(n) and irem(d(d(n)), n)=0:
select(q, [$1..30000])[]; # Alois P. Heinz, Oct 15 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 15 2021
STATUS
approved