OFFSET
0,6
LINKS
EXAMPLE
a(16) = 4 because we have [13, 3], [3, 13], [11, 5] and [5, 11].
MAPLE
b:= proc(n, i, p) option remember; `if`(n=0, irem(1+p, 2)*p!, (s->
`if`(s>n, 0, b(n, i+1, p)+b(n-s, i+1, p+1)))(ithprime(i)))
end:
a:= n-> b(n, 1, 0):
seq(a(n), n=0..70); # Alois P. Heinz, Dec 04 2020
MATHEMATICA
b[n_, i_, p_] := b[n, i, p] = If[n == 0, Mod[1 + p, 2]*p!, Function[s, If[s > n, 0, b[n, i + 1, p] + b[n - s, i + 1, p + 1]]][Prime[i]]];
a[n_] := b[n, 1, 0];
Table[a[n], {n, 0, 70}] (* Jean-François Alcover, Feb 26 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 04 2020
STATUS
approved