OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..5000
FORMULA
a(n) == n (mod 2).
EXAMPLE
a(6) = 18 = 1 + 1 + 1 + 1 + 1 + 1 + 2 + 2 + 1 + 1 + 6: (1)23, (1)32, 2(1)3, 23(1), 3(1)2, 32(1), (2)4, 4(2), (1)5, 5(1), (6).
MAPLE
b:= proc(n, i, p) option remember; `if`(i*(i+1)/2<n or i<1, 0,
`if`(i=n, i*p!, b(n-i, min(n-i, i-1), p+1))+b(n, i-1, p))
end:
a:= n-> b(n$2, 1):
seq(a(n), n=0..50);
MATHEMATICA
b[n_, i_, p_] := b[n, i, p] = If[i(i+1)/2 < n || i < 1, 0,
If[i == n, i*p!, b[n-i, Min[n-i, i-1], p+1]] + b[n, i-1, p]];
a[n_] := b[n, n, 1];
Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Jul 12 2021, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Aug 07 2020
STATUS
approved