OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
EXAMPLE
a(4) = 10: [2,2,2,2,1,1], [2,2,2,2,2], [3,2,2,2,1], [3,3,2,2].
MAPLE
b:= proc(n, i, t) option remember; `if`(n=0, `if`(t=0, 1, 0),
`if`(i<1, 0, b(n, i-1, t) +`if`(i>n, 0,
b(n-i, i, t -`if`(t>0 and isprime(i), 1, 0)))))
end:
a:= proc(n) option remember; local k;
for k from a(n-1) while b(k, k, n)<n do od; k
end: a(0):=0:
seq(a(n), n=0..50); # Alois P. Heinz, Jan 18 2014
MATHEMATICA
$RecursionLimit = 1000; b[n_, i_, t_] := b[n, i, t] = If[n == 0, If[t == 0, 1, 0], If[i<1, 0, b[n, i-1, t] + If[i>n, 0, b[n-i, i, t - If[t>0 && PrimeQ[i], 1, 0]]]]]; a[n_] := a[n] = Module[{k}, For[k = a[n-1], b[k, k, n] < n, k++]; k]; a[0] = 0; Table[a[n], {n, 0, 61}] (* Jean-François Alcover, Jan 27 2014, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
J. Stauduhar, Jan 18 2014
EXTENSIONS
More terms from Alois P. Heinz, Jan 18 2014
STATUS
approved