seq(sum(stirling2(n, 3+4*k), k=0..(n-3)/4), n=3..26); #
Emeric Deutsch, Dec 15 2004
# Alternative:
with(combinat):
b:= proc(n, i, m) option remember; `if`(n=0, `if`(m=3, 1, 0),
`if`(i<1, 0, add(multinomial(n, n-i*j, i$j)/j!*
b(n-i*j, i-1, irem(m+j, 4)), j=0..n/i)))
end:
a:= n-> b(n$2, 0):