OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..200
Project Euler, Problem 427: n-sequences
FORMULA
a(n) = Sum_{k=1..n} k*A228154(n,k).
a(n) ~ (2-exp(-1)) * n^n. - Vaclav Kotesovec, Sep 10 2014
EXAMPLE
a(2) = 6 = 2+1+1+2: [1,1], [1,2], [2,1], [2,2].
MAPLE
a:= proc(n) option remember; local b; b:=
proc(m, s, i) option remember; `if`(m>i or s>m, 0,
`if`(i=1, n, `if`(s=1, (n-1)*add(b(m, h, i-1), h=1..m),
b(m, s-1, i-1) +`if`(s=m, b(m-1, s-1, i-1), 0))))
end; forget(b);
add(m*add(b(m, s, n), s=1..m), m=1..n)
end:
seq(a(n), n=0..30);
MATHEMATICA
a[n_] := a[n] = Module[{b}, b[m_, s_, i_] := b[m, s, i] = If[m>i || s>m, 0, If[i==1, n, If[s==1, (n-1) Sum[b[m, h, i-1], {h, 1, m}], b[m, s-1, i-1] + If[s==m, b[m-1, s-1, i-1], 0]]]]; Sum[m Sum[b[m, s, n], {s, 1, m}], {m, 1, n}]];
Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Dec 20 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Aug 17 2013
STATUS
approved