OFFSET
1,2
COMMENTS
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
FORMULA
a(n) ~ sqrt(3) * zeta(3) * exp(Pi*sqrt(2*n/3)) / Pi^2. - Vaclav Kotesovec, Oct 20 2024
MAPLE
b:= proc(n, i) option remember; `if`(n=0 or i=1, [1, n],
b(n, i-1)+(p-> p+[0, p[1]*i^2])(b(n-i, min(n-i, i))))
end:
a:= n-> (b(n$2)-b(n-1$2))[2]:
seq(a(n), n=1..40); # Alois P. Heinz, Feb 23 2022
MATHEMATICA
b[n_, i_] := b[n, i] = If[n == 0 || i == 1, {1, n},
b[n, i-1] + Function[p, p + {0, p[[1]]*i^2}][b[n-i, Min[n-i, i]]]];
a[n_] := (b[n, n] - b[n-1, n-1])[[2]];
Table[a[n], {n, 1, 40}] (* Jean-François Alcover, Apr 25 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Omar E. Pol, Feb 08 2012
STATUS
approved