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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A144036 Shifts left when Euler transform applied 4 times. 2
0, 1, 1, 5, 19, 89, 410, 2052, 10440, 54874, 293549, 1597621, 8807766, 49107289, 276358791, 1567866228, 8957204966, 51486464912, 297548288251, 1727856600935, 10076859047404, 58996263573440, 346614270372761, 2042929868812385, 12076076910981403 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

LINKS

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

N. J. A. Sloane, Transforms

MAPLE

k:=4: with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:='a': b[1]:=etr(a): for t from 2 to k do b[t]:= etr(b[t-1]) od: a:= n-> `if`(n<2, n, b[k](n-1)): seq(a(n), n=0..30);

MATHEMATICA

etr[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[ j]}]*b[n-j], {j, 1, n}]/n]; b]; A[n_, k_] := Module[{a, b, t}, b[1] = etr[a]; For[ t = 2, t <= k, t++, b[t] = etr[b[t-1]]]; a = Function[m, If[m == 1, 1, b[k][m-1]]]; a[n]]; Table[A[n, 4], {n, 0, 30} ] (* Jean-Fran├žois Alcover, Mar 05 2015, after Alois P. Heinz *)

CROSSREFS

4th column of A144042.

Sequence in context: A149802 A149803 A149804 * A110210 A244899 A147139

Adjacent sequences:  A144033 A144034 A144035 * A144037 A144038 A144039

KEYWORD

eigen,nonn

AUTHOR

Alois P. Heinz, Sep 07 2008

STATUS

approved

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

Content is available under The OEIS End-User License Agreement .

Last modified May 27 11:17 EDT 2015. Contains 257872 sequences.