login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A258851 The pi-based arithmetic derivative of n: a(p) = pi(p) for p prime, a(u*v) = a(u)*v + u*a(v), where pi = A000720. 40
0, 0, 1, 2, 4, 3, 7, 4, 12, 12, 11, 5, 20, 6, 15, 19, 32, 7, 33, 8, 32, 26, 21, 9, 52, 30, 25, 54, 44, 10, 53, 11, 80, 37, 31, 41, 84, 12, 35, 44, 84, 13, 73, 14, 64, 87, 41, 15, 128, 56, 85, 55, 76, 16, 135, 58, 116, 62, 49, 17, 136, 18, 53, 120, 192, 69, 107 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

COMMENTS

The pi-based variant of the arithmetic derivative of n (A003415).

LINKS

Alois P. Heinz, Table of n, a(n) for n = 0..10000

Wikipedia, Arithmetic derivative

Index entries for sequences computed from indices in prime factorization

FORMULA

a(n) = n * Sum e_j*pi(p_j)/p_j for n = Product p_j^e_j with pi = A000720.

a(A258861(n)) = n; A258861 = pi-based antiderivative of n.

a(a(A258862(n))) = n; A258862 = second pi-based antiderivative of n.

a(a(a(A258995(n)))) = n; A258995 = third pi-based antiderivative of n.

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

(Scheme) (define (A258851 n) (if (<= n 1) 0 (+ (* (A055396 n) (A032742 n)) (* (A020639 n) (A258851 (A032742 n)))))) ;; Antti Karttunen, Mar 07 2017

CROSSREFS

Column k=1 of A258850, A258997.

First differences give A258863.

Partial sums give A258864.

Cf. A000720 (pi), A003415, A020639, A032742, A055396, A056239, A258861, A258862, A258995, A278510.

Sequence in context: A007203 A110412 A266285 * A194031 A064357 A191735

Adjacent sequences:  A258848 A258849 A258850 * A258852 A258853 A258854

KEYWORD

nonn,look

AUTHOR

Alois P. Heinz, Jun 12 2015

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 15 11:13 EST 2019. Contains 329144 sequences. (Running on oeis4.)