OFFSET
0,4
LINKS
EXAMPLE
a(6) = 10 because we have [5, 1], [4, 2], [3, 2, 1], [3, 1, 2], [2, 4], [2, 3, 1], [2, 1, 3], [1, 5], [1, 3, 2] and [1, 2, 3].
MAPLE
N:= 50: # for a(0)..a(N)
P:= select(isprime, [2, seq(i, i=3..N, 2)]):
PP:= sort([1, seq(seq(p^j, j = 1 .. ilog[p](N)), p=P)]):G:= 1:
for s in PP do
G:= G + series(G*x*y^s, y, N+1);
od:
G:= convert(G, polynom):
T:= add(coeff(G, x, i)*i!, i=0..N):
seq(coeff(T, y, i), i=0..N); # Robert Israel, Jun 28 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 01 2020
STATUS
approved