login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A258652 Sum of the k-th arithmetic derivative of n-k for k=0..n. 2

%I #12 Jun 01 2018 04:04:42

%S 0,1,2,4,5,9,11,16,14,25,36,59,99,209,419,860,1730,3862,9464,21868,

%T 74371,244648,727345,3098351,13469007,56269849,281642632,1406177909,

%U 9597415332,58891421656,411673964638,3406742649805,24202753250241,176482943622608

%N Sum of the k-th arithmetic derivative of n-k for k=0..n.

%H Alois P. Heinz, <a href="/A258652/b258652.txt">Table of n, a(n) for n = 0..100</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Arithmetic_derivative">Arithmetic derivative</a>

%F a(n) = Sum_{k=0..n} A258651(n-k,k).

%p d:= n-> n*add(i[2]/i[1], i=ifactors(n)[2]):

%p A:= proc(n, k) option remember; `if`(k=0, n, d(A(n, k-1))) end:

%p a:= proc(n) option remember; add(A(h, n-h), h=0..n) end:

%p seq(a(n), n=0..40);

%t d[n_ /; n>1] := n*Sum[i[[2]]/i[[1]], {i, FactorInteger[n]}]; d[_] = 0;

%t A[n_, k_] := A[n, k] = If[k == 0, n, d[A[n, k-1]]];

%t a[n_] := a[n] = Sum[A[h, n-h], {h, 0, n}];

%t Table[a[n], {n, 0, 40}] (* _Jean-François Alcover_, Jun 01 2018, from Maple *)

%Y Antidiagonal sums of A258651.

%Y Cf. A003415.

%K nonn

%O 0,3

%A _Alois P. Heinz_, Jun 06 2015

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 14:13 EDT 2024. Contains 371254 sequences. (Running on oeis4.)