OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
B. A. Huberman and T. Hogg, Complexity and adaptation, Evolution, games and learning (Los Alamos, N.M., 1985). Phys. D 22 (1986), no. 1-3, 376-384.
N. J. A. Sloane, Transforms
FORMULA
Euler transform applied 4 times to all-1's sequence.
MAPLE
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: b[0]:= etr(1): for k from 1 to 2 do b[k]:= etr(b[k-1]) od: a:= etr(b[2]): seq(a(n), n=0..25); # Alois P. Heinz, Sep 08 2008
MATHEMATICA
i[ n_, m_ ] := 1 /; m==1 || n==0; i[ n_, m_ ] := (i[ n, m ]=1/n Sum[ i[ k, m ] Plus @@ ((# i[ #, m-1 ])& /@ Divisors[ n-k ]), {k, 0, n-1} ]) /; n>0 && m>1
(* Second program: *)
A[0|1, _] = A[_, 1] = 1; A[n_, k_] := A[n, k] = Sum[DivisorSum[j, A[#, k-1] * #&]*A[n-j, k], {j, 1, n}]/n;
a[n_] := A[n, 5];
Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Jun 16 2018, after A290353 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
EXTENSIONS
More terms from Christian G. Bower, Aug 15 1998
STATUS
approved