OFFSET
1,1
EXAMPLE
Arithmetic derivative of 1728 is 6912. Prime factors of 1728 are 2^6, 3^3; prime factors of 6912 are 2^8, 3^3 and 6*1/2 + 3*1/3 + 8*1/2 + 3*1/3 = 9.
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(12)-a(15) from Giovanni Resta, Mar 20 2014
STATUS
approved