OFFSET
0,4
COMMENTS
The pi-based variant of the arithmetic derivative of n (A003415).
LINKS
FORMULA
a(n) = n * Sum e_j*pi(p_j)/p_j for n = Product p_j^e_j with pi = A000720.
a(0) = a(0*p) = a(0)*p + 0*a(p) = a(0)*p for all p => a(0) = 0.
a(p) = a(1*p) = a(1)*p + 1*a(p) = a(1)*p + a(p) for all p => a(1) = 0.
a(u^v) = v * u^(v-1) * a(u).
MAPLE
with(numtheory):
a:= n-> n*add(i[2]*pi(i[1])/i[1], i=ifactors(n)[2]):
seq(a(n), n=0..100);
MATHEMATICA
a[n_] := n*Total[Last[#]*PrimePi[First[#]]/First[#]& /@ FactorInteger[n]]; a[0] = 0; Array[a, 100, 0] (* Jean-François Alcover, Apr 24 2016 *)
PROG
(PARI) A258851(n)=n*sum(i=1, #n=factor(n)~, n[2, i]*primepi(n[1, i])/n[1, i]) \\ M. F. Hasler, Jul 13 2015
KEYWORD
nonn,look
AUTHOR
Alois P. Heinz, Jun 12 2015
STATUS
approved