OFFSET
1,1
EXAMPLE
Arithmetic derivative of 2401 is 1372. Prime factors of 2401 are 7^4; prime factors of 1372 are 2^2, 7^3 and 1/7 + 1/7 + 1/7 + 1/7 + 1/2 + 1/2 + 1/7 + 1/7 + 1/7 = 2.
MAPLE
with(numtheory); P:= proc(q) local a, b, c, n, p;
for n from 2 to q do a:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
b:=ifactors(a)[2]; c:=ifactors(n)[2]; if type(add(c[k][2]/c[k][1], k=1..nops(c))+add(b[k][2]/b[k][1], k=1..nops(b)), integer) then print(n); fi; od; end: P(10^9);
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Paolo P. Lava, Mar 20 2014
EXTENSIONS
a(15)-a(21) from Giovanni Resta, Mar 20 2014
STATUS
approved