OFFSET
1,1
EXAMPLE
The first eight arithmetic derivatives of 75 are 55, 16, 32, 80, 176, 368, 752, 1520 and d(55) + d(16) + d(32) + d(80) + d(176) + d(368) + d(752) + d(1520) = 4 + 5 + 6 + 10 + 10 + 10 + 10 + 20 = 75.
MAPLE
with(numtheory): P:=proc(q) local a, b, k, n, p; for n from 1 to q do a:=0; k:=1; b:=n;
while a<n do b:=b*add(op(2, p)/op(1, p), p=ifactors(b)[2]); if b>0 then a:=a+tau(b); else break; fi; od;
if n=a then print(n); fi; od; end: P(10^6);
PROG
(PARI) ad(n) = if (n<1, 0, my(f = factor(n)); n*sum(k=1, #f~, f[k, 2]/f[k, 1]));
isok(n) = {ss = 0; kn = n; while (ss < n, der = ad(kn); if (der == 0, break); ss += numdiv(der); kn = der); ss == n; } \\ Michel Marcus, Apr 08 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Apr 06 2016
EXTENSIONS
a(35)-a(45) from Amiram Eldar, Jun 23 2023
STATUS
approved