|
%I #19 Jan 24 2024 17:08:03
%S 1,1,2,2,4,4,6,8,10,11,12,14,17,21,25,29,33,37,43,49,54,59,66,72,80,
%T 89,98,106,116,126,137,148,161,174,187,200,216,232,248,266,284,302,
%U 321,344,367,391,414,440,465,493,523,556,584,616,650,689,726,768,808
%N a(0) = a(1) = 1; a(n+1) is the number of partitions of n into parts a(k), 0 <= k <= n, each k occurring at most once.
%H Alois P. Heinz, <a href="/A057601/b057601.txt">Table of n, a(n) for n = 0..10000</a>
%e a(6) = 6 because 5 = a(0) + a(4) = a(0) + a(5) = a(1) + a(4) = a(1) + a(5) = a(0) + a(2) + a(3) = a(1) + a(2) + a(3).
%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<0, 0,
%p b(n, i-1)+`if`(a(i)>n, 0, b(n-a(i), i-1))))
%p end:
%p a:= proc(n) a(n):= `if`(n<2, 1, b(n-1, n-1)) end:
%p seq(a(n), n=0..70); # _Alois P. Heinz_, May 26 2013
%t b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 0, 0, b[n, i - 1] + If[a[i] > n, 0, b[n - a[i], i - 1]]]];
%t a[n_] := If[n < 2, 1, b[n - 1, n - 1]];
%t a /@ Range[0, 70] (* _Jean-François Alcover_, Nov 10 2020, after _Alois P. Heinz_ *)
%K nonn
%O 0,3
%A _Leroy Quet_, Oct 06 2000
%E More terms from _Naohiro Nomoto_, Oct 28 2001
%E Extended beyond a(45) by _Alois P. Heinz_, May 26 2013
|
