%I #13 Oct 31 2020 12:45:58
%S 1,1,0,0,0,0,0,0,1,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,2,0,0,0,0,0,
%T 0,2,6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,2,0,0,
%U 0,0,0,0,2,6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,6,0,0,0,0,0,0,6,24
%N Number of compositions (ordered partitions) of n into distinct cubes.
%H Alois P. Heinz, <a href="/A331845/b331845.txt">Table of n, a(n) for n = 0..65536</a>
%H <a href="/index/Com#comp">Index entries for sequences related to compositions</a>
%H <a href="/index/Su#ssq">Index entries for sequences related to sums of cubes</a>
%e a(36) = 6 because we have [27,8,1], [27,1,8], [8,27,1], [8,1,27], [1,27,8] and [1,8,27].
%p b:= proc(n, i, p) option remember;
%p `if`((i*(i+1)/2)^2<n, 0, `if`(n=0, p!,
%p `if`(i^3>n, 0, b(n-i^3, i-1, p+1))+b(n, i-1, p)))
%p end:
%p a:= n-> b(n, iroot(n, 3), 0):
%p seq(a(n), n=0..100); # _Alois P. Heinz_, Jan 30 2020
%t b[n_, i_, p_] := b[n, i, p] = If[(i(i+1)/2)^2 < n, 0, If[n == 0, p!, If[i^3 > n, 0, b[n-i^3, i-1, p+1]] + b[n, i-1, p]]];
%t a[n_] := b[n, Floor[n^(1/3)], 0];
%t a /@ Range[0, 100] (* _Jean-François Alcover_, Oct 31 2020, after _Alois P. Heinz_ *)
%Y Cf. A000578, A023358, A032020, A032021, A032022, A218396, A219107, A279329, A331843, A331844, A331846, A331847.
%K nonn,look
%O 0,10
%A _Ilya Gutkovskiy_, Jan 29 2020